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Summary 


Gamma-ray bursts (GRB) are extreme events taking place at cosmological distances. Their origin 
and mechanism have been puzzling for decades. They are crudely classified into two groups based 
on their duration, namely the short bursts and the long bursts. 

Such a classification has proven to be extremely useful to determine their possible progenitors: 
the merger of two compact objects for short bursts and the explosion of a (very) massive star for 
long bursts. Further classifying the long GRBs might give tighter constraints on their progenitor 
(initial mass, angular momentum, evolution stage at collapse) and on the emission mechanism(s). 

The understanding of several aspects of GRBs has greatly advanced after the launch of the 
Swift satellite, as it allows for multi-wavelength observation of both the prompt phase and the 
afterglow of GRBs. On the other hand, a world collaboration to point ground optical and radio 
telescopes has allowed many breakthroughs in the physics of GRBs, for instance with several 
detections of supernova in the late afterglow phase. 

In my thesis, I present evidence for the existence of a sub-class of long GRBs, based on their 
faint afterglow emission. These bursts were named low-luminosity afterglow (LLA) GRBs. I 
discuss the data analysis and the selection method of these bursts. Then, their main properties 
are described (prompt and afterglow). Their link to supernova is strong as 64% of all the bursts 
firmly associated to SNe are LLA GRBs. This motivated the study of supernovae in my thesis. 

Finally, I present additional properties of LLA GRBs: the study of their rate density, which 
seems to indicate a new distinct third class of events, the properties of their host galaxies, which 
show that they take place in young star-forming galaxies, not different from those of normal long 
GRBs. 

Additionally, I show that it is difficult to reconcile all differences between normal long GRBs 
and LLA GRBs only by considering instrumental or environmental effects, a different ejecta 
content or a different geometry for the burst. Thus, I conclude that LLA GRBs and normal long 
GRBs should have different properties. 

In a very rudimentary discussion of the possible progenitor, I indicate that a binary system is 
favored in the case of LLA GRB. The argument is based on the initial mass function of massive 
stars, on the larger rate density of LLA GRBs compared to the rate of normal long GRBs and 
on the type of accompanying SNe. 

Such a classification of GRBs is important to constrain their emission mechanisms and pos¬ 
sible progenitors, which are still highly debated. However, more multi-wavelength observations 
of weak bursts at small redshift are required to give tighter constraints on the properties of both 
the burst and its accompanying supernova if present. 
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Introduction 


1.1 Diversity of GRBs Progenitors 

Gamma-ray bursts (GRBs) are extraordinary events appearing randomly in the sky. They were 
discovered in 1967 [22], by the US Vela military satellites which were designed to monitor nuclear 
tests in the atmosphere and in the outer space after the Limited Nuclear Test Ban Treaty signed 
between the USSR and the USA. Within three years, 16 events were recorded by the satellites 
Vela 5 and 6 (see Figure O)- However, it was clear very quickly that these events originated 
from the sky and not from the Earth: since there were more than one satellite, it was possible 
to measure the time delay between two satellites to get the direction of the events. 



Figure 1.1: The Vela 5b satellite before its launch. Image Credit: NASA 
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A GRB appears as a sudden burst of high energy photons (keV - GeV), lasting from a few 
tens of milliseconds up to several minutes. Thanks to the results of BATSE [T it is known 
that they are distributed isotropically in the sky (see Figure 11.101) which is an indication of 
their cosmological origin, later confirmed by redsliift measurements. Since they are located at 
cosmological distances [25J J2U , the recorded flux at the Earth implies that they are powerful 
emitters and release an enormous isotropic equivalent energy ranging from 10 48 to 10 54 ergs. 
The variability observed in the radio emission (hours to weeks after the prompt emission) shows 
that the emission region has a stellar size in the order of 10 17 cm [2S], see Figure 9 of |2B]. When 
these two pieces of information are combined, they point towards a phenomenon running extreme 
physics near a compact object (stellar mass black hole or young magnetar). 

For a very long time, the amount of energy needed for this phenomenon puzzled the theoreti¬ 
cians, and several models were built without considering it. The models were operating a central 
(compact) engine and the questions were about how to provide the energy. Nowadays, several 
objects are supposed to be able to trigger a gamma-ray burst. 

Massive stars (for example Wolf-Rayet stars) are the first possibility. Indeed, they can collapse 
to form a black hole. The remaining of the star is accreted by a newly born black hole which 
powers the GRB [27] (see Figure fL2l) . 



Figure 1.2: Artist close view of a GRB, with the remaining of the progenitor being accreted by 
the central engine. Image Credit: NASA. 


Another possible progenitor is a binary system formed of two compact objects (Neutron 
Star (NS)-NS, Black Hole (BH)-NS, BH-BH), merging together as a result of orbital angular 
momentum lost when gravitational wave radiation is emitted 128] . 

The last possibility is a Magnetar which is a neutron star with a large magnetic field and a 
high rotation speed. The rotational energy is extracted by the rotating magnetic field, slowing 
down the neutron star. As a result it can eventually collapse into a black hole m- 

All these objects are formed by the death of stars: GRBs seem to be linked to the final act 
of the stellar evolution. This idea is even strengthened by the observation of a GRB-Supernova 
(SN) association. Other important observations link GRBs and massive stars in the place of 
birth of long GRBs. 
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1.2 The Last Stages of Stellar Evolution 

1.2.1 How do stars die? 

The first stars were created 13.8 billion years ago, shortly after the Big Bang [501151]. The 
lifetime of the most massive star is only a few million years (for example, it is 3 millions years for 
a 60 Mq), while it can be up to trillions of years (10 12 ) for the least massive stars [32]. The star 
evolution towards its end depends on the physical parameters of star: initial mass, metallicity, 
mass loss rate, rotation speed, etc. Below, the fate of stars is presented as a function of increasing 
initial mass. 

Brown and red dwarfs represent extremely low massive (0.072 Mq < M < 0.8 Mq) stars, which 
are mainly observed in binary systems [33[ 131]. They are thought to be formed by a collapsing 
cloud of gas and dust like normal stars, and they begin to burn their hydrogen. However, in the 
case of brown dwarfs, their surface temperature and luminosity never stabilize since their mass 
is not big enough to maintain thermonuclear fusion. Thus, they cool down as they age. [35] . On 
the other hand, for red dwarfs, helium (He) is produced and constantly remixed with hydrogen 
(H) in the whole volume of the star, thus avoiding the creation of a He core. They develop slowly 
for a long time, expected to be much larger than the age of universe. Thus, advanced red dwarfs 
are not observed. 

The low mass stars between 0.8 and 8 Mq take a few billion years to burn out their hydrogen 
[36] . At the end of hydrogen burning, the size of the star decreases (because the radiation pressure 
decreases) and the central density and temperature increase. Under this condition, helium fusion 
begins (see [33(35] for reviews) converting it to oxygen and carbon and increasing the radius of 
the star: this is the red supergiant phase. In the next contraction stage the density becomes so 
high that the size of the star stabilizes due to the electron degeneracy pressure: electrons are 
Fermi particles, so two of them cannot be at the same energy level, leading to a pressure-like 
strength at very high density. During this stage the outer layers of the star are expelled by strong 
stellar winds and form a planetary nebula, composed of hot gas, ionized by ultraviolet radiation 
from the core of the star, see one example in Figure 11.31 At the end, a white dwarf is formed, 
which cools down to become a black dwarf. The life cycle of a low massive star is summarized in 
Figure fU4l al. However, this kind of compact object cannot create a GRB, because the accretion 
onto the white dwarf, when it is forming, is very small. Thus the accretion disk cannot power 
a GRB and a black hole cannot be created since white dwarf cannot reach the Chandrasekhar 
limit (1.4 Mq). 

The massive stars with initial mass larger than 8 Mq cannot reach such a stable level: they 
start the CNO cycle, and they can burn all heavier elements starting with carbon and oxygen, 
up until the iron which is the most stable element (see Figure IT751 for various burning stages). 
They acquire an onion structure with the lighter elements composing the outer layers. 

When the iron core reaches its Chandrasekhar mass, electron degeneracy pressure cannot 
balance the gravitational pull. The core collapses, and in the meantime, the iron is photo- 
disintegrated into free nucleons and alpha particles. Finally, the collapse stops as the central 
density reaches 10 14 g.cm -3 , which allows for neutron degeneracy pressure, which is similar to 
electron degeneracy pressure. At the end, a proto-neutron star is formed, of dimensions in the 
order of some tens of kilometers. The sudden deceleration launches a shock wave 20 km from 
the center which travels upstream through the in-falling matter left of the iron core, which is 
accreted onto the proto-neutron star after being shocked. The shock loses energy as it fights its 
way out through the inflow, but neutrinos emitted by the cooling of the proto-neutron star pushes 
it on. For a few seconds, the proto-neutron star emits an enormous amount of energy mostly 
in neutrino and gravitational waves. A tiny fraction of this energy (approximately 10 51 ergs) 
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Figure 1.3: The Ring Nebula in Lyra, of size about 1.3 light year (to be compared to 7.3 x 10 -8 
light year for the size of the Sun), is a giant shell of gas surrounding a central star. Image Credit: 
NASA/Hubble Heritage 


is sufficient to make the star explode as a supernova (SN) by the shock wave which propogates 
through the outer layer of the star at thousands of kilometers per second ( e.g. |T for a review). 
Below initial masses of ~ 25 M 0 , neutron stars are formed. Above that, black holes form, either 
in a delayed manner by fallback of the ejected matter or directly during the iron-core collapse 
(above ~ 40 M 0 ) as it can be seen in Figure ITT>1 pQ. The life cycle of a very massive star is 
summarized in Figure iLTT bh 

However, the end of extremely massive stars is different. For example, above 100 M 0 , stars 
suffer electron-positron pair instability after carbon burning. This begins as a pulsational in¬ 
stability of the helium cores of ~ 40 M 0 . As its mass increases, the pulsations become more 
violent, ejecting any remaining hydrogen envelope and an increasing fraction of the helium core 
itself. An iron core can still form in hydrostatic equilibrium in such stars, but it collapses to a 
black hole. The pair-instability supernovae can produce all elements from oxygen to iron through 
nickel. Depending on the initial mass function, the creation of the Ni 56 isotope might be possible. 
Indeed, it requires the synthesis of at least five solar mass of Nickel (Ni 58 ), which is produced by 
the silicon burning process pQ. All of these massive stars (8 M 0 < M < 260 M 0 ) are thought 
to produce a cataclysmic event at the end of their life, which can be classified according to 
their observed properties. In brief, deaths of non-rotating massive stars is summarized in Table 
11.11 PQ. These discussions are summarized on Figure ITiil [i^ which also displays the type of the 
corresponding expected SNe. 

All these SNe will leave a remnant: a nebula composed of the outer layers of the massive star 
expelled during the collapse, with a compact object (black hole or neutron star) at its center. 
For instance, one of the first observed and recorded supernova was SN 1054 from the year of 
its explosion. In its remnant, called the Crab nebula, there is a young neutron star: the Crab 
Pulsar, see Figure [FT] 
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Nebula phase 
Red supergiant 
Planetary nebula 
White dwarf 
Black dwarf 



Red supergiant 
Supernova 
Neutron star 
Black hole 


Figure 1.4: The life cycle of a low mass star (a) and a high mass star (b). Both scenarios start 
from the nebula phase at the center of the figure. Depending on the initial mass, two evolution 
paths can be followed, which end by a black dwarf for small initial mass (towards the left) and 
a neutron star or a black hole for large initial mass (towards the right). Image Credit: NASA 


Table 1.1: Deaths of non rotating massive helium stars [I 


Main Sequence Mass 


He Core 


Supernovae Mechanism 


8 M 0 < M < 95 M 0 2.2 Mq < M < 40 M 0 


95 Mq < M < 130 M q 
130 M q < M < 260 M q 


40 Mq < M < 60 Mq 
60 Mq< M < 137 Mq 


Fe core collapse to neutron 
star or a black hole 
Pulsational pair instablity 
followed by Fe core collapse 
Pair instability supernovae 


M > 260 Mq 


M > 137 Mq 


direct collapse 
to a black hole 
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Figure 1.5: Illustratration of the internal shell structure of a supergiant on its last days as an 
example of a massive star which has synthesized elements up to iron. Image Credit: CSIRO 


1.2.2 Supernovae (SNe) classification 

Supernovae (SNe) are violent explosions associated with the death of stars. They are charac¬ 
terized by a large increase in brightness up to -20 magnitudes. Observationally, by identifying 
different elements in their spectra, they can be divided into different subclasses [3]. This clas¬ 
sification depends mainly on the composition of the star, so on its initial mass pfl. An early 
classification of SNe was proposed, in which type I and type II are differentiated by the absence 
(type I) or presence (type II) of hydrogen (H) lines [39]. Figure [L8] displays the spectra (left) 
and light-curves (right) of different types of SNe. 

Type II SNe are usually mentioned in the literature with four subclasses: type IIP (plateau in 
their light-curves), type IIL (linear decline in their light-curves) which constitutes the majority 
of type II SNe ( e.g. [40] 1. type Iln with narrow line in their spectra, and type lib which is 
an intermediate type of SNe, with early features of type II SNe which are replaced by type lb 
features at late times. 

Moreover, type I SNe are divided into three subclasses. Type la SNe show silicon lines in 
their spectra and are thought to originate from the thermonuclear explosion of an accreting 
white dwarf (e.g. 05)- Type lb SNe show helium lines but do not show silicon lines in their 
spectra. Finally, type Ic SNe are very similar to type lb SNe: they only lack He lines. The basic 
properties of SNe and their classification are shown in Figure 11.91 Type lb and Ic SNe have 
many similarities, therefore they are sometimes jointly called Type Ib/c supernovae. Together 
with type II SNe, they are referred to as core-collapse SNe, because they are believed to be 
produced by the core collapse of massive stars (e.g. 02]). However, type Ib/c and type II SNe 
have different progenitors (differently evolved stars). 

Type Ib/c SNe have lost their hydrogen (lb) and helium (Ic) envelopes before the explosion; 
thus they are called stripped-envelope SNe. The mass loss is thought to arise from the increased 
mass of the progenitor relative to the progenitors of type II SNe. This kind of progenitors show 
stronger stellar winds which can blow away the outer layers of the progenitor star. Moreover, the 
higher mass loss can also be explained by line-driven winds due to the increased metal content 
of the progenitor [121 01] 05]. Thus, the progenitors of type Ib/c SNe could also have higher 
metallicity than type II SNe and their most likely progenitors are believed to be massive W-R 
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Figure 1.6: Initial and final mass function of non-rotating stars of solar composition, from pQ. 
The blue curve shows the mass of the star without its envelope (if lost) while the red curve shows 
the final mass of the remnant. Heavy elements are produced in the dark and light green region. 
Below 8 Mq (initial mass), in the asymptotic giant branch (AGB), a star loses its envelope and 
leaves a CO (carbon-oxygen) or NeO (neon-oxygen) white dwarf. Above 8 Mq and below 25 Mq, 
the H envelope and part of the He core are ejected. Then, He is converted to C and O. Finally, 
heavier elements are synthesized, leading to the creation of a neutron star, as the star explodes 
as a type II SN. Above 25 Mq and below 33 Mq, a black hole forms either by fallback of the 
ejecta from the star or directly by the iron core collapse, leading to a type Ib/c SN explosion. 
Above 33 Mq, this explosion (type Ib/c) is the result of a Wolf-Rayet (W-R) star evolution. In 
this region, there are two possibilities. First, the low W-R mass loss rate, represented by the 
blue and red tick dashed lines. It can synthesize heavy elements, leading to different type of 
W-R stars: WNE (early nitrogen in its spectra), WC (carbon in its spectra), WO (oxygen in its 
spectra). The collapse of these W-R stars gives a faint SN. Second, there are high W-R mass 
loss rate, represented by the thin dashed lines (blue and red). In this region, W-R stars collapse 
and form a bright SN. The final remnant is a black hole. 


stars surrounded by a dense wind envelope [46]. Such winds are at the origin of broad emission 
lines in the spectrum of the stars and change the surface composition, reflecting the presence of 
heavy elements created by nuclear burning (see Figure ITTBTl . These kinds of stars are thought 
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Figure 1.7: The remnant SN 1054 (the Crab nebula in visual band and X-ray with clearly seen 
crab pulsar). Image Credit: NASA 
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Figure 1.8: Left: Comparison of the spectra of different SNe types: SN 1994D for type la, 
SN 1999em for type II, SN 2004gq for type lb, SN 2004gk for type Ic. The broad line Ic SN 
( e.g. SN 1998bw) does not show H or He lines in its optical spectra and its lines widths are 2-3 
times larger than the normal type Ic SN which is indicative of high velocity ejecta [2]. Right: 
The light-curve of different type of SNe [5]. 


to be at the origin of type I SNe because they have little H (as in WNL: nitrogen in late type 
spectra) which can even be completely absent, like in WC and WO. These N, C, and O subtypes 
of W-R stars indicate the presence of strong lines of nitrogen, carbon, or oxygen in their spectra 

SZl- 

In addition, most massive stars are in binary systems. Because of the gravitational interac¬ 
tions within the binary, H can be removed from one of the star, which will end as a type Ib/c 
SN. 
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Figure 1.9: Schematization of supernovae classification 


1.2.3 Extreme case : Gamma-Ray Burst (GRBs) 

Interestingly enough, all GRBs which are associated to SNe were identified to the stripped- 
envelope SNe Ib/c. If the released isotropic energy of GRBs is between 10 48 -10 54 ergs, the 
emission is often assumed to be collimated in jets in order to decrease the required energy. 
In the collapsar model, the GRB is assumed to be emitted by a massive star experiencing a 
gravitational collapse m- Therefore, very special conditions are required for a star to evolve all 
the way to a gamma-ray burst: it should be very massive (probably at least 40 M 0 on the main 
sequence) and rapidly rotating to form a black hole and an accretion disk of mass around few 
0.1 M 0 . GRBs can be powered from the black hole by at least two different ways: 

• neutrino annihilation: u e + v~ = e + + e~ |481 149j . In this process, neutrinos are efficiently 
produced in the accretion disk because of its high temperature and high density. These 
neutrinos escape the accretion disk and annihilate in the polar region of the black hole. As 
large amounts of neutrinos are produced, the created electron positron plasma is optically 
thick and it is at the origin of the fireball. 

• magnetically dominated instability (Blandford-Znajek mechanism) [501 . In this process, 
the magnetic field of the black hole and of the accretion disk extracts the rotational energy 
of the black hole. 

In other models, the merger of two compact objects (two neutron stars or a neutron star with 
a black hole) is involved [28(. In these models, the energy can be extracted thanks to neutrino- 
antineutrino annihilation produced in the hot and dense torus and this energy available for the 
relativistic jet is sufficient to power a short GRB if a modest beaming is assumed. 
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1.3 Observational Properties of GRBs 

1.3.1 Spatial properties 

After the discovery of GRBs and before the discovery of their afterglows, they were explained 
by more than 100 different theoretical models m, the reason being that their distance was not 
known. The first clue about the distance came with the result of BATSE: it was observed that 
GRBs are distributed isotropically on the sky (see Figure 11.101 for all observed burst by the 
BATSE satellite). It implies that GRBs originate at cosmological distances [22] or from a large 
halo around our galaxy. 



+90 


+180 


-180 



Fluence, 50-300 keV (ergs cm' 2 ) 


Figure 1.10: The distribution of all 2704 GRBs detected by BATSE satellite: they are clearly 
isotropically distributed [J]. Image Credit: NASA/BATSE 


This controversy was solved by the redsliift measurements of GRB 970228 at redshift z = 
0.835 (53] which was made possible by the BeppoSAX satellite. It provided the localization of 
the burst with less than an arcminute uncertainty in radius |23j . which allowed a follow-up at 
optical [24] and other wavelengths. 

1.3.2 Temporal properties 

Based on the observations, two main times of emission can be defined, namely the prompt phase 
and the afterglow. The prompt phase can be associated to the central engine activity; the 
afterglow is a long-lasting emission, gradually decreasing, coming after the prompt phase. It 
is thought to be the result of the interaction of a relativistically expanding plasma with the 
environment of the progenitor, with the possible contamination of emission from late central 
engine activity. 

The prompt phase is characterized by a high flux of gamma-ray photons (keV - MeV). In 
this phase, each source shows different behavior and trend in its light-curve: single peak, double 
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peaks, multiple peaks, smooth or spiky light-curves, see Figure 11.111 The observed variability 
is high, down to 10 miliseconds [53]. The duration of the prompt phase called Tgo corresponds 
to the time during which 90% of the energy is emitted in the keV - MeV range. In 1993, with 
the data collected by BATSE, It was found that the Tgo distribution of all sources shows two 
separated groups, see Figure [LT71 Specifically, the short bursts have a duration Tgo< 2 s with 
the average observed value around 0.73 seconds while the long bursts have a duration Tgo> 2 s 
with the average observed value around 17 seconds [55] (551 [57]. 



TrXigy T 4Q6 



Trigger 1974 







Figure 1.11: Different light-curves from the GRBs prompt phase, as observed by BATSE (5, O'. 


1.3.3 Spectral properties of the prompt emission 


The spectrum of the prompt phase is well-fitted by the Band model 
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where the four parameters are the amplitude A , the low and high energy spectral indexes, a and 
/3 respectively, and the spectral break energy Eq. This function is made of two broken power law 
smoothly joined. 

The mean values of the spectral indexes are a = - 0.92 ± 0.42 and [3 = - 2.27 ± 0.01 for the 
long GRBs and a = - 0.4 ± 0.5 and ,6 = - 2.25 [55] for the short GRBs [101157]. For some bursts, 
the low energy spectral index is larger than - 2/3, which is not compatible with the synchrotron 
emission theory (the so-called synchrotron line of death) pT| . 

There is also another difference between short and long GRB spectra. The peak energy of 
short GRBs is on average larger (harder spectrum) than the peak energy of the long ones (mean 
values 398 keV and 214 keV respectively [57]). 
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BATSE 4B Catalog 



Figure 1.12: Duration distribution of short-long GRBs observed by BATSE [7 . Image Credit: 
NASA/BATSE 


The spectra of the prompt phase have been fitted with other models EH H3], for instance 
band plus black body, to understand and separate the different emission mechanisms (thermal 
or non-thermal). 

Very high energy photons (GeV) have been detected in several GRBs. The spectrum is well- 
fitted by a power law and a break in energy can sometimes be seen. The main properties of that 
feature are that it can be delayed for some seconds compared to the prompt phase and it lasts 
much longer. The emission mechanisms producing these photons are still puzzling and strongly 
debated. 

1.3.4 Afterglow 

The afterglow is observed at all wavelengths: X-ray [23], optical [M], IR- and radio (2T • Thanks 
to its low variability and observed time range (from minutes to weeks after the GRB event), 
a canonical X-ray light-curve for the afterglow was defined from the result of Sm/t/BAT-XRT 
instruments. It is displayed on Figure 11.131 the 0 symbol indicates the prompt phase, and the 
four remaining segments, with their corresponding temporal indexes, are associated two by two 
and identified as early and late afterglow [MIESIES]: i and II (respectively the steep and shallow 
decay), and III and IV (respectively a standard afterglow and a jet break) [67]. Part I and III, 
marked by solid lines, are most common and the other three components, marked by dashed 
lines, are only observed in a fraction of all bursts. Part I is thought to be associated with the 
prompt phase 0[§E] when the central engine is still active; the rest of the afterglow is due to 
the dynamics of the interaction between the jet and the surrounding medium. 

The spectra in the keV range of all phases but the prompt are well-fitted with a simple 
power-law model: 



N(E) = A 


100 keV 


( 1 . 2 ) 
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where A is the amplitude and j3 is the spectral index. 

The spectral index is constant throughout parts II, III, IV and is in the order of /3 ~ 1. 
However, it is softer for parts I and V (flares which are thought to be the extended emission of 
the prompt phase [55]). 



Figure 1.13: The schema of the various afterglow phases in long GRBs shown in the log flux-time 
plane [8]. The prompt phase (0) is generally followed by a steep decline afterglow (I) which can 
then break to a shallower decline called plateau (II), a standard afterglow phase (III), comes 
after and possibly, a jet break (IV). Sometimes X-ray flares can be seen (V) [8]. 


Even if the detection rate is high [70j . the number of well-sampled light-curves is smaller in 
other wavelengths: ~300 in X-ray mm ~68 in optical m and in radio. These light-curves 
can be compared with the X-ray light-curves to see if they share common properties. The most 
important result from the optical photometric observations is the optical bump in the light-curve 
between 10-25 days after the GRB explosion, which is interpreted as an SN explosion and indicate 
the connection between GRBs and SNe. On the other hand, the optical spectral observations 
provide the redshift by the measurement of absorption and emission lines. 


1.3.5 GRB-SN association 

The first direct evidence for a GRB-SN association was made when GRB 980425 was spectro¬ 
scopically and photometrically linked with type Ic-broad line SN 1998bw [74], as seen in Figure 
PI This connection was also predicted theoretically by the collapsar model m- 

This connection was further confirmed in 2003 with the spectroscopic and photometric as¬ 
sociation between GRB 030329 and SN 2003dh [75] ES], and GRB 031203 and SN 20031w [77) . 
These three bursts provided clear evidence that the progenitors of some, if not all long GRBs, 
are associated to the explosion of a massive star. Since then, other SNe were observed, either 
spectroscopically or identified as an optical bump in the late optical light-curve around 10 days 
after the burst. 
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1.3.6 Different kinds of GRBs 

Short bursts vs. long bursts: The study of the duration of the prompt phase has shown 
the existence of two classes of GRBs, namely the short and the long bursts, see subsection ll.3.21 
They are also differentiated by the hardness ratio [E3IS2 defined as the ratio of the bolometric 
fluence measured in two different bands: Sioo-aookev/Sso-iookeV) see Figure H.14I The hardness 
ratio of short GRBs is slightly larger than that of long GRBs. 



Figure 1.14: Hardness ratio of short - long GRBs obtained by BATSE and BeppoSAX [7]. Image 
Credit: NASA/BATSE 


In addition, Swift has shown that the redshift distributions of short and long bursts are not 
consistent m, which strengthen the theory that at least two different progenitors are respon¬ 
sible for GRBs. Another strong hint comes from the associated host galaxies. Long bursts are 
associated to star forming regions with high star forming rate in spiral or irregular galaxies with 
low mass and low metallicity. On the other hand, short bursts seem to be associated to regions 
with low star formation rate, either inside a galaxy or a low star-forming elliptical galaxy char¬ 
acterized by high mass and high metallicity or even have evaluated end left the central part of 
their galaxy. These findings strongly suggests old stars or stellar remnants as the progenitors of 
short bursts. 


Ultra long burst: Within the diversity of GRBs, one of the extreme case is the so-called ultra 
long GRBs of typical duration in the order of thousands of seconds (e.g. 25000s for GRB 111209A 
my They are thought to originate from the collapse of a low-metallicity blue supergiant in 
rapid rotation US [ED] IHII- 


X-ray flashes (XRF): A class of sources, called X-ray flashes or XRFs, was characterized 
by the HETE-2 (High Energy Transient Sources Experiment) satellite in 2001 jSJj. They have 
similar soft spectrum and duration to that of long GRBs [E5j. Their total energy is small in the 
prompt phase and they are brighter in the X-ray emission. They appear to be correlated with 
supernova explosions. As a result, they might have the same origin as GRBs. Possibly, they 
are either a lower energetic subset of the same phenomena or the differences are due to their 
orientation relative to our line of sight. 
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Dark bursts: Another kind of GRBs is composed of the so called dark GRBs. They are defined 
by a large absorption (gas) and extinction (dust) in their host galaxies: they are characterized 
by a hydrogen column density larger than 10 23 cm -2 and by an extinction larger than 2.6 mag. 

Long bursts without SNe: Moreover, several long/soft GRBs ( e.g. GRB 060614) without 
an associated supernova have been discovered by Swift |84) . These events open the door for a 
so far unknown third class of GRBs, which challenges the idea of collapsing massive stars and 
binary mergers being the only progenitors of GRBs and suggesting that the tidal disruption of 
a star by a black hole would be an ideal way to power a long duration GRB [S3]. 

To conclude, the diversity of GRBs can reflect multiple progenitors and different types of 
interactions with their environments. 

1.3.7 Progenitors of GRBs 

Since the discovery of GRBs, many possible progenitors have been proposed. However, only a 
few remain, mainly because of the enormous energy budget required. Among the massive star 
progenitor models, two are popular: the first one is the collapsar model [II [IH EH [49], and the 
second one is the millisecond magnetar model ESI Em EH [S3 - 

Collapsar model 

A collapsar is a fast-rotating massive star about the collapse. It lost its H layer during its 
evolution by stellar winds. As it collapses onto a black hole, it creates an accretion disk. The 
nickel at the origin of the SN is created by the accretion disk, while the jet is created by accretion 
of the matter of the disk onto the black hole, and is collimated by the material of the star. 

Not all SNe and not even all most energetic ones (hypernovae) produce GRBs. This could be 
explained by the requirements on metallicity and rotational speed or by the fact that the jet is 
not pointing towards the Earth. However, it is yet unclear if all long GRBs are accompanied by 
an SN. For example, a coincident supernova for GRB 060614 has to be 100 times fainter than a 
standard SN, as imposed by the flux limits. This may indicate that long GRBs are composed of 
different populations. 

Millisecond Magnetar model 

A magnetar is a fast-rotating neutron star m • Its high magnetic field and high angular momen¬ 
tum provide the energy for the GRB and the SN. As the produced outflow is highly magnetized, 
a high radiative efficiency is expected [SSI . 

Most GRB light-curves show a plateau and some times flares. They are interpreted by energy 
injection in the outflow long time after the prompt phase. The millisecond magnetar model is 
able to explain this energy injection by accretion of matter onto the BH. 

Binary of Compact Object 

Two compact objects lose their energy and orbital angular momentum by gravitational wave 
radiation and merge as a result [28]. These binaries are composed by Neutron Star (NS) - NS, 
Black Hole (BH)-NS, BH-BH gHl El US HO] • 

Enough energy to power a GRB can be created in a binary merger. The duration of the GRB 
is comparable to the lifetime of the accretion disk. So in the binary mergers, the duration of the 
resulting GRB is expected to be some miliseconds while it is longer in the collapsar model. 
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1.4 Emission Theory of Gamma-Ray Bursts: The Fireball 
Model 

1.4.1 General description 

Nowadays, the standard model for the emission mechanisms is the so-called fireball model [su¬ 
it considered a large amount of energy (10 51 — 10 54 ergs) released in a very small region (r 0 ~ 
10 6 or 10 7 cm) as implied by the observed variability in the prompt phase. 

From the volume and the energy, it is possible to show that the region of emission is optically 
thick for pair creation: this is the so-called compactness problem m- It can be solved by 
considering that the emitting object is moving relativistically towards the observer. This is 
possible when the energy is much larger than the rest mass energy. The released energy can be 
dominated either by thermal energy or by magnetic energy, leading to two different acceleration 
mechanisms |931I94| . The Lorentz factor of the outflow is increasing up to the saturation radius. 
In the case of radiation-dominated outflows, it is defined as r s = tq x rj, where rj = Eq/(Mqc 2 ) 
is the ratio between the total energy Eq and the rest mass energy Mqc 2 of the ejecta. The 
coasting Lorentz factor is estimated to be as large as T = 100 - 1000. Once all energy has been 
converted to kinetic energy, it is necessary to find a way to convert it back to radiation with a 
high efficiency. This is achieved by shocks. 

This model is summarized on Figure fl.lSI ffor the typical radius and associated emission mech¬ 
anisms) and Figure IT. 161 (for the corresponding evolution of the Lorentz factor). The different 
possible processes are discussed in details below. 

Photospheric emission 

Under the conjugate action of the density decrease and the increase of the Lorentz factor, the 
opacity for Compton scattering is decreasing. As it drops below the unity, all photons initially 
trapped within the outflow are set free to reach a distant observer. This is called the photospheric 
emission. The outflow becomes transparent typically a radius of some 10 12 cm (see Figure IT. 1511 . 

The efficiency is between 5 and 30% of the total energy of the burst, and can even be larger 
in some specific cases (for example GRB 090902B jjJ5]). The spectrum of this emission is thermal 
(it has a nearly black-body shape) and is characterized by its normalization and temperature 
T 0 b s . However, a black-body spectrum is too narrow to account for most observations. Sub- 
photospheric (below the radius at which the plasma becomes transparent) dissipation were also 
studied as a way to broaden the expected thermal spectrum, and further increase the efficiency 
of the photosphere. 


Internal shock model 

As the emission of the outflow at the center is not steady, different parts of the plasma are 
moving with different velocities. When two parts of the plasma with different speeds collide, 
internal shocks form. They accelerate part (if not all) of the electrons in a power law and locally 
increase the strength of the magnetic field, at the expense of the kinetic energy. The electrons, 
which are accelerated by the 2nd order Fermi mechanism, radiate synchrotron emission in the 
induced magnetic field. Additionally, the synchrotron photons can be inverse Comptonized to 
even higher energy (synchrotron self-Compton). 

The typical radius for internal shock, rj s ~ rj 2 ro, is in the order of 10 14 cm (see Figure fT. 1511 . 
The main problem of the internal shock model is that the efficiency is low (5% to 20%), while 
the observations indicate an efficiency larger than 50% for most bursts [961 . 
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Figure 1.15: Standard fireball model. Image Credit: NASA /Swift 
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Figure 1.16: Schematic evolution of the jet Lorentz factor and examples of symbolic locations of 
radius: the saturation radius r s , photospheric radius r p h, internal shock radius Ti S and external 
shock r es [6]. 


External shock model 

The remaining kinetic energy (which was not dissipated by internal shock) is released when the 
outflow is decelerated by the interstellar medium (ISM), which can be either of constant density 
n ~ 1 cm~ 3 or the wind of the progenitor, n oc r~ 2 , where r is the radius. This deceleration 
creates the reverse shock (which propagates back into the ejecta [ST]) and the forward shock 
(which propagates into the external medium): these compose the so-called external shock model. 
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These two shocks are created as soon as the outflow starts to expend at the central radius 
?'o, but they become efficient only when the kinetic energy can be efficiently reduced, that is to 
say when the swept ISM mass is large enough. Such a condition is fulfilled at typical radii 10 16 
- 10 17 cm [97] (see Figure I1.15D . As for the internal shocks, the magnetic field is increased by 
plasma instabilities and accelerated particles radiate synchrotron in the induced magnetic field. 

GeV emission in External shock 

The GeV emission which is seen in the spectrum of some GRBs can be explained by different 
models. The first possibility is inverse Compton emission of synchrotron photon produced in 
the external shock and scattered by accelerated electrons (synchrotron self-Compton) [5%] . The 
second possible model is the pair loading of the surrounding medium of the burst by prompt 
photons. These photons are scattered by cold electrons of the ISM, then they interact with 
other photons of the prompt emission to create electron-positron pairs. The local ISM density 
is increased by a large factor [515] . These pairs are accelerated to relativistic speeds by the huge 
radiative pressure and inverse Comptonize the prompt MeV photons. Another explanation might 
be hadronic processes e.g. proton-proton interaction, neutron decay, proton-photon interaction 

ESS. 

Jet break 

In addition, GRBs are assumed to be collimated into jets. The relativistic beaming produces 
a visibility cone with an opening angle of 1/T, which increases as the outflow decelerates (see 
Figure [Ll6l) . When the opening angle is equal to the jet half-opening angle (1/r = 9j) the 
afterglow light-curve should break and decline more rapidly: this is called the jet break. 

The jet half-opening angle can be computed from the time of the jet break by assuming the 
standard model for the afterglow. There are two possibilities: 

• the afterglow interacts with the ISM of constant density, 

• the afterglow interacts with the wind of the progenitor whose density decreases propor¬ 
tionally to r -2 , where r is the radius. 

For the ISM of constant density, the half-opening angle is given by ITUTlfTU21 : 

9(tb, E i80 ) = 0.161 f ) n 1/8 (-^— ) , (1.3) 

\1 + zj \E iso , 52 J 

where the standard values for the number density of the medium n = 1 cm -3 is assumed. It 
also takes into account the radiative efficiency of the prompt phase, r/ 7 . Ei SO is the isotropic 
equivalent energy. In this study, each quantity expressed as X n is such that X = X n 10™. z is 
the redshift of source. Finally, the break time R is in days, 

1.4.2 Detailed afterglow theory 

Blast-wave evolution 

The following derivations can be found with more detailed in m- GRB afterglows can be suc- 
cesfully explained by the interaction between the outflow and the ISM, resulting in the decrease 
of the outflow speed [64]. The kinetic energy lost by the outflow is converted to kinetic energy of 
the shocked ISM and to internal energy. A fraction es of this internal energy is assumed to be 
in the form of magnetic field while a fraction e e is assumed to be in the form of kinetic energy of 
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electrons. Blandford and McKee solved the relativistic hydrodynamic equations of motion con¬ 
sidering adiabatic and radiative relativistic blast waves m ■ When the internal energy created 
in the collision is completely emitted, the evolution is said to be fully radiative, otherwise, it is 
adiabatic. The evolution of the Lorentz factor in a constant ISM for the adiabatic case is given 
by: 


T(t) 2* 5.4 


1 + * 
tday 


3/8 



(1.4) 


where t day is the time measured by the observer in days after the GRB, no is the ISM proton 
number density and K 52 is the total explosion energy. The evolution of the Lorentz factor in a 
constant ISM for in the radiative case is given by m- 




where T is the initial Lorentz factor of the outflow. 


(1.5) 


Electron distribution 


The shocks, created by the interaction of the relativistic outflow with the ISM, accelerate the 
electrons in a power law between the Lorentz factors 7 m i n and 7 max ■ The comoving distribution 
function of the shocked electrons is assumed to be: 


dNe{ 7 ) _ p 

-7- oc 7 

07 


7 min ^ 7 — 7 maxi 


( 1 . 6 ) 


where p is the electron injection index. 

Since the fraction of internal energy given to the electrons e e is prescribed, and since the elec¬ 
tron distribution function gives the kinetic energy of the electrons by integration, the minimum 
electron Lorentz factor can be computed and is given by: 


_ Wp 

r )min — 

m e 


v -2 

P~ 1 


(1.7) 


when assuming that all electrons are accelerated, 7 1 , p > 2 and 7 max 7 mm, and where m v 
and m e are the mass of a proton and of a electron respectively, and T is the Lorentz factor of 
the outflow. 

The electrons gain energy by the second-order Fermi acceleration process while they lose their 
energy by radiating synchrotron photons. To compute 7 max, the acceleration rate is set equal to 
the radiation-loss rate, which yields: 


Imax 


= 2 x 10 s 


e 


1/2 

max 


1/4 


n 0 T 1 / 2 ’ 


( 1 . 8 ) 


where e max is a constant in the order of the unity, (see Equation 11.44 in [1031 1. 

Finally, the cooling Lorentz factor y c corresponds to the Lorentz factor for which the rate of 
energy lost by synchrotron emission is equal to the rate of energy lost by adiabatic cooling: 


9 m e (1 + z) 

128 m p a T e B n 0 c F 3 t day 


(1.9) 


where or is the Thompson scattering cross section, c is speed of light and t is the observer frame 
dynamical time. 














20 


CHAPTER 1. INTRODUCTION 


Taking into account all different cooling rates, the electron distribution function is approxi¬ 
mated by: 


dN 8 _ x 


7i 


7i 


-(p+i) 


for 7 o< 7 < 7 i) 
for 7 l < 7 < 7 max- 


( 1 . 10 ) 


In the so-called slow-cooling regime , the magnetic field is weak so that the electrons do not 
cool below ^min by emitting synchrotron radiation: the cooling affects only the electrons in the 
high-energy tail of the distribution. In this case, 70 = 7 m j n , 71 = 7 c and s = p. Conversely, 
if the magnetic field is strong enough, the electrons are cooled by synchrotron radiation to 
Lorentz factor smaller than 7 m in on the dynamical time scale of the system. This regime, 
called fast-cooling regime, is characterized by 70 = 7 C , 71 = 7 min and s = 2. Both regimes are 
summarized on Figure [1.171 




Figure 1.17: Electron distribution function for the slow (left) and fast-cooling (right) assuming 
p = 2.5. 


Synchrotron emission 

The observed characteristic synchrotron frequency is given by [ 1051 : 

e B y 2 r 

v = ----, 

2 7 r m e c (1 + z) 


( 1 . 11 ) 


where e is the elementary charge, 7 is the Lorentz factor of the electron, B is the comoving 
magnetic field and T is the Lorentz factor of the outflow. 

Using Equation 11.111 the cooling frequency, v c , is obtained from the corresponding cooling 
Lorentz factor 7 C for an adiabatic blast wave: 

i/ c = (9x 10 12 Hz)( 1 + z)- 1 / 2 eB 3 / 2 n“ 1 E 5 " 2 1 / 2 t” 1/2 , (1.12) 


The typical synchrotron frequency (or injection frequency), v m , is calculated by following the 
same way and considering 7 m in- 


v m = (6 x 10 15 Hz)(l + 2 


*>‘ /2 (^f y 


2 1/2 ^1/2 -3/2 

e e t B rj b2 b d > 


(1.13) 
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Another characteristic frequency is introduced: v a is the transition frequency below which the 
photons are absorbed by synchrotron self-absorption. It is called the synchrotron self-absorption 
frequency. Its evolution is given by v a oc t -1 / 2 for an adiabatic blast wave and v a oc t~ 4 / 5 for a 
radiative blast wave in the fast-cooling regime. However, the slow-cooling regime is characterized 
v a oc const in the adiabatic case and v a oc t ~ 3 / 35 in the radiative case. As an example, for 
adiabatic expansion in the slow-cooling regime, v a is given by: 

v a = (2 x 10 9 Hz)(1 + ^) _1 e7 1 4 /5 « 3/5 - E 52 5 , (1.14) 


The synchrotron spectrum is composed of four segments separated by the typical frequencies, 
v c , v m and v a . In the slow-cooling regime corresponding to v m < v c , the spectrum is given by 
0 : 


' (v a /v m ) 1/3 (v/v a ) 2 
(v/v m ) 1/3 

w^r (p ~ 1)/2 

. {v c /V m )~ ( ' P ~ 1),2 {y/V c )~ p/2 


V a < V < V m , 
V m <V< V m , 
Vc kk V V-uiax, 


(1.15) 


where F Utmax = (20mJy)(l + z)e l J j ~E §2 2 2 $ and v max is the maximum synchrotron frequency 
computed from 7 max . It is also usually assumed that v a < min(i/ m , v c ). 

In the fast-cooling regime corresponding to v m > v c , the spectrum is given by 


' (v a /v c ) 1/3 {v/v a ) 2 
(■ v/v c ) 1 /3 
(v/v c )~ 1/2 

. ( v m /v c y 1/2 (v/v m )~ p/2 


V a < V < V c , 
v c <v< v m , 

Vm A V y V m ax • 


(1.16) 


The spectra for the two regimes are illustrated in Figure H. 181 



Figure 1.18: Synchrotron spectrum in the two regimes: slow-cooling (left) and fast-cooling (right) 

13- 


In the dynamical evolution, the outflow goes from initially fast to slow-cooling regime. The 
light-curves at a given frequency v * can be computed by considering the time to which is the 
transition time between the fast and slow-cooling regimes. It is computed by setting v c = v m 
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and is given by (Tj: 

{ 210 e 2 B C e E 5 2 n\ days, adiabatic, 

4.6 e 7 J 5 e 7 J 5 days, radiative. 

The corresponding frequency v o is such that uq = v c (to) = v m {ta) and is given by p3]: 

f 1.8 x 10 11 e~ R 5 ^ 2 e~ 1 E^n^ 3 ^ 2 Hz, adiabatic, 

radiative. 


(1.17) 


i/ 0 = 


8.5 x 10 


10 


Hz. 


(1.18) 


Neglecting synchrotron self-absorption, there are two possibilities: v* > vq or z/* < vo- The 
first possibility (referred to as high-frequency light-curve) implies to > t m > t c , where t m and t c 
correspond to the times when z/* = v m and z/* = v c respectively. The light-curve at frequency z/* 
is composed of four segments, respectively: 

• z/* < v c (labeled B, fast-cooling), 

• u c < z/* < u m (labeled C, fast-cooling), 

• v*> v m > v c (labeled D, fast-cooling), 

• z/* > v c > Pm (labeled H, slow-cooling: because v c > z/ m ). 

The second possibility, zq, < zz 0 , referred to as low-frequency light-curve implies to < t m < t c . In 
a similar way, the light-curve at frequency v* is composed of four segments, respectively: 

• z/* < v c < v m (labeled B, fast-cooling), 

• z/* < v m < v c (labeled F, slow-cooling), 

• v m < < v c (labeled G, slow-cooling), 

• z'm < z'c < v* (labeled H, slow-cooling). 

The light-curve is represented in Figure H. 191 taken from [9]. 



Figure 1.19: Synchrotron light-curves in the two regimes: low frequency (left) and high frequency 
(right) E] . 


In all calculations above, it is assumed that the surrounding medium is a uniform ISM. In 
the case of a wind, which corresponds to n oc r - , all equations can be established in a similar 
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way. The shape of the spectrum is unchanged, only the time evolution of the cooling, injection 
and absorption frequencies are changed ITMI : 


v c = 1.1 x 10 14 EH 2 A: 2 e~ 3l l T] /2 Hz, 


'-£,-2 


(1.19) 


where E $2 is the isotropic energy in units of 10° 2 ergs, is the initial number density in the 
wind, Td is the time expressed in days after the burst. 


= 0.6 x W 13 El y { 2 el_, e~ 3/ _\ T d 1/2 Hz, 

3.5 x 10 12 E~ 2 2/5 Al 1/5 e 6 J 5 _ 2 T~% 5 Hz 
9.3 x 10 9 E^ 2/5 A 6 J 5 e~ 1 _ 1 e^_ 2 T~ 3/5 Hz 

Being able to distinguish between two environments is imp( 
the progenitor of GRBs. 

Closure relations: The broadband spectrum and corresponding light-curve of synchrotron 
radiation can be derived from a power-law distribution of electrons accelerated by relativistic 
shocks as explained above. The power-law parameters of the spectrum and the light-curve can 
be obtained by modeling the flux F v as: 



radiative, 

adiabatic. 


( 1 . 20 ) 

( 1 . 21 ) 


ortant to determine the nature of 


F v oc t~ a p-P, (1.22) 

where t is the observer time, v is the observer frequency, a and /? are respectively the temporal 
and the spectral indexes, see Equation [TTT5] and Equation II. 161 

These spectral and temporal decay indexes can be combined into several closure relations for 
synchrotron emission ( 1071 9 1108 : 1091 2Tj. These relations can be used to compare theoretical 
models with the afterglow observations in order to investigate the geometry of the burst, its 
surrounding medium, the microphysics of the fireball, and its cooling state. In Table 11.21 the 
closure relations are recalled for an homogeneous interstellar medium, for a wind and for a jet 
in the slow and fast-cooling regimes, taking into account p >2 (standard decay index for the 
accelerated electrons). 
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Table 1.2: Closure relations, temporal index a and spectral index p in various afterglow models 
where the convention F„ oc t a v& is adopted ESI- And validation of closure relations in X-ray 
and optical band. 



p >2 

X-ray 

Optical 

ISM, slow-cooling 
v a < v < v m 

a(/3) 

a = 1/2, /3 = 1/3 

no 

yes 

v m < v <v c 

a = 3/3/2 

yes 

yes 

v > v c 

a = (3/3 + l)/2 

yes 

no 

ISM, fast-cooling 

v a < v <v m 

a = 1/6, p = 1/3 

no 

yes 

v > v m 

a = (3/3 + l)/2 

yes 

no 

Wind, slow-cooling 

v a < v < v m 

CO 

i—1 

II 

'Q. 

o' 

II 

S 

no 

yes 

v m <v <v c 

a = ( 3/3 - l)/2 

yes 

yes 

v > v c 

a = (3/3 + l)/2 

yes 

no 

Wind, fast-cooling 

Va < v < V m 

a = -2/3, P = 1/3 

no 

yes 

v > Vm 

a = (3/3 + l)/2 

yes 

no 

Jet, slow-cooling 

V a < V < V m 

a = -1/3, P = 1/3 

no 

yes 

V m < V < V c 

P 

II 

to 

to 

1— L 

yes 

yes 

V > v c 

a = 2P 

yes 

no 


1.5 Purpose of the Thesis 

The work presented in this thesis attempts to strengthen our understanding of the diversity of 
GRBs which can result from their progenitors or their environments. 

It is known and widely accepted that GRBs are associated to SNe. Unlike SNe, which are 
precisely classified into different types, mainly depending on the composition of the star, GRBs 
are only grossly classified into two sub-types depending on the duration. However, a more detailed 
classification would allow to better understand the diversity of GRBs and better constrain their 
environment(s) and progenitor(s). 

It is thought that the long GRBs originate from the core collapse of a massive star. Core 
collapse SNe are very diverse, depending on the mass of the progenitor, its size and the remaining 
elements in its outer layers. Thus, different kinds of GRBs might be determined based on the 
same criteria: the mass of the progenitor, its size and the remaining elements in its outer layers. 
Indeed, many differences can be found, which have led to the creation of groups: Ultra-long 
GRBs (based on the duration), X-ray flashes (based on the spectral hardness of the prompt 
emission) and dark bursts (based on the absorption and extinction of the afterglow light). These 
differences (and these groups) are also seen in the subclass of GRBs firmly associated to SNe. 

Even with deep observations of bright GRBs, their origin is still unknown. In this thesis, 
I present the results of my work, dedicated to a better understanding of the GRB diversity. 
In particular, I studied the faint GRBs. I defined a new sample of GRBs, based on the X-ray 
afterglow luminosity. The bursts in that sample are called low-luminosity afterglow (LLA) GRBs. 
The highlight was set on contrasting the properties of LLA GRBs to those of normal long GRBs. 

Indeed, using data from BeppoSAX, it was discovered that the X-ray light-curves of long 
GRBs define several well-separated groups once the distance effects are corrected [ITO] . This 











1.5. PURPOSE OF THE THESIS 


25 


was later on confirmed by extending to the samples from the XMM-Newton and Chandra data 
m- Three groups were determined based on the empirical properties of the X-ray afterglows, 
namely group I (mean flux: 7.0 x 10 _12 erg.s _1 cm~ 2 at one day), group II (mean flux: 3.1 x 
10“ 13 erg.s _1 cm _2 at one day) and group III (all other bursts less luminous than group II events) 
m- The focus was put only on groups I and II, in order to explain the observed clustering. 
Indeed, group III events were too few in order to perform a meaningful statistical study. However, 
thanks to the Swift satellite jll2j . the sample of available bursts has grown and it was possible 
to resume the study of group III namely low-luminosity afterglow (LLA) GRBs [72]. 


Hereafter, I give the outline of my thesis. 

In Chapter 2, the data reduction processes of Type lib SN 2004ex on the photometric and 
spectroscopic data (in the optical band) are presented. The results are compared with the 
prototypical spectra and light-curves of type II SNe (SN 1993J, SN 2008ax) and with the light- 
curves of type lb (SN 2007Y) and type Ic (SN 19941) SNe. The spectral lines of SN 2004ex are 
identified by using the spectral similarities with SN 1993J. The velocity, mass and kinetic energy 
of the ejected H are calculated. Finally, the physical properties of SN 2004ex are computed by 
using the similarities with the light-curve of SN 2008ax. 

In Chapter 3, the data processing of GRBs in the X-ray band is explained. The selection 
method of the LLA GRBs sub-sample is given. Their spectral and decay indexes are computed. 

In Chapter 4, the statistical study of LLA GRBs is presented. First, the redsliift distribution 
is analyzed and it is found that LLA GRBs are in average closer than normal lGRBs. Second, 
different selection effects are discussed and it is proved that LLA GRBs do not suffer from them. 
Third, their afterglow properties are compared by using the closure relations in both X-ray and 
optical bands. Finally, the prompt properties are discussed in light of the Arnati correlation. 

In Chapter 5, I discuss many of the other properties of LLA GRBs: rate density, host galaxy 
and possible progenitor. The conclusions of the thesis and the perspectives follow in Chapter 6. 
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Chapter 2 


Observation and Data Reduction 
Applied to SN 2004ex 


In this chapter, I will present the photometric and spectroscopic analysis of the type lib supernova 
2004ex during the transition (30 days < t < 60 days) and nebular (t > 60 days past explosion) 
phases. The purpose of this work was to complete the analysis of the data from SN 2004ex taken 
by the 1.82 m Copernico Telescope of Mt. Ekar, 3.5 m Telescopio Nazionale Galileo (TNG) and 
2.2 m MPG/ESO telescope. 

The results are compared with the prototypes of type lib: SN 1993J and SN 2008ax respec¬ 
tively. I find that the light-curve of SN 2004ex is very similar to the two prototypes light-curves, 
however it is closer to the one of SN 2008ax. On the other hand, its helium lines (He I A = 4394, 
6561, 6914, 7168 A) are still well-detected 4 weeks after the explosion. However, the hydrogen 
(H a ) at 6286 A is weaker than the one found for SN 1993J at a similar period. 


2.1 The intermediate characteristics of type II SNe 

Different kinds of SNe can be identified by performing spectroscopic analysis and considering 
the elements synthesized during the SN evolution and explosion. An other criterion is about the 
brightness of the SN which can be obtained by performing a photometric analysis. 

SNe were classified into type I and type II according to the absence or the presence of H 
lines respectively [39). Type II SNe only occur in spiral galaxies. They are thought to be the 
result of the core collapse of massive young stars [113] . The progenitors of core-collapse SNe 
are massive stars, either single or in binary systems, which have completed the nuclear burning 
stage. Different types of core-collapse SNe (lib, IIP, IIL, lb, Ic) have been identified based on 
their spectroscopic and photometric properties, forming a sequence explained by the progenitor 
mass loss history. 

Recently, several SNe with intermediate characteristics have been discovered, suggesting a 
smooth transition and requiring the introduction of the hybrid classes type lib, Ic and Ib/c SNe, 
the latter being associated with GRBs [1141 . The spectra of lib SNe are dominated by H lines 
at early times while He lines become prominent at late times. On the one hand, spectra of type 
Ic SNe show neither He I lines nor H a , while the spectra of lb SNe show Hel lines and no H a . 

Furthermore, the similarity between the light-curves of type lib with those of type lb events 
[115] , in addition to the known spectral similarities at late times and the similar peak radio 
luminosities also suggests that these two types of events might come from similar progenitor 
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systems. Their association with GRBs also supports the idea that SNe Ic and Ib/c are related. 
Thus, the intermediate characteristics of type-IIb SNe supports the idea of a continuous sequence 
of SNe having different envelope mass (i.e. types II, lib, lb, Ib/c, Ic [llfij l. 

However, it was investigated with the type lib and type lb that they could be associated with 
interacting binaries inn- This idea was also studied by several other authors: Blinnikov et al. 
m showed the relevance of this possibility by the modelling of SN 1993J and Bersten et al. 
m tested this idea on SN 201 ldh by using single and binary progenitors modellings of type 
lib SN. This issue is also widely discussed by Dessart et al. |T5ailT5Tl who simulated binary-star 
models for the production of SNe Ilb/Ib/Ic and proposed that the progenitors of SNe lib and lb 
should have main-sequence masses smaller than 25 M 0 and be in a binary system whose stars are 
close to each other. However, SNe Ic should be the result of a more massive single star because 
of the lack of He I lines in their spectra. 


2.2 Introduction of SN 1993J and SN 2008ax: comparative 
sample 

2.2.1 Properties of SN 2008ax 

SN 2008ax has been discovered just 18 days before its maximum peak in the NGC 4490 at 
redshift 0.001855. It shows a rapidly declining light-curve up to 40 days. The amount of 56 Ni 
synthesized in the explosion is between 0.07 and 0.15 M 0 . The kinetic energy of SN 2008ax 
estimated to be 1 — 6 x 10 51 ergs, while its total ejecta mass is 2 — 5 M 0 . The ejecta velocity is in 
the range of 23 000 - 26 000 km.s -1 . The ejected element H a has very high early-time velocity 
of 13 500 krn.s^ 1 , which is rapidly decreasing until the day 14 in velocity evolution (for more 
information, see 1T1.1122 j '). 

From the spectroscopic and photometric evolution of SN 2008ax, the progenitor constrained 
to be a W-R star of type WNL, which is either a single massive W-R star or a lower-mass W-R 
star (main-sequence mass 10 -14 M 0 ) in an interacting binary system HHH22I- 

2.2.2 Properties of SN 1993J 

The first ever SN lib event identified was SN 1987K m; another more recent example is 
SN 1993J which is a nearby event observed in the M81/NGC 3031 galaxy at redshift 0.000113 
[124] , The light-curve of SN 1993J was unusual with a narrow peak (shock break out) followed by 
a secondary maximum (optically thick region), which is similar to the one of SN 1987A. After a 
rapid luminosity decline around 50 days after the explosion, the light-curve of SN 1993J showed 
almost well-fitted exponential tail (optically thin region) with a decline rate faster than normal 
SNe II and similar to that of SNe la. It is remarkable that the late time spectra, except for the 
presence of He I lines, appear similar in all core-collapse SNe types. 

The progenitor of SN 1993J was a supergiant star of spectral type K m- It was the massive 
member of a binary system and its mass was in the range 12 to 17 M 0 . In this binary system 
the two stars had comparable main sequence masses j l261 . 

2.3 Supernovae: SN 2004ex 

SN 2004ex discovered in the galaxy NGC 182 on October 10 th , 11.34 UT, 2004. The SN was 
firstly discovered by Tenagra II, a 0.81 m telescope, with an unfiltered magnitude of about -17.7. 
It was later observed by the LOSS telescope on October 13.33, UT 2004. Its celestial coordinates 


















2.4. PHOTOMETRY 


29 


are a= 00h38ml0.19s, 5= +02h43ml7.2s (equinox 2000.0). It is located 33” West and 25.3” 
South of the nucleus of its host galaxy [127], which is a SBa barred spiral galaxy, at redshift 
2 = 0.01755 obtained by the host galaxy spectral observation, as it can be seen in Figure 12.11 
(left). 


Susp. 00h 38m 10.187s 02“ 43' 17.18" (arrowj^^ 

Hi 

NGC182 

\ ; 

!• 

From star 1 : • • 

9.159s West 
92.84" North 

(J2000 and Astroart II Astrometry) 



Figure 2.1: Optical image of SN 2004ex with its host galaxy, NGC 182, (left) and the map of 
possible standard stars (right) used for the magnitude calibration. 

The optical photometric datasets (in BVRI bands) of SN 2004ex were acquired in 6 post¬ 
explosion epochs, while spectroscopic observations were performed for 4 post-explosion epochs. 
They were obtained using the following instruments: 

(i) Device Optimized for the LOw RESolution (DOLORES), with a scale of 0.252 arcsec.px -1 
and a field of view of 8.6 x 8.6arcmin 2 , at the 3.58 m Telescopio Nazionale Galileo (TNG) (La 
Palma, Spain), 

(ii) The Wide Field Imager (WFI), with a scale of 0.238 arcsec.px -1 and a field of view of 
34 x 33arcmin 2 , at the 2.2 m European Southern Observatory Telescope (La Silla, Chile), 

(iii) The Asiago Faint Object Spectrograph and Camera (AFOSC), with a scale of 0.46 arcsec.px - 
a field of view of 7.8 x 7.8arcmin 2 and range 355-780 nm, resolution 2.4 nm, mounted on the 
1.82 nr Copernico Telescope of Mt. Ekar (Asiago, Italy). 

The observations and the data reductions are presented in the subsection 12.41 for photometry 
and m for spectroscopy. The results are compared with two other type lib SNe: SN 1993J and 
SN 2008ax. The conclusion follows. 


2.4 Photometry 

2.4.1 Data analysis 

The photometric datasets were collected in the wavelength range from 3650 to 8060 A which 
corresponds to the Johnson-Bessel filters UBVRI at all 6 epochs. Table 1 shows the results of 
the observations as well as the telescopes used. 

The datasets were pre-reduced according to the classical prescriptions ( i.e. overscan, trim, 
bias, flat-field correction etc.) in order to remove the instrumental effects. The photometric 
analysis was performed by using the Queen’s University supernova Belfast Archive (QUBA) 
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semi-automatic pipeline Point Spread Function (PSF) module based on the Image Reduction 
and Analysis Facility (IRAF) 11281 . The dataset acquired for the photometric night of December 
17 th , 2004 obtained by the WFI instrument, was calibrated by using the Stetson extension of 
the Rul49 Landolt catalog [1291 [1301 - A list of five local standard stars was built. They are 
labeled by 1, 2, 4, 10, 12 in Figure [ATI (right). This list was used for the calibration of the 
other non-photometric nights and to calculate the magnitude of the supernova by considering 
the contribution of the galaxy to the background. It should also be taken into account that the 
contribution of the galaxy is spacially rapidly changing when the SN is on one of its arm (see 
Figure HOT) . Indeed, the contribution of the galaxy is larger if the SN is on one of its arm than 
when it is not. 


2.4.2 The light-curve of SN 2004ex 

The output of the PSF-fitting with the QUBA pipeline on the Johnson-Bessell (UBVRI) pho¬ 
tometry of SN 2004ex is shown on Figure 12721 and is presented in Table [All together with the 
error on the magnitude, estimated by the PSF-fitting technique. 



JD+2453300 


Figure 2.2: The 4-band optical light-curve of SN 2004ex. The x-axis is in JD days and the y-axis 
is the apparent magnitude. 


Date 

JD + 2453300 

B 

V 

R 

I 

Inst. 

Now. 16 th 

26.39 


17.71i0.002 

17.86i0.04 

17.39i0.05 

Ekar 

Dec. 9 th 

49.52 

20.45±0.04 

19.17T0.03 

18.55i0.02 

17.96i0.02 

WFI 

Dec. 12 th 

54.55 

20.52±0.03 

19.24i0.03 

18.7i0.03 

18.08i0.02 

WFI 

Dec. 17 th 

57.54 

20.51±0.03 

19.28i0.03 

18.7i0.03 

18.09i0.02 

WFI 

Jan. 1 th 

75.28 


19.07i0.08 

19.55i0.07 

18.33i0.07 

Ekar 

Jan. 17 th 

88.21 

21.38±0.15 

19.91i0.08 

19.41i0.07 

18.57i0.06 

Ekar 


Table 2.1: The photometric observations. The Johnson-Bessell UBVRI photometric results were 
obtained with the PSF-fitting technique. 
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Since data used for this analysis was taken during the transition (30 days < t < 60 days) 
and the nebular (t > 60 days past explosion) phases, it is not possible to obtain important 
properties: luminosity peak, width of the peak, the nickel mass, energy of the ejecta. However, 
the light-curve of SN 2004ex is complete in the R-band when adding the results from Arcavi et al. 
m obtained during the photospheric phase which is in the optically thick region (t < 30 days 
past explosion 11311 1 as shown in Figure 1^31 As a result, the peak magnitude is —17.2 and the 
brightness of the source decreases rapidly after 21.5 days. The explosion date is 53289 MJD 
(October 11 th , 2004). 
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R-band in this work ♦ 

R-band in I. Arcavi 2012 ♦ 
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Figure 2.3: The optical light-curve of SN 2004ex in the R-band. The photometric results from 
Arcavi et al. m are in red, while the results of this analysis are in blue. 


The light-curves are usually fitted by a polynomial function of order 4. Here it is performed 
using the Legendre polynomials base defined as: 


f n (x) = ai F’i(x) + a 2 F 2 {x) H-b a(n) F n (x ) (2.1) 

where the a n are coefficients to be fitted for, the F n (x) are the Legendre polynomials containing 
terms of order x( n_1 ). 


Fi(x) = 1, 

F 2 (x) = x, (2.2) 

F n {x) = [(2 n - 1) F n _ i(x) x - (n - 1) F n _ 2 (x)] /n. 

The fourth and third-order Legendre polynomials fit parameters are presented in Table f2~?l 
The best fit to the light-curves are shown in Figure HOI The third-order polynomial was used 
because the fit obtained by the fourth order was not good and especially the peak could not 
be properly reproduced. The maximum magnitude of SN 2004ex (—17.27 mag) obtained by the 
third-order Legendre polynomial is compatible with the result (—17.2 mag) obtained by 1171 . 
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Table 2.2: The fourth and third-order Legendre polynomial fit parameters and maximum mag¬ 
nitude; 


Legendre 

ai 

a 2 

a 3 

a\ 

Magnitude mox 

polynomial 





mag 

fourth order 
third order 

-1.85T0.5 

-I.23i0.12 

0.01T0.002 

0.05A0.0005 

(-4.3il.2)xl0 -5 

(-1.1±0.12) x 10 -5 

(9.4i2.7)xl0 -y 

-17.6 

-17.27 


2.4.3 The comparison with other type lib SNe 

In order to better constrain the source, the light-curve is compared with the ones of two other 
sources: Figure [53] presents the comparison of SN 2004ex with the prototypical type lib SN 1993J 
and the well-known type lib SN 2008ax. These light-curves were corrected by the distance 
modulus which is given by 


fi = M — m = 5logio(d) — 5 + Ay (2-3) 

where M and m are the absolute and the apparent magnitudes respectively and d is the distance. 
Ay is the absorption and it is given by: Ay = 3.2 x E(B — V ) where E(B — V ) interstellar 
reddening and it is 0.022 for SN 2004ex. A distance estimation gives d = 70.55 Mpc for NGC 182 
(SN 2004ex) from the redshift calculation with the recession velocity of NGC 4490 (SN 2008ax) 
corrected for Local Group infall into the Virgo Cluster of 5240 km.s' - 1 m by adopting a 
cosmological model with H 0 = 70 km.s - .Mpc -1 . The distance modulus for NGC 4490 is // = 
29.92±0.29 mag (d = 9.6 Mpc) [TT] and jj, = 28.06 mag (d = 4.11 Mpc) for M81/NGC 3031 
(SN 1993J). The time was also rescaled by defining the peak time as t p = 0 for each SN. The 
comparison is made in the R-band. 

As it is shown that SN 2004ex is similar to SN 1993J regarding its spectral properties and 
SN 2008ax and when it comes to its photometric properties, the initial mass of its progenitor 
can be between 10 - 17 M 0 . SN 2004ex is a stripped-envelope core-collapse SN. Because of the 
many similarities with SN 2008ax, it is natural to consider that the same amount of 56 Ni was 
synthesized in the explosion (between 0.07 and 0.15 M 0 ). An upper limit for the kinetic energy 
of SN 2004ex can be then be estimated as being the kinetic energy of SN 2008ax: 1 — 6 x 10 51 ergs. 




JD+2453000 JD+2453000 


Figure 2.4: The fitted optical light-curve of SN 2004ex. Optical R-band light-curve of SN 2004ex 
(red dots) together with the best fit from a Legendre function (green continuous line). To the 
left (right), the fourth-order (third-order) Legendre polynomial is considered. 
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Figure 2.5: The optical light-curves of SN 2004ex compared with the light-curves of type lib 
SN 1993J [TO] and SN 2008ax [TT] in the R-band. 


Table 2.3: The physical parameters of three SNe (SN 1993J, SN 2008ax and SN 2004ex). 


SN 

Type 

M.B ,max 

mag 

E(B-V)tot 

be Ni mass 

Mq 

Ejecta mass 

Mq 

E kinetic 

10 51 erg 

Ref. 

SN 2008ax 

lib 

-17.32i0.50 

0.4i0.1 

0.07- 0.15 

2 - 5 

1 - 6 

.133 

SN 2004ex 

lib 

-17.71i0.002 

0.022 

~ 0.07 — 0.15 

- 2 - 5 

~ 1-6 

this work 

SN 1993J 

lib 

-17.23i0.50 

0.2 

0.10 

1.3 

0.7 

CO 


2.4.4 The comparison with other types of SNe 

In Figure 12.61 the light-curve of SN 2004ex is compared to the light-curves of type Ic SN 19941 
[fT2l and Ib/c SN 2007Y [13] 0. The distances 6.2 1133 and 19.31 Mpc [13] were adopted for each 
of them respectively. The light-curve of SN 2004ex is similar to that SN 2007Y, as expected 
since type lib SNe and lb SNe have similar light-curves. 


2.5 Spectroscopy 

2.5.1 Data analysis 

Spectroscopic datasets have been acquired at 4 epochs, respectively 36, 52, 60 and 86 days after 
the explosion, by using the 1.82 m Ekar and the 3.5 m TNG telescopes (see Table im . The 
datasets consist of 2 spectra collected in the transition phase and 2 during the nebular phase. 

Spectra were pre-reduced and calibrated according to standard methods using IRAF routines 
( onesdspec, ccdred and specred packages ): the raw images were de-biased, overscan-corrected, 
trimmed, and flat-field-corrected after normalization of the flat-field image along the dispersion 
axis before the extraction of the spectra. During de-bias, firstly the master bias image was created 


1 The type of SN 2007Y is still debated: type lb [131 or type Ib/c 11151 
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-19 

-18 

-17 


-14 

-13 


w 


•• \ 


1993J lib 
2008ax lib 
2004ex lib 
2007Y lb 
19941 Ic 




20 40 60 

Days since maximum 


2004ex in this work 
2004ex in I. Arcavi 2012 
2007Y lb 
19941 Ic 




20 40 60 

Days since maximum 


Figure 2.6: Left: The optical light-curves of SN 2004ex (blue diamonds)(blue diamonds) com¬ 
pared with four light-curves of type lib (SN 1993J, (black squares) and SN 2008ax (magenta 
circles), type Ic SN 19941 (red crosses) [TJ] and type Ib/c SN 2007Y (cyan squares) [13j in the 
R-band. Right: The optical light-curves of SN 2004ex (red and blue diamonds) compared with 
two light-curves of type Ic SN 19941 (black square) ii 1(2] and Ib/c SN 2007Y (magenta crosses) 
m in the R-band. 


Date 

JD 

Day 

Inst. 

Grism 

Phase 

Nov. 16 th 

2453326.39 

36 

Ekar 

Gr4 

transition 

Dec. 1 th 

2453341.39 

52 

TNG 

Gr4 

transition 

Dec. 9 th 

2453349.52 

60 

Ekar 

Gr4 

nebular 

Jan. 4 th 

2453375.32 

86 

Ekar 

Gr4 

nebular 


Table 2.4: Spectroscopic observations of SN 2004ex. 


by combining all available bias images in order to improve the statistic and reduce the random 
noise; then it was subtracted to the source images. The overscan correction was made firstly by 
determining the mean bias level of each image in the overscan region and then by removing it 
from the all images. After the overscan correction, the trim was applied to cut off the affected 
region from the edge. At the end, the flat-field correction was made by using a master flat which 
is a combination of all flat images, in order to increase the signal-to-noise. Lastly, sky-flats were 
used to remove artifacts such as cosmic rays and stars. 

The one-dimensional spectra were obtained by mean of optimized extraction across the dis¬ 
persion and by subtracting the galaxy contribution with the apall package. After the definition 
of both the aperture and the background region at one wavelength, the task traces the position 
of the aperture on the bidimensional image. With this procedure, the night sky lines are also 
removed. The extraction of lamp spectra (usually Ne-HgCd, He-Ar lamps), obtained in the same 
instrumental configuration, allows the wavelength calibration with the apsum package. In most 
cases, the error on the wavelength calibration is less than 2 A. 

The next step is the flux calibration. The response curve for the given instrumental configu¬ 
ration has been obtained by spectroscopically observing standard stars with the same telescopes. 
The flux calibration is typically accurate within 20%. 

In my analysis, the spectra were corrected for the redshift (z = 0.01755) of the host galaxy, the 
reddening was eliminated using the deredden task in IRAF program and the external extinction 
value (host galaxy reddening, E(B-V) = 0.022) has been calculated according to the celestial 
coordinates and observation time [T3'6 j. The spectral evolution is shown on Figure [2771 at each 
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epoch (phases 36, 52, 60, 86). In these figures, the evolution of the elements can be followed 
during each phase. The dominant elements, especially He I, are present until the phase 86. 



Figure 2.7: The spectroscopic datasets of SN 2004ex at four phases. 

The identifications of spectral lines of SN 2004ex were originally performed the absorption 
lines and some emission lines by using the sarith , identify and splot tasks which are in the specred 
package in IRAF. 

2.5.2 The spectral result of SN 2004ex 

Then the position of the lines are compared with the typical identifications given T24I for 
SN 1993J and for SN 2008ax (TT [ 1133] . These lines are presented in Table 12751 and shown in 
Figures 12.81 12.91 12.101 12.111 for each phase respectively, together with a comparison with the 
spectrum of SN 1993J. Finally, Figure f?T51 shows a comparison of the spectra of SN 2004ex and 
SN 2008ax for all phases. 

2.5.3 The velocity and mass of hydrogen 

The spectra of SN 2004ex show both prominent He I lines and a relatively faint H a line, all 
with P-Cyg profiles nnzi- The wavelengths are shown in Table [53] Phase 52 of the SN 2004ex 
spectrum closely resembles that of SN 1993J at 41 days after the explosion m- This therefore 
defined the spectroscopic classifications of SN 2004ex as type lib supernova and confirmed the 
result of [137]. 

The velocity of the element can be calculated by using P-Cyg profiles or by considering balmer 
a transition (used in my computation): 

^ (c\ A \ 

V — %dopp x C — “ X c (2.4) 

where A = 6563 A is the wavelength of the balmer a transition, X a is the measured absorption 
wavelength, c is the speed of light and Zdopp is the doppler redsliift. The wavelength of the Hq, 
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absorption line at phase 52 is 6286 A. The velocity of the hydrogen envelope is estimated to be 
as a result 12662 km/s in this phase. However, in the early observation times, the H Q absorption 
minimum velocity found to be around 14000 km/s 1138 . which is in the range 10000-15000 km/s 
expected for type lib SNe. 

The mass of the ejected element can be calculated by Ek in = | m v 2 . When assuming that 
the kinetic energy of SN 2008ax is a lower limit for the total energy of SN 2004ex, the mass of 
the hydrogen envelope is estimated to be around 0.6 M 0 (where Eu n is 10 51 erg and vu a is 
12662 km/s). 


Epochs 

A 

36. 

2004ex 

37. 

1993J 

41. 

1993J 

52. 

2004ex 

52. 

1993J 

60. 

2004ex 

62. 

1993J 

86. 

2004ex 

89. 

1993J 

O I 

7254 

7144 



7439 


7288 


7215 


Ca II 

7202 




7238 


7114 




He 1 

7283 

6717 



7116 

7132 

6967 

7161 

6788 


He I 

7065 

6573 

6882 

6844 

6878 

6892 

6765 

6904 

6660 


He I 

6678 

6204 

6523 

6540 

6510 

6551 

6412 

6549 

6335 


H a 

6563 

6004 

6343 

6346 

6286 

6353 

6187 

6358 

6154 


He I+Na I 

5876&5892 

5519 

5736 

5732 

5688 

5737 

5591 

5743 

5547 


O I 





5482 


5372 




Fe II 

5269 



5077 

5084 

5083 

5013 

5086 


5090 

He I+Fe II 

5015&5018 

4985 

4921 


4925 

4923 

4834 

4915 

4769 

4925 

He I+Fe II 

4921&4925 


4828 

4922 


4831 


4837 



Ha 

6861 


4720 

4833 

4715 

4724 

4631 

4337 


4724 

Mg 

4571 




4612 






He 1 

4471 

4619 

4429 


4394 

4395 

4309 

4397 

4337 

4410 

He 1 

4437 


4116 

4723 







H 7 

4340 


4828 


4225 

4237 


4218 


4222 

Ca 11 



4720 

4219 

3821 







Table 2.5: The spectral lines of each element obtained for all different phases of SN 2004ex and 
SN 1993J. The intrinsic wavelengths of the line were also given for some elements. 


2.5.4 The spectral comparison of type lib SNe 

The overall characteristics of SN 2004ex reminds those of SN 1993J and SN 2008ax, except for a 
likely smaller mass for the hydrogen envelope (0.6 M 0 , 1.3 M 0 , 2-5 M 0 respectively). Figures 
12.8112.9112.101 and 12.111 show a comparison of the spectrum of SN 2004ex and SN 1993J. Figure 
12.121 shows the comparison of spectra between SN 2004ex and SN 2008ax for all phases. In all 
phases He I lines are dominant and the H a line is weak. These features are compatible with a 
type lib supernova. 


2.6 Discussion 

The optical light-curves of SN 2004ex displays a temporal decline. From my results, it is not 
possible to extract information about the Ni mass, or the total energy of the ejecta. Indeed, the 
life-time of 56 Ni is around 8.8 days after the explosion; as the first spectrum is at 36 days after 
the explosion I cannot constrain 56 Ni. 

According to the classification of SN light-curves m, SN 2004ex has a rapidly declining 
light-curve which is a property of type lib SNe. And it is similar to the light-curves shapes of 
type lb SNe, suggesting that they have similar progenitors. 
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Figure 2.8: The spectrum at +36 days is shown with the identifications of the lines to compare 
the chemical composition of SN 2004ex and SN 1993J. 



4000 6000 8000 10 4 


X(A) 

Figure 2.9: The spectrum at +52 days is shown with the identifications of the lines to compare 
the chemical composition of SN 2004ex and SN 1993J. 


The comparison shows that the light-curve of SN 2004ex is very similar in shape to that of 
SN 2008ax, rather different that of the prototypical type lib SN 1993J. Indeed, SN 2004ex shows 
a faster decline rate at early phase after the peak and lacks the prominent narrow early-time 
peak which is seen in the light-curve of SN 1993J. Moreover, it shows a faster decline rate at late 
phase, similar to the light-curve of SN 2008ax. It is known that H-poor core-collapse supernovae 
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Figure 2.10: The spectrum at +60 days is shown with the identifications of the lines to compare 
the chemical composition of SN 2004ex and SN 1993J. 


ph 74(snl993J) - 



4000 6000 8000 


Figure 2.11: The spectrum at +86 days is shown with the identifications of the lines to compare 
the chemical composition of SN 2004ex and SN 1993J. 


display a wide range of behaviors in their light-curves evolution, as seen in the light-curves of 
all type lib SNe. Even if the light-curve of SN 2004ex is quite similar to the light-curve of 
SN 2008ax, it is significantly brighter than SN 2008ax and fainter than SN 1993J during the 
whole evolution. 

Type II SNe are generally characterized by their spectral properties rather than their photo- 
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Figure 2.12: The spectra of SN 2004ex at all phases are compared with those of SN 2008ax. 

metric properties. As it can be seen in all phases of SN 2004ex, the evolution of He features are 
prominent at late times. The nebular phase develops some months after the SN explodes, when 
the ejecta become transparent and the decrease of density allows for the formation of forbidden 
emission lines. 

A spectrogram of SN 2004ex obtained at early times shows a blue continuum with relatively 
broad, P-Cygni H lines similar to the spectra of young type II supernovae. To conclude, the 
comparison with SN 1993J and SN 2008ax shows that SN 2004ex strongly is a type lib. 

2.6.1 Possible progenitor 

The nebular phase of an SN is important to obtain information about the inner region of the 
explosion and its optically thin region. The modeling of this phase can provide information on 
the ejected mass, the kinetic energy, abundances, geometry (which can hardly be obtained by the 
other methods). Datasets for this phase were obtained for some type lib SNe which are listed in 
[139] ( e.g. SN 1994 and SN 2008ax). SN 2004ex can be added to the list with the results of this 
work. 

As it is discussed by Maurer et al. 1391 . H a lines at late observation times come from the 
interaction between the ejecta and the wind, it means that type lib SNe should be surrounded by 
massive stellar winds which can be from the companion star if the binary system is considered. 
As a result SN 2004ex could be a member of a binary system. 

The aim of the study in this Chapter [2] is to figure out the emission progress of SNe in the 
optical band and to compare with the different sub type of SNe. This could give the idea that 
if different kinds of SNe are connected to each other why mainly one type of SN (type Ic SN) 
observed to be associated to GRBs. 

This work is also encouraged to look for information about associated SNe to the GRBs which 
were found out during the sampling of low luminosity afterglow GRBs as seen in the following 
chapters. 
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Chapter 3 

Data Reduction Applied to X-ray 
observations of GRBs 


In this chapter, after introducing the Swift satellite and describing the data analysis techniques, 
the selection method of LLA GRBs is explained and discussed. A brief description of the sample 
follows. All errors are quoted at the 90% confidence level, and I used a standard flat ACDM 
model with fl m = 0.3, Ha = 0.7 and Hq = 70 km s -1 Mpc -1 . 

3.1 The Swift Satellite 

3.1.1 Detection techniques for X-rays 

In 1948, a system was proposed by Kirkpatrick and Baez to focus X-rays. It forms real images 
and consists of a set of two orthogonal parabolas, on which incident X-rays reflect successively, 
see Figure [TT~T1 > (left). This system satisfies the Abbe sine condition. It should be satisfied by 
every optical system to form clear images of an object at infinity, see Figure 13.1 I t. 



Principal Surface 



Figure 3.1: a) Schematic view of the Abbe sine condition. An optical system forms an image of 
an infinitely distant object as each ray in the parallel beaming emanates from its surface. The 
quantity / = h/sin{9 ) should be constant to form a clear image of the object. Here h is the 
radial distance of the ray from the optical axis and 6 is the angle between the final ray direction 
and initial ray direction, b) Kirkpatrick-Baez X-ray telescope simplified design, c) Wolter type 
I X-ray telescope design. Image Credit: NASA 


The right part of Figure 13.11 b shows the parallel mirrors (usually cylindrically symmetric) 
which increase the surface area. The most commonly used system is the Wolter Type I system, 
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which is represented in Figure 13.1b . It is a simple mechanical configuration and it provides the 
possibility of nesting several telescopes inside one another in order to increase the effective area. 

The X-ray mirrors of Wolter type I systems provide two things: 

• the ability to determine the location of the arrival of an X-ray photon in two dimensions, 

• simultaneously possessing a reasonable effective area. 

These instruments can be made of different materials, gold (the most common) or iridium 
which can reflect X-ray photons. For instance, the mirrors of the Swift/X RT are made of gold 
which is cooked on Nickel. 


3.1.2 Instrumental properties of Swift 

The Swift satellite made a “big bang” effect in the GRBs area, as it updated (and is still 
updating) all observational knowledge about GRBs. Since its launch on November 20 th , 2004, it 
has detected more than about 1000 GRBs and their associated afterglows. Additionally, it has 
allowed for multi-wavelength observations of both the prompt and the afterglow emission of the 
bursts in great detail. The description which follows is taken from }112j . 

The Swift satellite has three instruments: namely the Swift’s Burst Alert Telescope (BAT) 
observing in the gamma-ray band, the X-Ray Telescope (XRT), and the Ultraviolet/Optical 
Telescope (UVOT). They are shown on the spacecraft in Figure EO a [112) . BAT is dedicated 
to observing the prompt emission of GRBs and has an energy range of 15 — 150 keV. It has the 
capacity to detect weak bursts with its two-dimensional coded aperture mask and large area solid 
state detector array, as well as to detect bright bursts with its large field of view (1.4 steradians). 
BAT is able to locate a burst within an arcminute positional accuracy, allowing the satellite to 
point XRT and UVOT in the direction of the burst in ~100 s, in order to observe the afterglow. 
It provides spectra and light-curves at X-ray, ultraviolet and optical wavelengths allowing for a 
concurrent multi-wavelength examination of each burst. 

The Swift /XRT is a focusing X-ray telescope with a 110 cm 2 effective area, 23.6 x 23.6 arcmin 
field of view (FOV), 18 arcsec resolution (half-power diameter), and 0.2 - 10 keV energy range. 
The XRT can locate GRBs to 4 arcsec accuracy within 10 seconds of target acquisition for 
a typical GRB. It uses grazing incidence Wolter I mirror to focus X-rays onto a CCD. The 
XRT instrument is specifically designed to study X-ray counterparts of GRBs providing spectra 
and light-curves over a wide time range beginning 20 - 70 seconds after the burst trigger and 
continuing for days to weeks thence, covering more than seven orders of magnitude. The layout 
of the XRT is shown in Figure [221 b. 

The Swift /UVOT instrument has a unique capacity for afterglow studies. It is designed to 
capture the early UV afterglow, which is not possible from the ground because of the absorption 
by the atmosphere; as well as optical photons from the afterglow in the 170 — 600 nm band. 
Long term observations are usual. This instrument has a modified Ritchey-Chretien optical 
configuration with a 30 cm primary mirror which provides a good localization of the events with 
the error radius around 2 arcsec. It has UV/optical grisms and broadband color filters, a 4x 
magnifier; a clear white-light filter, and a blocking filter. 

An important improvement is provided by the quick ground-based follow up. When there is a 
possible GRB, Swift slews to the burst direction and quickly sendsthe position of the event to the 
ground to allow observations by other observatories. Position is relayed by the GRB Coordinate 
Network [ ITO] . 
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Figure 3.2: a: Schematic view of the Swift satellite with its three instruments, b: Overview of 
the XRT instrument. Image Credit: NASA/GSFC 


3.2 X-ray data analysis 

3.2.1 Extraction of light-curves and spectra 

For each source, I analyzed the XRT light-curve in flux unit, which was retrieved from the 
online Swift light-curve repositonQ [1411 [71], see one example in Figure |.'l.:!l However, I did 
not directly use these light-curves, as the flux calibration is too sensitive for an automated 
analysis. Because I cannot use directly the data downloaded from the online repository, and 
needed to estimate independently the spectral index and the count-to-flux conversion factor. I 
first applied, when needed, the latest available calibration to the data. First of all, the early 
XRT observation^, considering only photon-counting (PC) mode, were reprocessed at the given 
positior0 of each source using the xrtpipeline tool, which is part of the XRT software, distributed 
with the HEASOFT package]. 

Secondly, I analyzed the data using FTOOL^ taking the following steps. After the image 
extraction, I have chosen the source region (30 pixels radius) and the background region (60 
pixels radius), see Appendix C (see Figure [33] (left)). One difference with the detailed method 
followed in the Appendix [C] is the consideration of pile-up level. Pile-up occurs when multiple 
photons enter a pixel of the CCD within the same temporal frame (2.5 sec in the case of PC 
mode) and are read by the detector as a single photon of energy equal to the sum of all incident 
photon energies. It generally happens for early observations of bright sources. In this case, the 
source region has to be reduced by introducing an annulus to remove the region in which the 
count rate is larger than 0.6 counts.s -1 , see Figure [33] (right). 

However, since most sources in my sample are too weak to be effected by the pile-up, I simply 
removed data in that time interval when pile-up occurs instead of reducing the region of interest. 

1 http: / / www. swift .ac.uk/xrt .curves 

2 http: / / www. swift .ac.uk/swift .portal / 

3 http: / / www. swift .ac.uk/grb_region/ 

4 http://heasarc.gsfc.nasa.gov/heasoft/ 

5 http://heasarc.gsfc.nasa.gov/docs/software/ftools/ftools_menu.html 
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For that, after filtering the source region, I binned the light-curve by 100 or 500 seconds per bin, 
depending on the brightness of the source. Then, I extracted the light-curve on the full time 
interval and I restricted it to the time range for which there is no pile-up, i.e. when the count 
rate within a 30-pixel radius circle region centered on the source was below 0.6 counts.s -1 in 
PC mode [1421 . This limit is represented by the vertical line in Figure [331 while the horizontal 
line shows the time limit (which varies for each source). After that, I extracted the spectrum in 
the time interval free of pile-up. The background light-curve and spectrum are extracted from a 
region free of pile-up and prompt contamination m- 

I created an Ancillary Response File (ARF) using the full frame event lists in xrtmkarf and I 
used the Response Matrix File (RMF) which can be found in the Calibration Database0to fit the 
data (in Appendix O I show the full analysis of one GRB step by step with the corresponding 
explanation). 


Swifl/XRT data of GRB 081007 
blue: WT - red: PC 



Time since BAT Irigger (s) 


Figure 3.3: Example of a light-curve (here of GRB 081007) from the Swift repository. 


3.2.2 Fitting the spectra 

I fitted the extracted spectrum with Xspec using a power-law model, absorbed twice in order to 
take into account the host galaxy and the Galaxy contributions at low energy E < 2 keV. The 
redshift of the burst was set to the value obtained by optical spectroscopy and the galactic column 
densities ( Nn, ga i) was set to the value from the Leiden/Argentine/Bonn (LAB) Survey of Galactic 
HI |143] (which was calculated by the NH FTOOI0). However, the host column density ( Nh,x ) 
was fitted for. I generally used the \ 2 statistic when there were at least 20 — 25 counts per bin. A 
good fit for this statistic is obtained when 0.8 < y 2 < 1.1. Otherwise, I used the cash statistic for 
which at least one count per bin required. This statistic can provide parameters without giving 
information about the goodness of the fit. As an example, the fit result of the bright GRB 061007 
can be seen in Figure [331 Finally, using the fit, I calculated the energy correction factor of each 
GRB to reduce the energy range of observation from 0.3 — 10 keV to 2 — 10 keV. When this 
spectral model was found different from the one used in the automated pipeline, I corrected the 


6 CALDB; http://swift.gsfc.nasa.gov/docs/heasarc/caldb/swift/ 

7 http: //heasarc.gsfc.nasa.gov/lheasoft/ftools 
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Offset = 54015 10:44:27.5360 (SC time: 181910603.352004) 
Binsize = 100.000 s 

/Users/husnedereli/swift XRT/data strange GRBs/GRB061007/downlo 


rj T T 1 



% ... ! 


0 5000 10 4 1.5x10 4 2x10 4 


TIME s 


Figure 3.4: Light-curve of a bright burst (here, GRB 061007) during the data extraction. The 
horizontal line shows the threshold limit 0.6 counts.s , such that there is no pile-up. 



Figure 3.5: Left: The source and background regions are shown by the green circles. Right: The 
annulus in light green delimits the piled-up region. 


light-curve according to the parameters (spectral index and energy correction factor) from my 
analysis (for all these processes, one example can be seen in Appendix ICl). 
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data and folded model 



0.5 1 2 5 

Energy (keV) 


Figure 3.6: Example of an obtained spectrum (here, GRB 061007). The data and the best fit 
are represented together. 


3.3 Selection of the global sample and corrections of the 
light-curves 

3.3.1 Selection of long GRBs 

In order to define my sample, I first started with the list of bursts maintained by J. GreineiH. 
I took into account all bursts with a measured redshift observed before February 15 th , 2013, 
without consideration for the detector triggered by the event and/or observing it. Being inter¬ 
ested in the X-ray afterglow, I discarded all bursts not observed by an X-ray instrument (Bep- 
poSAX, Chandra, XMM-Newton, INTEGRAL, Swift). They are GRB 050408, GRB 051022, 
GRB 060712, GRB 061217, GRB 120716A. This leads to a first sample of 283 sources which 
contains long and short bursts. 

I then restricted it to long GRBs only. Searching all sources in the literature, I discarded 
12 short GRBs which have Tgo < 2 s. I also took into account cosmological time dilation and 
removed every intrinsically short burst with (Tgo, res t = 7go/(l + z) < 2 s) [21]. All discarded 
short GRBs are presented in Appendix [A] It leaves only 254 bursts in the sample; among these, 
61 were already investigated by [14]. GRBs of the global sample are presented in Appendix 
El with some of their parameters: redshift, energy correction factor (ECF), log(T a ), (where T a 
is the time at which the plateau ends) and spectral index. This work does not contain any 
update made in the Greiner webpage or in the Lancaster UK webpage more recent than the 
February 15 th , 2013. Moreover, GRBs with uncertain redshift measurements are not considered 
in this statistical study. Indeed, all redshifts were determined by spectroscopic observations in 
the optical band. 


ht t p: / / www. mp e. mpg. de/ ~ j eg/grbgen. ht ml 
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3.3.2 Selection of the late-time afterglow 

Once the flux calibrations have been checked (and when appropriately corrected, see section 
13.2.2D . I selected the late-time afterglow in the light-curves. Indeed, it was clearly stated that any 
prompt emission cause a broadening of the clustering, preventing any analysis m- I followed the 
same method, removing from the light-curves all emission present before the end of the plateau 
phase using the time T a at which the plateau ends. T a is obtained for each GRB by several 
methods described below. 

I mainly used the fit from the Swift webpag^f] when there is clear plateau and a break. One 
example is given in Figure 13.71 (left) for this condition. When no conclusion could be reached 


Swift/XRT data of GRB 060614 
blue: WT - red: PC 


Swift/XRT data ot GRB 050223 




Time since BAT trigger (s) 


5000 10 4 2x10 4 5x10 4 

Time since BAT trigger (s) 


Figure 3.7: Two examples for the determination of T a . Left: There is a clear plateau breaking 
at 44400 s. Right: when the break is not visible in the light-curve or there is not enough data 
to provide a meaningful fit, T a is imposed to be 10 4 seconds. 


with the light-curve obtained from XRT only, a comparison between the BAT and XRT light- 
curves0 was performed. When the decay measured by XRT corresponds to the decay of the 
prompt phase measured by BAT, this part of the light-curve up to the next change in the slope 
is considered to be associated to the prompt emission. Since the slope of the plateau is expected 
to be close to zero, the slope of the next segment in the light-curve determines if it is the plateau 
or not. If so, T a is taken as the time at which the next change in the slope happens. Conversely, 
an upper limit on T a is obtained. An example is given in Figure [TE1 

When no plateau can be identified, since I am interested only in the late afterglow emission, 
I have to identify segment 3 and eventually 4, see Figure 11.131 For that, two smoothly joined 
power laws (given in Eauation l3.14ll are fitted to the two last segments of the light-curves. If the 
slope of the very last segment is close to —2 and the slope of the previous segment is close to 

— 1.2, both parts of the light-curves are kept for the analysis and identified as standard afterglow 
and jetted afterglow respectively. Conversely, only the last part is considered. An example is 
given in Figure 13.91 and in Table 13.11 for GRB 090424. As the last decay index is found to be 

— 1.35, not compatible with the jetted afterglow, only the last part of the afterglow is considered 
in the analysis. 

Lastly, when the fit is not constrained (mainly for weak in X-ray emission), I imposed T a ~ 
10000 s, see one example in Figure [3?71 


9 http: / / www. swift .ac.uk/xrt_live_cat/ 

10 http: / / www. swift .ac.uk/burst .analyser / 
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BAT-XRT data of GRB 060912A 
BAT: Black — XRT: WT: Blue; PC: Red 



Figure 3.8: Left: Fit of the light-curve obtained from the Swift webpage. The break is at around 
600 s. Right: BAT and XRT light-curves together. The first observations by XRT seem to be 
associated to the decay of the prompt emission. After the break at 2000 s, the slope is not close 
to zero, thus T a is smaller than 2000 s. 



Figure 3.9: One example of an X-ray afterglow, fitted by smoothly joined power laws. As the 
decay index at late times is not compatible with a jetted afterglow, only the last segment of the 
light-curve is considered in the analysis. 


I finally removed all flaring emissions from the light-curves, which are thought to be part of 
the prompt emission [55k In addition, GRB 071122 was also removed from the sample because 
it was not observed enough. 
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Table 3.1: The fit result of GRB 090424. 


GRB z 

First 

index 

Second 

index 

x 2 

(reduced) 

d.o.f 

Tbreak 

(s) 

090424 0.544 

1.16T0.04 

1.35T0.03 

1.10 

311 

To 5 


3.3.3 Cosmological Scaling 

In all previous works on the clustering of X-ray afterglow, rather than using the luminosity, the 
authors used the flux at a common distance (set at z = 1 ) in order to reduce uncertainties on 
the Hubble constant and spectral parameters. Indeed, each measurement has its own associated 
error, and the k correction is very sensitive to the spectral index of the afterglow. Using a flux 
unit has the advantage of reducing the amount of correction applied, and therefore reduces the 
statistical scattering inserted in the data by the distance correction. To choose the faint sources 
from the global sample, I followed the same approach for the same reasons. 

Firstly, I applied the standard time dilation scaling: 


, tobs 

trest — , 5 

1 + ~0 

(3.1) 

^obSjZ—l — ^rest(l H - 1)? 

(3.2) 

j, _ ^ofcs(l + 1 ) 

tobs,z= 1 — i , > 

1 + z 0 

(3.3) 


where t res t is the time measured in the rest frame of the source, t 0 b s is the observed time and zq 
is the redshift of the burst. 

Using a flux unit has the advantage of reducing the amount of correction applied, and con¬ 
sequently this reduces the statistical scattering inserted in the data by the distance correction. 
The total luminosity: 


L = 4tt d\ F (3.4) 

= 4tt d 2 L z=1 F z = i, (3.5) 

where the last equality is obtained by imposing a constant observed luminosity (here L and F 
are the bolometric luminosity and flux with unit ergs.s -1 .cm~ 2 ). Then, the flux which would 
have been observed if the burst redshift was z = 1 is: 


F z= i = 


dl 


d l,z =1 


x F z . 


Finally, is the luminosity distance and it is computed by |144j : 


(3.6) 


d L = —{1 + z) 
no 


ri(a,n 0 ) = 2 \J s 3 + 1 


? 7 (l,Uo) -r? —— ,U C 

\1 + 2 j 

1 s s ^ s ^ 

-r - 0.1540^- + 0.4304— + 0.19097— + 0.066941s 4 
a 4 a 6 or a 


- 1/8 


s 


3 


1 — flo 
Uq 


(3.7) 


which is a fit of the exact integral expression. It gives a result precise at 4%. Finally, the 
parameters are FI o = 70 km.s _ 1 .Mpc _1 , FIm = 0.3. 
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3.3.4 Energy correction 

Then I restricted the energy range from 0.3 — 10 keV to 2 — 10 keV band. This allows to neglect 
absorption effects, so any correction is needed to take into account the absorption by the ISM. 
To do so, I used the late afterglow spectrum parameters (if not available I used the time averaged 
spectrum parameters): Nn^ai, Njj,int > z, photon index obtained from the automatic data anal- 
vsiJ^l I calculated the energy correction factor (ECF) by using a Heasarc tool: WebPIMM§F\ 
The conversion rate was taken to be one since a flux light-curve was used. Finally, the ECF is 
applied to the observed flux: 


Fcorr = Fobs x ECF (3.8) 

The k-correction takes into account cosmological correction and the down-shift of the gamma- 
ray energy from the burst to the observer’s reference frame since a higher energy component of 
the source spectrum is redshifted into the sensitivity band of the detector for high redshift sources 

PS] - It is: 


k{z) 


fH E N(E) dE 

J[1 dE1 


(3.9) 


where [ei, e^] is the sensitivity band of the instrument, N(E) is the source photon spectrum 
fitted with a simple power law for the afterglow of GRBs, see Eauation ll.2l Using Equations 1 1.2 1 
and 13.91 the k correction can be computed and is applied to all afterglow light-curves: 


(1 + zf 
~ ( 1 + z ) 2 ’ 


(3.10) 


where /3 is the spectral index, 
correction is: 


Rescaling all bursts at redshift one implies that the total 


k(zo) _ (1 + *o)^(l + l) 2 
"fcOT - (l + z 0 )2(l + l)/3’ 


(3.11) 


The k correction is very sensitive to the measured spectral index of the afterglow. As an 
example, with a redshift of 4 and a precision of 1.0 ± 0.3 for /?, the uncertainty on k is 5 ± 3, i.e. 
60%. Rescaling to z = 1 leads to 1.5 ± 0.7, i.e. an uncertainty of 16.8% : this method reduces 
the scattering induced by the uncertainties on the measurements, allowing for a more precise 
selection of the sample. 

After taking into account the k correction and distance correction, the flux at redshift z = 1 
is given by: 

F z ~! = — K z =i x F obs x ECF. (3.12) 

d L,z =1 


3.4 Selection of the LLA sample from the global one 

All corrected 254 light-curves are displayed in Figure l3T0l In 2008, [13] found a clear dichotomy 
in the flux of late afterglow light-curves considering 61 GRBs. The two groups were named group 
I (high luminous) and group II (middle luminous) and their properties determined. However, 
now the statistical significance of the complete sample is so high that several bursts can lie within 


11 http: //www.swift.ac.uk/xrt_spectra/ 

12 http: //heasarc. gsfc.nasa.gov/Tools/w3pimms.html 
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the separation zone of each group. In order to select only group III events (that are neither in 
group I nor in group II), I had to define the group I and II mean values and standard deviations. 
For this purpose, I found the flux values for the sources which have measurements before and 
after at one day after burst then I extracted the flux distribution at z=l at that time, and fitted 
it with two Gaussian distributions, assuming that each group is represented by a normal law. 
The Gaussian distributions are shown in Figure 13.Ill and their parameters summarized in Table 
13.21 Subtracting group I and group II leaves only the faint bursts. This way, I determined a flux 
threshold of 10~ 13 erg.s -1 .cm -2 (or 1.4<r from the mean value of group II, the 90% confidence 
level) at one day. 



io -3 io 2 io 1 io° io' io 2 io 3 

Time since burst (days) 


Figure 3.10: The light-curves of all 254 sources, corrected for distance effects (see text) and 
rescaled at a common redshift 2 = 1. The decay lines 1.2 and 1.4 are presented in black and blue 
respectively. 


However, not all bursts have a flux measurement at the time I performed the fit. I thus finally 
used the standard decay index a = 1.2, expected from the late time of the fireball evolution in 
a homogeneous medium with the power-law index of the accelerated electrons fixed to p = 2.3, 
to construct a time-varying upper limit for group III. In Figure 13.101 the black line is drawn 
by using the flux threshold of 10 -13 erg.s -1 .cm -2 at one day and a decay index of 1.2. When 
this limit was not enough to choose the sample, I used the mean decay index value of group II 
events (i.e. 1.4) which is presented in blue in the Figure 13.101 I finally discarded all bursts which 
significantly cross these limits, leaving only 31 LLA GRBs in sample. They are represented by 
blue diamonds in Figure [TT] and listed in Table PXTTl The bursts which are not in this sample are 
used as a control sample and they are represented with red points in Figure [TT] 


Table 3.2: Gaussian fit parameters of the flux distribution at one day for group I and II. 


Groups 

Mean flux 

Standard deviation 


(10 -12 c.g.i.) 

(log (10~ 12 c.g.i.)) 

I 

6.3 

0.40 

II 

0.6 

0.48 


As a result, to characterize a sample of GRBs with low luminosity X-ray afterglows, I selected 
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20 

18 



1e-15 1e-14 1e-13 1e-12 1e-11 1e-10 

Flux (ergs.crrf 2 .s _1 ,2-10 keV) 

Figure 3.11: The flux distribution at one day fitted by two Gaussians. 


a sample consisting of the 12% faintest X-ray afterglows from the total population of long GRBs 
with known redshift. 


3.5 The spectral index and temporal decay index calcula¬ 
tions 

The spectral and temporal decay indexes are obtained from the spectrum and the light-curve 
after T a , fitted by considering Eauation ll.221 To compute the spectral index, I followed the same 
method described in section 13. 2. II f extraction of the light-curve and the spectrum). Observations 
were added to the first one up to the point that a good statistic (x 2 or cash) was obtained. The 
results of the analysis are presented in Table [3731 with their % 2 or cash statistics. Moreover, the 
values of the column density Nh,x and the spectral indexes (which is equal to the photon index 
minus one) are summarized in columns 4 and 8 of Table BZD 

For the computation of the temporal decay index, the fit is performed either with a simple 
power law: 

f(x) = a x b (3.13) 


or a smoothly broken power law given by: 


f(x) = a x b exp(—(x/xb) n ) + c x d (1 — ea :p(—(x/xb) n )), 


(3.14) 
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Table 3.3: Spectral fit results of LLA GRBs. When the intrinsic column density Nh.x was found 
very small its value was set to zero and the fit performed again. The real values of Nh,x are 
given in parentheses on column 4. 


GRB 

z 

photon 

N H , x 

X^ 

d.o.f 



index 

10 22 cm -2 

(reduced) 


GRB 980425 

0.0085 





GRB 011121 

0.36 





GRB 031203 

0.105 





GRB 050126 

1.29 





GRB 050223 

0.5915 





GRB 050525 

0.606 

2.1T0.4 

0-38++ 

1.25 (cash) 

9 

GRB 050801 

1.38 

q 1 +0.39 
°- ± -0.36 

(0.0) 

0.81 (cash) 

108 

GRB 050826 

0.297 





GRB 051006 

1.059 

o c+0.44 
z -°-0.46 

(0.0) 

0.93 (cash) 

89 

GRB 051109B 

0.08 





GRB 051117B 

0.481 

2.23 

(0.0) 

1.38 

2 

GRB 060218 

0.0331 





GRB 060505 

0.089 

(1.8) 

(0.0) 



GRB 060614 

0.125 





GRB 060912A 

0.937 

1.63A0.18 

0.0 (0.06) 

1.07 

8 

GRB 061021 

0.3463 

2.02A0.06 

0.06±0.02 

0.93 

158 

GRB 061110A 

0.758 

1.42+6 

0.0 (9.51) 

0.88 

6 

GRB 061210 

0.4095 

(1.8) 

(0.0) 



GRB 070419A 

0.97 





GRB 071112C 

0.823 

1-83+3 

0.096+4 

0.66 

6 

GRB 081007 

0.5295 

1-99+4 

0-97++ 

0.213 

1 

GRB 090417B 

0.345 

2-34++ 

9 9q+0.31 

z.zo_ 0< 27 

0.78 

56 

GRB 090814A 

0.696 

T89++ 

0.0 (6.2e-10) 

0.9 (cash) 

35 

GRB 100316D 

0.059 

1 ^4+0-45 
i - O4 -0.47 

0.0 (0.84) 

1.25 (cash) 

44 

GRB 100418A 

0.6235 

1.87A0.26 

0.0 (2.6e-02) 

0.96 (cash) 

101 

GRB 101225A 

0.847 

(1.8) 

(0.0) 



GRB 110106B 

0.618 

9 q9+0.67 
z.oz_ 0 32 

(0.0) 

0.04 

2 

GRB 120422A 

0.283 

1.42±0.3 

0.0 (8.8e-07) 

1.08 (cash) 

85 

GRB 120714B 

0.3984 

9 ci +0.61 
z - J1 -0.66 

0.0 (0.0) 

1.24 (cash) 

14 

GRB 120722A 

0.9586 

1 17+4.34 

4*4 ' —1.98 

28.04+f 9 4 

0.29 

1 

GRB 120729A 

0.8 

1.79 ±0.23 

0.0 (4.16e-7) 

0.57 
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where x t, is the position of the break. Since I am interested in a sharp transition and as n is 
poorly constrained by the fit, it is set to the value n = 3. The results are presented on column 
7 of Table 14.11 and in Figure 13.121 Additionally, the y 2 and the degrees of freedom (d.o.f) are 
given in Table I7T~T1 


3.6 Conclusion 

X-ray instruments are important to work on the early and late afterglow of GRBs. Using the X- 
ray data, the environment of the bursts can be identified and the different emission mechanisms 
can be constrained. The X-ray afterglow light-curves show different properties: flares, injection 
of energy in the plateau phase and jet physics. Thanks to the increased number of observed 
GRBs in the past few years, statistical studies can now be performed. In the next chapter, I will 
present the statistical study which was performed on the properties of GRBs late afterglow. 
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Table 3.4: The temporal fit results of LLA GRBs. 


GRB 

Z 

First 

index 

Second 

index 

X 2 

(reduced) 

d.o.f 

GRB 980425 
GRB 011121 

0.0085 

0.36 


0.26A0.062 

0.0025 

2 

GRB 031203 

0.105 


0.42A0.037 

0.0004 

6 

GRB 050126 

1.29 

2.23A0.34 

0.94A0.09 

3 

0.32 

GRB 050223 
GRB 050525 

0.5915 

0.606 


0.91A0.025 

0.045 

1 

GRB 050801 

1.38 

l.liO.21 

1.25A0.13 

14 

1.73 

GRB 050826 

0.297 

0.31T0.19 

1.64A0.19 

0.27 

6 

GRB 051006 

1.059 

2.8±3.69 

1.69±0.13 

11 

2.88 

GRB 051109B 

0.08 

0.0 

0.96±0.1 

1.34 

13 

GRB 051117B 

0.481 

1.41T0.62 

1.03±0.50 

1 

1.51 

GRB 060218 

0.0331 


1.15A0.03 

0.75 

37 

GRB 060505 

0.089 


1.91±0.2 

2 

0.39 

GRB 060614 

0.125 

0.0 

1.73±0.05 

1.4 

147 

GRB 060912A 

0.937 

1 .02±0.1 

1.07A0.08 

1.15 

23 

GRB 061021 

0.3463 

1.38A0.33 

1.05A0.01 

0.94 

267 

GRB 061110A 

0.758 

1.79T0.55 

0.68A0.072 

1.12 

9 

GRB 061210 

0.4095 


1.67±0.85 

1 

0.81 

GRB 070419A 

0.97 


0.56±0.0 

0 

0 

GRB 071112C 

0.823 

2.86A0.49 

1.58A0.07 

0.57 

13 

GRB 081007 

0.5295 

0.47A0.31 

1.23±0.05 

42 

0.94 

GRB 090417B 

0.345 

1.62±5.7 

1.36±0.03 

1.03 

75 

GRB 090814A 

0.696 

2.3±0.27 

1.3±0.15 

0.72 

11 

GRB 100316D 

0.059 


1.42±0.06 

1.6e-28 

1 

GRB 100418A 

0.6235 

0.13A0.24 

1.42±0.08 

1.74 

22 

GRB 101225A 

0.847 

2.8±0.37 

3.05A0.87 

63 

2.05 

GRB 110106B 

0.618 

0.63A0.06 

1.35+/0.06 

52 

1.27 

GRB 120422A 

0.283 

0.27A0.16 

1.17A0.16 

0.98 

6 

GRB 120714B 

0.3984 


1.89±0.02 

0.18 

5 

GRB 120722A 

0.9586 

0.27±0.6 

1.99±0.66 

1.63 

1 

GRB 120729A 

0.8 

1.17T0.06 

2.8A0.23 

1.02 
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Figure 3.13: Continue of Figure 13.121 
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Chapter 4 


Low-Luminosity Afterglow GRBs 


In this chapter, I present the result of my study of the intrinsic properties of LLA GRBs (spectral 
index, decay index, distance, luminosity, isotropic radiated energy and peak energy), to assess if 
LLA GRBs and the other brighter GRBs belong to different populations. 


4.1 Properties of Low-Luminosity Afterglow GRBs 

The 31 LLA GRBs are represented by blue diamonds in Figure ITTl They represent about 12% 
of all bursts with known redshift. The other bursts, 223 GRBs, represented by red points, are 
used as a control sample. 



10" 3 10 2 10 1 10° io' io 2 10 a 

Time since burst (days) 


Figure 4.1: The afterglow light-curves of all GRBs, rescaled at a common redshift 2 = 1. LLA 
events are shown by blue diamonds while the control sample is shown with red dots. 


The LLA GRBs are also listed in Table [T~T1 which displays the GRB name, redshift, galactic 
and host absorption, Nh, galactic and host extinction, Ay , the afterglow temporal and spectral 
indexes, the isotropic and peak energies, and the time during which 90% of the energy of the 
prompt is emitted (Tgo). 
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Table 4.1: LLA GRBs and their main characteristics (see text). The spectral and temporal indexes of GRBs which happened before 


August 2006 are taken from D . 


GRB 

z 


N h 


Av 

Afterglow 

logT a 

E iso 

Ep,i 

Tgo 

Ref. 



Gal 

Host 

Gal 

Host 

Temporal 

Spectral 

(s) 

(10 52 erg) 

(keV) 

(s) 




(10 

21 cm- 2 ) 


(mag) 

index 

index 






GRB 980425 

0.0085 

0.428 


0.071 

1.73 

0.104=0.06 

(0.8) 


(1.3=t0.2)xl0 -4 

554=21 

18 

(1), (13) 

GRB 011121 

0.36 

0.951 


0.061 

0.38 

1.34=0.03 

(0.8) 


7.974=2.2 

10604=275 

28 

(1), (14), (15) 

GRB 031203 

0.105 

6.21 


0.117 

0.03 

0.54=0.1 

0.84=0.1 


(8.24=3.5)xl0 -3 

1584=51 

40 

(1), (14) 

GRB 050126 

1.29 

0.551 

(0.0) 

0.182 


1 -1+0.6 
-0.5 

0.74=0.7 


[0.4 - 3.5] 

>201 

24.8 

(17) 

GRB 050223 

0.5915 

0.729 

(0.0) 

0.078 

>2 

0.914=0.03 

1.44=0.7 


(8.84=4.4)xl0 -3 

1104=55 

22.5 

(2), (18) 

GRB 050525 

0.606 

0.907 

0.38]to.38 

0.221 

0.364=0.05 

1.44=0.1 

1.14=0.4 

3.8 

2.34=0.5 

1294=12.9 

8.8 

(5), (19) 

GRB 050801 

1.38 

0.698 

(0.0) 

0.989 

0.34=0.18 

1.254=0.13 

1 84"*" 0-56 

-L.Oft_o.53 

3.2 

[0.27 - 0.74] 

<145 

19.4 

(5), (17) 

GRB 050826 

0.297 

2.17 

st 3 

2.398 


1.134=0.04 

1.14=0.4 

4.04 

[0.023 - 0.249] 

>37 

35.5 

(17) 

GRB 051006 

1.059 

0.925 

(0.0) 

2.345 


1.694=0.13 

1 r+0.44 
- L -°-0.46 

2.77 

[0.9 - 4.3] 

>193 

34.8 

(17) 

GRB 051109B 

0.08 

1.3 

<2 

0.3 


1.14=0.3 

0.7 4=0.4 

3.14 



14.3 


GRB 051117B 

0.481 

0.46 

(0.0) 

0.321 

< 1.4 

1.034=0.5 

(0.8) 


[0.034 - 0.044] 

<136 

9.0 

(11), (17) 

GRB 060218 

0.0331 

1.14 

6±2 

0.437 

0.54=0.3 

1 S'*" 11 
±-O -0.6 

0.514=0.05 

5.0 

(5.44=0.54) xlO -3 

4.94=0.49 

~2100 

(1), (20) 

GRB 060505 

0.089 

0.175 

(0.0) 

0.209 

0.634=0.01 

1.914=0.2 

(0.8) 


(3.94=0.9)xl0 -3 

1204=12 

~4 

(1), (21) 

GRB 060614 

0.125 

0.313 

0.5±0.4 

0.068 

0.11±0.03 

2 Q+ 0 - 3 

0.84=0.2 

4.64 

0.224=0.09 

554=45 

108.7 

(3), (21) 

GRB 060912A 

0.937 

0.420 

(0.0) 

1.436 

0.54=0.3 

1.014=0.06 

0.64=0.2 

3.3 

[0.80 - 1.42] 

>211 

5.0 

(4), (17) 

GRB 061021 

0.3463 

0.452 

0.6±0.2 

0.185 

< 0.10 

0.97±0.05 

1.024=0.06 

3.63 


46.2 

(3) 

GRB 061110A 

0.758 

0.494 

(0.0) 

<0.10 

< 0.10 

1.14=0.2 

0.44=0.7 

3.68 

[0.35 - 0.97] 

>145 

40.7 

(3), (17) 

GRB 061210 

0.4095 

0.339 

(0.0) 

0.489 


1.674=0.85 

(0.8) 


[0.10 - 0.33] 

>105 

85.3 

(17) 

GRB 070419A 

0.97 

0.24 

< 10 

0.081 

<0.8 

0.564=0.0 

(0.8) 


[0.20 - 0.87] 

<69 

115.6 

(5), (17) 

GRB 071112C 

0.823 

0.852 

<5 

0.203 

0 20' 1 " 0 ' 05 
u - zu _ 0 .04 

1.434=0.05 

n q+0-5 

U -°-0.4 

3.0 



15 

(4), (17) 

GRB 081007 

0.5295 

0.143 

n 97+6.9 
' -0.97 

0.196 

0 sfi+ 006 

u.oo_o 04 

1.234=0.05 

n 99+ 0 - 88 

0.43 

4.5 

0.184=0.02 

614=15 

10 

(7), (22) 

GRB 090417B 

0.345 

0.14 

22±3 

0.083 

0.84=0.1 

1.444=0.07 

1.34=0.2 

3.54 

[0.17 - 0.35] 

>70 

>260 

(6), (17) 

GRB 090814A 

0.696 

0.461 

(0.0) 

0.15 

<0.2 

1.04=0.2 

(0.8) 

3.5 

[0.21 - 0.58] 

<114 

80 

(7), (17) 

GRB 100316D 

0.059 

0.82 

(0.0) 

0.088 

2.6 

1.344=0.07 

0.54=0.5 


(6.94=1.7)xl0 -3 

20±10 

292.8 

(8), (23) 

GRB 100418A 

0.6235 

0.584 

(0.0) 

0.623 

0.0 

1.424=0.09 

0.94=0.3 

4.82 

[0.06 - 0.15] 

<50 

7.0 

(9) (17) 

GRB 101225A 

0.847 

0.928 

(0.0) 

0.311 

0.75 


(0.8) 

4.65 

[0.68 - 1.2] 

<98 

1088 

(12), (17) 

GRB 110106B 

0.618 

0.23 

(0.0) 

0.032 


1.354=0.06 

1 oo + 0.67 

4.04 

0.734=0.07 

1944=56 

24.8 

(24) 

GRB 120422A 

0.283 

0.372 

(0.0) 

1.241 

0.0 

1.34=0.3 

0.44=0.4 

5.07 

[0.016 - 0.032] 

<72 

5.35 

(10), (25) 

GRB 120714B 

0.3984 

0.187 

(0.0) 

0.077 


1.894=0.02 

(0.8) 


0.08±0.02 

694=43 

159 

(17) 

GRB 120722A 

0.9586 

0.298 

350+??g 

0.555 


1.24=0.4 

1.24=1.2 


[0.51 - 1.22] 

<88 

42.4 

(17) 

GRB 120729A 

0.8 

1.4 

(0.0) 

0.112 


2.84=0.2 

0.8 4=0.3 

3.9 

[0.80 - 2.0 ] 

>160 

71.5 

(17) 


Note: for A v values: (1) P3Sj; (2) 1 J7|: (3) [SB]; (4) 1 if) : (5) SO : (6) [ST ; (7) [S3; (8) PS3-; (9) [SI]; (10) [Si ; (11) [S51 : 


(12) 157 ; for E iso & E Pji values: (13) |158| : (14) [159] : (15) [160] : (17) in this work; (18) [161] : (19) [162] : (20) [163] : (21) [164] : (22) 

[THH] : (23) [S3; (24) [SS]; (25) [S3 


O 

CO 
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4.2 Statistical Properties 

4.2.1 The redshift distribution 

The redshift distributions of LLA GRBs is shown in Figure 14.21 together with that of all lGRBs 
in the global sample. From a rapid examination of Figure 14.21 it is seen that LLA GRBs are 
closer than normal lGRBs, whose mean value is z ~ 2.2 ( e.g. [1681 1169 ) ). Table \L2 \ displays the 
statistical parameters of the two distributions. A Kolmogorov-Smirnov test is performed on the 
two datasets which shows that the probability for the two distributions to be based on the same 
population is 1.1 x 10 -14 . 



Redshift 



Redshift 


Figure 4.2: Left: Redshift distribution of LLA GRBs (blue) compared to the normal lGRBs 
(red). Right: Cumulative distribution of the same samples. 


Table 4.2: Statistical parameters of the cumulative redshift distributions. 


Parameter 

LLA GRBs 

All lGRBs 

mean 

0.5454 

2.1785 

median 

0.5295 

1.98 

standard deviation 

0.3746 

1.187 


The difference between the redshift distribution of LLA GRBs and that of lGRBs can be due 
to a selection bias, or alternatively, it is an intrinsic property. Faint events are more difficult to 
detect than brighter ones. In section 14.41 I will discuss the prompt properties of LLA GRBs. 
However, these events are barely more energetic than the detection threshold, see e.g. m, and 
the measure of the redshift is based on spectroscopic observations at optical wavelengths. As the 
afterglow decays with time, the probability to obtain a successful measurement decreases as well. 
If the afterglow is faint from the very beginning of the event, then it is initially difficult to measure 
a redshift. Several of these bursts could be distant GRBs with a low luminosity afterglow and 
one hypothesis is that the large difference between the two distributions is due to distant events. 
In order to investigate this hypothesis, the distribution of LLA GRBs and of normal lGRBs are 
restricted to bursts with redshifts z < 1. At this distance, all bursts are detected independently 
of their intrinsic luminosity. I have recomputed the cumulative redshift distributions for this 
sub-sample (see Figure fOl) . The difference is still large, and from a Kolmogorov-Smirnov (K-S) 
test the probability that the two distributions are drawn from the same population is 9.4 x 10~ 4 . 
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I can therefore conclude that the observed LLA GRBs are nearby events because of their low 
intrinsic luminosity. 



Figure 4.3: Cumulative distribution of the redshift of LLA GRBs (blue) compared to that of all 
normal lGRBs (red) with redshifts z < 1. 


4.2.2 Absorption and Extinction 

I constrained absorption (gas) and extinction (dust) on the line of sight of LLA GRBs. The 
X-ray absorption has little (if none) effect on my selection criteria since I used the flux in the 
2.0 - 10.0 keV band, where absorption can be neglected m- The intrinsic hydrogen column 
density Nh can be linked to the host properties m ■ In addition, the optical extinction can bias 
the distribution (for instance the well-known problem of dark bursts, e.g. (1731). It additionally 
decreases the probability of measuring the redshift. 

Milky Way Galaxy 

For consistency, I checked that the distribution for the Milky Way Galaxy values of Ay,Gai and 
A r H,Gai (i.e. the optical extinction and X-ray absorption parameters) is consistent with the 
whole sample of normal lGRBs. The Nh,Go.i was calculated as explained in section 13 . 2.21 The 
optical extinction was calculated using the NASA/IPAC extragalactic database^ for the Landolt 
V-band measured by (136 . for all bursts but except GRB 060904B and GRB 061110A. These 


1 http: / / ned.ipac.caltech.edu/forms / calculator.html 
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two bursts are seen in projection on the galactic disk ( i.e. with a galactic latitude between —5 
and 5 degrees), where the measures of [136] are highly variable with the position. For these two 
events, I will rely on the most accurate measurements of m and m respectively. The results 
are reported in Figure 14.41 and Figure 14.51 As the two samples are consistent for Nh, ga i only, I 
will consider in the following that the absorption in the Milky Way Galaxy has not introduced 
a bias in the LLA GRBs sample. 



-1 - 0.5 0 0.5 1 

iog N H Qa| (10 21 cm" 2 ) 

Figure 4.4: Galactic HI column density distribution for LLA GRBs sample (blue) and normal 
long GRBs (red). 


Host galaxy 

I obtained the Nh,x values for the intrinsic host absorption from the X-ray data analysis of all 
LLA GRBs by using Xspec. They are displayed on Figure 14.61 Most of them are compatible 
with little or no intrinsic absorption. Comparing Figure 14.41 and Figure 14.61 I can see that for 
the bursts with a non-zero Nh,x, the absorption of the host galaxy is on average a factor 10 
larger than in the Milky Way Galaxy, as already noted by m ■ At low redshift this effect was 
attributed to the gas in the host galaxy. 

Optical extinction values are displayed on Figure 14.71 The extinction for most sources can 
be fitted by the Small Magellanic Cloud (SMC) model, with the exception of GRB 060912A for 
which the Milky Way (MW) model was applied. I cannot see strong differences with the typical 
values reported for the Milky Way Galaxy, suggesting that bursts are not located in a dusty 
medium. 
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Figure 4.5: Histogram of the optical extinction Ay in the Milky Way Galaxy for LLA GRBs 
sample. GRB-SN associations are represented by cyan bars while other LLA GRBs are shown 
with blue bars. 


4.3 Afterglow properties 

The late time afterglow emission corresponds to the decaying part after the plateau phase in the 
canonical X-ray afterglow light-curve EZj, see also segments III and IV of Figure 11.131 In this 
regime the flux F v is modelized by 


F v oc t a v /3 , (4.1) 

where a and /3 are respectively the temporal and the spectral indexes. 

4.3.1 Temporal decay index 

The distribution of the temporal decay indexes for LLA GRBs is presented in Figure ITT51 In 
order to have a reference sample for comparison, I used a reference sample of bursts (higher 
luminous) listed in P~5| and not members of the LLA GRB subclass. Note that the decay index 
of GRB 060607A reported in this last article is incorrect and is not considered in the comparison. 
As it can clearly be seen, the two samples seem similar. This is confirmed by a K-S test (p = 
0.88), that confirm the hypothesis that the two samples are drawn from same populations of 
bursts. 
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Figure 4.6: Host galaxy HI column density, N^.x■, distribution for LLA GRBs sample. The blue 
bars are fitted values while the white ones with red arrows are upper limits. 


4.3.2 Spectral index 

The distribution of the spectral indexes for LLA GRBs is presented in the left panel of Figure 
HU I used again the reference sample (higher luminous ones 114] ) for comparison. Here, the two 
distributions seem to be different and it is confirmed by a K-S test: the hypothesis of a single 
population is confirmed with a probability of 0.084. 

4.3.3 Closure relations in the X-ray band 

I used the closure relations [Mi muni Hnra, see also section [1. 4. 21 to investigate the burst 
geometry, the fireball micro-physics, its cooling state and the surrounding medium. They are 
displayed on Figure 14.101 together with the spectral vs. decay indexes of each burst. Clearly, 
within the measurement errors, most bursts can be explained by at least one of the closure 
relations, with two exceptions. 

GRB 120729A 

The non-jet closure relations are rejected for this event. In fact, this burst can be interpreted 
with jet in the slow-cooling phase, with the X-ray band located between the injection and cooling 
frequencies (be. v m < vxrt < ty), with an electron distribution index of p = 2.8 ± 0.2. This 
hypothesis is strengthened by the presence of a break in the light-curve at t(, =8.1 ks. Assuming 
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Figure 4.7: Histogram of the host galaxy optical extinction Ay for LLA GRBs. The green bars 
represent measured values of GRB-SN associations. On the other hand, the blue bars show 
measured values and the white bars with red arrows display upper and lower limits (depending 
on the direction of the arrow) burst with no SN detected. 


that the afterglow part located before this break is the true standard afterglow, I recomputed 
the temporal and decay indexes: the new values agree with the non-jet closure relations (green 
point in Figure IXIUl) . 

GRB 060614 

This event could be compatible with a jetted afterglow, with p = 2.25 ± 0.05. However, the 
errors bars are large enough to accommodate some non-jetted closure relations. Thus, no firm 
conclusion can be reached on the jet hypothesis for this source based on the closure relations 
alone. 


4.3.4 Closure relations in the optical band 

The temporal decay and spectral indexes of 7 LLA GRBs at optical wavelength (see Table 14.51) 
are calculated and they are presented along with closure relations of v a < v < and one of 
v m < v < v c , on Figure 14.111 The optical component of GRB 060614 can be explained by a 
jetted afterglow with v m < ^optical < v c , which strengthen the finding in the X-ray band. Finally, 
the five of LLA GRBs with optical measurements indicate a fireball in the slow-cooling regime 
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Figure 4.8: Distribution of the decay indexes of LLA GRBs (blue) compared to the reference 
sample (red). The decay index reported in (T3] for GRB 060607A is incorrect, and not shown on 
this figure. 




Figure 4.9: Left: Distribution of the spectral indexes of LLA GRBs (blue) compared to the 
reference sample (red). Right: Cumulative distributions of the same samples. 
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Figure 4.10: X-ray decay index versus spectral index of LLA GRBs. The purple filled circle 
and square represent GRB 120729A and GRB 060614 respectively. The green dot represent 
GRB 120729A before the break at = 8.1 ks. All closure relations, indicated by the lines, are 
computed for p > 2 in the slow-cooling phase. Solid and dash-dotted lines stand for v m < v < v c 
and v c < v respectively. Blue, red and black lines stand for ISM, wind medium, and jet effect 
respectively. As it can be seen from Table 11.21 the fast-cooling regime for ISM and wind v > 
u rn > v c has the same dependence with slow-cooling regimes of v m < v c < v. 


with v m < voptical < v c, which is expanding in a ISM. GRB 011121 is compatible with the slow 
cooled fireball expanding in a wind medium, compatible with constraints from the X-ray bands. 

GRB 060614 also showed the same properties in the jetted closure relation like seen in X-ray 
data v m < v < v c . Five of the GRBs with optical measurements indicated that v m < v < v c in 
the slow-cooling region of ISM. GRB 011121 indicated the same region but wind closure relation 
which is compatible with X-ray emission. 


4.4 Prompt properties 

The Swift /BAT and the Ferrm-GBM datasets are re-analyzed for the LLA GRBs listed in Table 
14.11 They are combined with previously published results from Konus-Wind and BeppoSAX. 
About half of the events have a firm measurement of the prompt parameters (E p ^ and Ei SO ), 
the other half present upper and lower limits. However, prompt parameters of some sources are 
missing in the table. Concerning GRB 051109B, it was only detection by Swift /BAT and its 
spectrum is fitted by a simple power-law, which does not allow to constrain E p ^. Concerning 
GRB 071112C, any values could not be derived because the dataset is incomplete. The prompt 
properties of GRB 061021 is given by Em] obtained from the Konus-Wind data for the first 
pulse, E Pt i is 1046 ± 485 and Ei SO is 0.46 ± 0.08 x 10 52 ergs, while it has also a weak tail. So this 
source was not taken to account in the study of E p ^-Ei SO correlation. 

During the calculations for the lower and upper limits for E p , the BAT spectrum was fitted 
by the power-law (PL) model. When the PL photon index is lower than —1.7, the low-energy 
index /3 is set to —2.3 and the fit is performed with Band function in order to constrain the 
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Figure 4.11: Optical decay indexes versus spectral indexes of LLA GRBs. Blue, red and black 
colors stand for ISM, wind medium, and jet effect respectively. Solid lines stand for slow-cooling 
regime with v m < v < u c . Big circles stand for slow-cooling regime with v a < v < v m . Finally, 
small circles stand for fast-cooling regime with v a < v < v m . In green optical decay vs. spectral 
indexes of LLA GRBs. The purple filled square represents GRB 060614. The cyan filled circle 
represents GRB 011121. 


Table 4.3: The temporal decay and spectral indexes of 7 LLA GRBs in the optical band. 


GRB 

name 

Redshift 

Spectral 

index 

Decay 

index 

GRB 011121 

0.36 

0.8±0.15 

1.9T0.1 

GRB 050525A 

0.606 

0.97T0.1 

1.3 

GRB 050801 

1.38 

1±0.16 

1.4T0.05 

GRB 060614 

0.125 

0.47T0.04 

1.9T0.01 

GRB 060912A 

0.937 

0.62 

0.93±0.04 

GRB 070419A 

0.97 

0.8 

1.3±0.03 

GRB 071112C 

0.823 

0.63±0.29 

0.9T0.02 
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E iso (erg) 


Figure 4.12: Location in the E pi — Ei SO plane of LLA GRBs sample compared to both short 
and normal lGRBs. The best fitting power law, E Pi i = a x E l ’ so , obtained by fitting the data 
accounting the sample variance. Derivation of E pp from the best fitting power law of ±2<r log E p ^ 
assuming that a log E p p is the central value of the 90 percent confidence interval, they are 
computed by using the Reichart method as reported in Ha- 


other parameters. This gives a lower limit on E p j. On the other hand, when the PL photon 
index is between —1.7 and —1.9, the situation is uncertain and any reliable limit on E p cannot 
be derived. However, in the opposite situation, when the PL photon index is larger than —2 the 
low-energy index a is set to -1, and the fit is performed with Band function in order to constrain 
the normalization, j3 and E p . This gives an upper limit of E p ^. Both upper and lower limits are 
at 90% confidence level and are computed using the error procedure in XSPEC. 

In order to estimate Ei SO and its uncertainty, first Ei so is computed by assuming the best fit 
PL. Then, Ei SO is computed by performing the fit using the Band function 11.11 with a , /3 and E p 
set at the values that they had when computing the lower or upper limit to E Pt i. In this way, 
a range of lower and higher value of Ei SO is obtained for each spectrum. Ei SO is taken as the 
average of lower and higher limits while the error is the higher value of Ei so minus the average. 

It is found that the E p .i values cluster broadly within the 40-200 keV range, as expected. 
Figure 14.141 shows the histogram of the E p ^ values for both LLA GRBs and normal lGRBs. 
However, the situation is different with Ei SO , There is a clear shift in the Ei SO axis with respect 
to most normal lGRBs. LLA GRBs are less energetic during their prompt phase compared 
to normal GRBs. Figure 14.141 shows the histogram of Ei SO values for both LLA GRBs and 
normal lGRBs. By taking into account the median redshift of LLA GRBs, the high energy limit 
of BAT, the E pp values, and the intrinsic scatter of the Arnati relation, measurements up to 
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Figure 4.13: Location in the E p ^ — Ei SO plane of LLA GRBs sample; on the left, compared to 
short events with the non-compatible lower limits on E p ^ and Ei SO ; on the right, same comparison 
with firm measurements and non-compatible lower limits on E p y only. 




1000 1500 2000 2500 3000 3500 

E .(keV) 


(10 52 ergs) 


Figure 4.14: Histogram of the E p ^ (left) and Ei SO (right) of LLA GRBs (blue) and of normal 
IGRBs (red). When only a lower and upper limits were given, the avarage values are considered. 


Eiso = 3 x 10 53 ergs would be expected, almost one order of magnitude larger than the BAT 
measurements listed in Table 14.11 Thus, I conclude that this effect is not due to a bias, but 
is an evidence that LLA GRBs are intrinsically less energetic, both during the prompt and the 
afterglow phases, compared to normal IGRBs. 


Amati Relation 

The peak energy (E Pt i) in the rest frame of the burst of the vFv spectrum is given: 

E Pt i = Epochs x (1 + z ), (4-2) 

where 2 is redshift, -E Pj0 bs is the observed peak energy: 


E Pi0 bs — (2 + o) x E 0 , 


(4.3) 
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where a is the low-energy index and Eq is the break energy. E p ^ 0 b s is the energy at which most 
of the power of the burst is emitted [55] . The isotropic equivalent energy ( Ei SO ) can be obtained 
by taking the integral of the spectrum and by using the redsliift. It is given by: 


F,- = 


47r ^l 

1 + 2 


^bol 7 


where cLl is the luminosity distance which is given by ESI: 

dz' 


Al = 


k^i 


1 0 \/flm(l + + flfc(l + 2 / ) 2 + 

Finally, Sboi is the bolometric luminosity computed by: 


(4.4) 


(4.5) 


„ „ J*EN(E)dE 

^bol ^obs /1 I _. \ i 

/e+li) EN (E)dE 


(4.6) 


where [E\ = 1 keV, E 2 = 10000 keV] is the bolometric gamma ray range, [ei, e 2 ] is the sensitivity 
band of the detector, N(E) is the source photon spectrum which is fitted by the Band function, 
see Equation 11.11 (551 . 

E pi and E iso are correlated: E pp oc E^„ ■ This was discovered in 2002 by using 11 BeppoSAX 
bursts m- This correlation was later confirmed and extended to X-ray rich GRBs (XRRs) and 
X-ray flashes (XRFs) based on HETE-2 data 1,7b . 

This correlation is used to constrain prompt emission mechanism, jet geometry and properties 
but also unify GRB and XRF etc. (see 179i for a review). 

The main implications of the correlation include prompt emission mechanisms, jet geometry 
and properties, GRB/XRF unification models, identification and nature of sub-classes of GRBs 
( e.g. sub-energetic, short) (see [179 1 for a review). I also used this relation to confirm that LLA 
GRBs are different from normal lGRBs when considering isotropy equivalent energy. 


4.4.1 Jetted afterglows 

GRB 120729A can be accounted for by the closure relation of a jet. The opening angle is given 
by m who extended the work of |l02j to account for the radiative efficiency of the prompt 
phase ? 7 7 : 

*(**• E ‘“ ) = °'i 61 (if;) " 1/8 (s&) ■ (4 - 7) 

In my computations, the standard values for the number density of the medium n = 1 cm' 3 and 
the radiative efficiency ?y 7 = 0.2 are used. I got 9 = 2.7°. This is indeed a jet-break as it is seen 
at the same time in both the X-ray and optical bands [ T5T)1 451 . 

GRB 060614 is compatible with a jet, according to the closure relations. This burst also 
displays an achromatic break (around 36.6 ks) in X-ray, optical and UV [182] . Before the jet- 
break, this burst features a plateau phase and not a standard afterglow. If I assume the presence 
of a jet, the corresponding jet opening angle is 6.3°. 

A statistical study by m shows that the mean opening angle of lGRBs is 9 = 4.7°. The 
results for both GRBs are consistent with this value, and the jet, when detected, is not different 
from that of normal lGRBs. 
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4.4.2 Environment of the bursts 

Most sources can be explained by either a wind environment or a constant ISM. As shown 
below, many of these sources are associated with SNe (see Table EB- This association would 
point towards a wind environment |184] . However, as shown by [TS5j . the termination shock can 
lie very close to the star, and I cannot conclude firmly on the surrounding medium. The absence 
of strong dust extinction also does not for a strong wind. 

One source deserves a more careful study: GRB 120729A. From the jet part of the light- 
curve, I derived p = 2.8 ± 0.2. This value is not compatible with both the spectral and temporal 
decay indexes (1.08±0.03 and 0.74±0.072 respectively) of the pre-break part of the light-curve. 
Only the spectral index is marginally consistent with this value of p , assuming v m < vxrt < 
and a constant ISM. The temporal decay is too flat (while expecting a value of at least 1.5). In 
order to reconcile all of these facts, I need to involve some late-time energy injection to flatten 
the light-curve m- This energy injection needs to be present during the pre-break part of the 
light-curve, but should stop during the post-break part. Let us note that the sampling of the 
X-ray light-curve is not good during the jet break and allows for some non-simultaneity. 


4.4.3 Microphysics of the fireball 

For LLA GRBs (see section 14751) . when a closure relation indicates a given status of the fireball, 
the X-ray band is located below the cooling frequency. This fact is not common. Indeed, most 
late GRB afterglows are compatible with the X-ray band located above the cooling frequency 
iTB7irraBi , see section IP1 

In the case of a constant ISM, the formula of the cooling frequency is given by Equation 1 1.1 91 
m and is reminded here: 


v c = 3.7 x 10 14 -K7j 1 / 2 n -1 (F + 1 )- 2 e~ 3 /_ 2 2 T~ ll2 Rz, (4.8) 

where E 33 is the isotropic energy in units of 10 53 ergs, n is the number density of the medium 
in cm -3 , Y is the Compton parameter, e b ,-2 is the fraction of internal energy in the magnetic 
field, Td is the observed time expressed in days after the burst. 

Instead in the case of a wind medium, the cooling frequency is given by [I06j: 

= 3.5 x 10 14 i4 / 2 A- 2 (y + l)- 2 e- 3 / 2 T] / 2 Hz, (4.9) 

where A* is a constant normalizing number density of the wind. 

I start by assuming that the fireball expands in the constant ISM. The X-ray band is from 
0.3 to 10.0 keV, which correspond to 7.2 x 10 16 Hz and 2.4 x 10 18 Hz respectively. I assume that 
v c is above 3.7 x 10 18 Hz for simplicity. Equation 14.81 simplifies to : 

10- 4 £- 1 / 2 e- 3 Z 2 2 < 1, (4.10) 

when assuming the standard density n = 1 cm -3 , the Compton parameter Y <C 1 and considering 
the observation time Td = 1. From the prompt parameters, I find that E 53 can be as low as 
10 -5 for GRB 980424, I then finally obtain : € 3,-2 > 0.1. A larger value of E 53 reduces the 
constraint on e 3 . In any case, a constraint of £b ,-2 > 0.1 is not really constraining, as typical 

values of € 3—2 should be in the order of 1 . I therefore conclude that, under the hypothesis of 

the fireball expanding in a constant ISM, the uncommon position of the cooling frequency is due 
to the small energy of the fireball. 

The situation is similar when assuming a wind medium, for which Equation 14.91 implies : 


10“ 4 E 


1/2 - 3/2 

53 e B -2 


< 1 , 


(4.11) 
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when assuming standard density and Compton parameter 7 < 1. and considering the obser¬ 
vation time Td = 1. From the prompt parameters, I find that £53 can be as high as 0.797 for 
GRB 011121, I thus finally obtain : e b ,-2 > 0.0012, which is again not constraining. A smaller 
value of £53 reduce the constraint on cb- 


4.4.4 Prompt properties of LLA GRBs 

In Figure 14.121 it is clearly seen that all outliers to the Arnati relation are LLA GRBs. Several 
explanations have been proposed to explain these events (see m and reference therein for 
details): GRB 060505 may be a short GRB (as its location in the £ Pi * - £j SO plane in Figure 
14.131 mav suggest); the E p l value of GRB 061021 refers to the first hard pulse, while a soft 
tail is present in this burst (so the true E p ^ may be lower); GRB 031203 may be much softer 
than measured by INTEGRAL/ISGRI as supported by dust echo measured by XMM. I notice, 
however, that the outliers are all located on the left part of the diagram, I will propose two 
alternative solutions. 


The Amati relation is valid for LLA events 

The left part of the E p ^ - Ei SO plane in Figure lT.lSI relates to soft and faint events. In this part of 
the diagram, the usual gamma-ray instruments are not well-suited to measure the prompt prop¬ 
erties. For instance, if only BAT measurements of GRB 060218 [489] are taken into account, this 
events is more similar to GRB 980425, i.e. a clear outlier. The reason can be lack of the softer 
part of emission due to limited energy band of the detector. However, the advantage of Swift 
compared to BeppoSAX is its ability to perform simultaneously X-ray and gamma-ray observa¬ 
tions. Combining the XRT and BAT measurements allows to make this event (GRB 060218) 
fully consistent with the Amati relation. 


The Amati relation is not valid for LLA events 

On the other hand, GRB 980425, GRB 060505 and GRB 050826 are not compatible at all with 
the Amati relation. Assuming that these measurements are correct, then the best-fit relation in 
the £ P) i - Ei SO plane changes dramatically, being way more flatter. A flatter Amati relation has 
been foreseen as early as 2003 m, using GRBs seen largely off-axis. This is not in contradiction 
with the results of the closure relations, as these relations apply only to events seen on-axis. For 
completeness, I also note that a similar explanation hold in case of the canonball model m- 
Being seen off-axis, these events are expected to be less luminous than normal lGRBs even during 
the afterglow (see e.g. l92j ). Thus, the discrepancies between the LLA GRBs and normal lGRBs 
are explained by the geometry of the events and not by the progenitor. 


4.5 Conclusions 

I present strong evidence that LLA GRBs belong to a population of nearby events. They represent 
about 12% of the total population of long bursts, and their afterglow emission are faint. It is 
also found that the afterglow spectral and temporal properties are similar to those of the general 
population. I show that there is no more selection effects due to absorption and extinction 
introduced by the Milky-Way and host galaxies than for normal lGRBs. 

However, these events are also faint during their prompt phase, and located in another area 
of the E Pt i~Ei SO plane in Figure 14.121 They also include all outliers of the Amati relation. In 
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addition, their redshift distribution is different from that of normal lGRBs. Indeed there is a lack 
of luminous lGRBs at small distances, which should have been detected, since they are luminous. 

The similar properties of LLA bursts and normal lGRBs point towards the same emission 
mechanism for the afterglow and might indicate that LLA GRBs are the tail of the luminosity 
distribution of normal long GRBs. However, it does not explain the lack of luminous GRBs in 
the local universe. This tension might be solved by further GRBs monitoring and the observation 
of luminous nearby sources. Instead this tension might increase with the observation of several 
additional low-luminous nearby sources, which might be the indication for a different kind of 
bursts. 
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Chapter 5 


Additional Properties of LLA 
GRBs 

5.1 Supernovae associated to LLA GRBs 


Nine LLA GRBs are firmly associated to SNe by spectral and photometric optical observations. 
There are also two other sources (GRB 070419A, GRB 100418A) that may be associated to SNe. 
The first one shows a faint bump in its light-curve similar to the one of GRB 980425 m while 
for the second one, a bright host galaxy may prevent to observe the signature of a faint SN, 
which cannot be brighter than —17.2, comparable to the magnitude of the faintest Ic SN [194] . 
However, two sources (GRB 060505, GRB 060614) are firmly not associated to SNe (while they 
are nearby). All referenced GRB-SN pairs are listed in Table IQ1 It also contains the pairs for 
which the burst is not an LLA GRB. As it can be seen for several well-known associations, the 
burst is an LLA GRBs. For example, GRB 980425/SN 1998bw, GRB 031203/SN 20031w, and 
GRB 060218/SN 2006aj. If the former was thought to be due to a hypernovae [l95| , the latter 
has been proposed to be a neutron star experiencing an SN [1961 . Being nearby objects, it is 
worth wondering why not all LLA GRBs are associated to SNe. The most trivial solution is a 
large optical extinction preventing the detection of the SN. However, when looking at Figure 14771 
which displays the distribution of host galaxy extinction Ay, the non-associated events are not 
more affected by a larger extinction than the associated ones. The same is true for the Galactic 
extinction, see Figure H~51 Other observational problems, such as a bright host, can mask out 
a faint underlying SN signal. However, a large fraction of LLA GRBs are associated with SNe 
(64%). This can lead to a potential problem: if LLA GRBs are intrinsically different (as seen in 
the E Pt i - Ei SO plane in Figure 14.121) from normal bursts, the SN-GRB connection might not be 
representative of all lGRBs. 

In general, types of SNe which are associated to GRBs are Ic, broad-line Ic and Ib/c, as seen in 
column 5 of Table I57T1 except SN 2001ke associated to GRB 011121 which is of type Iln: hydrogen 
can be find in its spectrum. This GRB and the other LLA GRBs have different properties. It 
is the brightest one (7.97 ± 2.2 x 10 52 ergs) with the highest peak energy (1060±275 keV). It 
is accommodated by the closure relation of a wind in both X-ray band (see Figure 14.101) and 
optical band (see Figure ITTT1) . 
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Table 5.1: GRB-SN events with some physical parameters. The last part of the table shows the 
firm associations of bursts which are not LLA GRBs. The isotropic luminosity Li SO is calculated 
from Ei SO by dividing it by Tgo, see Table 14.11 


GRB 

name 

redshift 

SN 

identification 

SN 

name 

SN 

type 

Host 

type 

LLA 

GRB 

Liso 

(GRB) 

10 49 erg.s -1 

Ref. 

GRB 980425 

0.0085 

spectral 

SN 1998bw 

BL-lc 

dwarf spiral 
(SbcD) 

yes 

0.033 

(1) 

GRB 011121 

0.36 

spectral 

SN 2001ke 

Iln 

no 

yes 

387 

(2) 

GRB 031203 

0.105 

spectral 

SN 20031w 

BL-lc 

Irr 

Wolf-Rayet 

yes 

0.56 

(3), (4) 

GRB 050525 

0.606 

spectral 

SN 2005nc 

~Ic 

no 

yes 

417 

(5), (6) 

GRB 060218 

0.0331 

spectral 

SN 2006aj 

BL-Ib/c 

dwarf Irr 

yes 

0.02 

(7), (8), 

(9) 

GRB 081007 

0.5295 

bump 

SN 2008hw 

Ic 

no 

yes 

30 

(10), (11) 

GRB 100316D 

0.059 

spectral 

SN 2010bh 

BL-lc 

Spiral blue 

yes 

0.056 

(12), (13) 

GRB 120422A 

0.283 

spectral 

SN 2012bz 

Ib/c 

no 

yes 

0.44 

(14), (15) 

GRB 120714B 

0.3984 

spectral 

SN 2012eb 

I 

no 

yes 

0.7 

(16), (17) 

GRB 070419A 

0.97 

bump 

no 


no 

yes 

6.6 

(18) 

GRB 100418A 

0.6235 

no detection 



dwarf blue 

yes 

24.4 

(19) 

GRB 060505 

0.089 

no association 



Spiral(SbcD) 

yes 

0.82 


GRB 060614 

0.125 

no association 



no 

yes 

2.2 

(21), 

(22), (23) 

GRB 021211 

1.01 

spectral 

SN 20021t 

~Ic 

no 

no 

969 

(25), (26) 

GRB 030329 

0.168 

spectral 

SN 2003dh 

BL-lc 

no 

no 

75 

(27) , 

(28) , (29) 

GRB 091127 

0.49 

bump 

SN 2009nz 

BL-lc 

no 

no 

345 

(30), (31) 

GRB 101219B 

0.55 

spectral 

SN 2010ma 

Ic 

no 

no 

29 

(32), (33) 

GRB 130215 

0.597 

spectral 

SN 2013ez 

Ic 

no 

no 

75 

(34), (35) 


Note: for GRB-SN associations: (1) 1741: C2i 11971; (3) [1981: (4) 11991: (5)J2jp; ffit IMT1: (7) Eg]; (8) 11631: f9l 

[2g];(10) gig]; (11) ggg]; (12) ggl]; (13) gig]; (14) gjg]; (15) gig]; (16) gig]; (17) gni]; (18) 11931 : (19) gnl; (20) 

[ 2 l 2 l : ( 21 ) [2l3l : ( 22 ) [2111 : (23) I 2 T 5 I: (24) l 2 l 6 l: [251 [2T71: ( 26 ) [2T51: (27) 12191: (281 l 220 l : (29) [T5] ; (30) I22T1 : (31) 

152^1 : (32) 12251 : (33) 12211 : (34) [2251 : (35) [2251 : 


5.2 The signature of an SN in the GRB light-curve 

For the GRB-SN associations, most SNe are identified by spectroscopic analysis as seen in column 
3 of Table 1-5.11 On the other hand, there is another way to detect an SN in GRB liglit-curves, 
which is to identify a bump in the optical afterglow light-curve. This bump is expected between 
10 to 20 days after the GRB explosion. For such an identification, the light-curve of GRB 980425 
is considered as a template. 

I used this template by transforming the U-band light-curve of GRB 980425 at redshift 0.0085 
to the redshift of GRB 060218 (z = 0.0331), to show the SN bump in its light-curve. This light- 
curve is compared in Figure 15.11 (left) with the one of GRB 060218. This association is also 
confirmed by spectral observations. 

To perform the transformation, hrst of all, I used the classical conversion of apparent magni¬ 
tude mi at redshift z\ to the apparent magnitude m 2 at redshift Zi'- 

ma = mi - 5 log >° (snSj) ’ (!U) 

where cLl is the luminosity distance and it is given by Equation 14.51 Finally, z\ and Z 2 are the 
redshifts of the GRB which used as a template and of the tested GRB respectively. 

Then, the time correction is applied: 


O — ti 


1 + Z\ 

l + z 2 ’ 


(5.2) 










5.3. LLA GRBS WITHOUT SNE IN THE SAMPLE 


79 


where t\ and O are times measured at redshift z\ and Z 2 respectively. 

However, the bump of the SN (assuming it exist) cannot be seen in most GRB light-curves 
which might be because of the bright afterglow. This is shown by GRB 061007 for example, see 
Figure [5J] (right). The template for the SN bump is transformed at the redshift GRB 061007 
and compared to each optical band (BVRI). 



"finaltUJ.b 
“grb061007 . 
TinaltR_1.txt 
"grb061007_R 

"grb06100O 

"grt>061007 " 


Figure 5.1: Left: U-band observations of GRB 980425 transposed to the redshift of GRB 060218. 
Right: The R and U-band observations of GRB 980425 transposed to the redshift of GRB 061007. 
They are compared with the multi-band observations of GRB 061007. 


There is also another way to figure out a possible bump and determine the time of the SN in 
the extrapolated light-curve of GRB without optical observation. To do it, I used the light-curve 
of GRB 090618, which shows an SN bump ten days after the explosion without a spectroscopic 
confirmation |227] . It is applied to GRB 970828 228 . I transformed the Rc-band light-curve 
of GRB 090618 at the redshift of GRB 970828 following the same method. This light-curve is 
compared with the limits which were given by |16| (purple) and 117) (blue) on Figure 15.21 the 
possible SN bump would have been visible from ~ 20 to ~ 40 days after the burst triggered 
by neglecting local absorption while optical observations were present for up to 7 days from the 
GRB triggers, reaching a R band magnitude limit ~ 23. 8 m and subsequent deeper images 
after ~ 60 days IT7I . 


5.3 LLA GRBs without SNe in the sample 

LLA GRBs are weak enough during the afterglow phase (by definition) to allow the detection of 
any coincident SN. However, only 9 GRBs out of 31 LLA GRBs are firmly associated to an SN. 
Several reasons can be given to explain this discrepancy. First, optical follow up observations 
may be in sufficient during the time at which SN appear depend on the weather conditions, 
position of the bursts and the capacity of the instruments. Second, the contribution of the host 
galaxy can be important, and allow the signal of a weak supernova to be hidden. For instance, 
the host galaxy of GRB 090417B is quite bright. Additionally, GRB 061210 is in a galaxy cluster 
which even prevents the exact identification of its host. On the other hand, I showed that even 
for bright bursts, the afterglow is often decaying fast enough to allow for the SN identification. 

The third argument is the distance. Among all GRB-SN pairs, the most distant one (GRB 021211) 
is at redshift z = 1.1 urn, while the most distant pair (GRB 050525) for LLA GRB is at com¬ 
parable redshift z = 0.606. Evidence of GRB-SN associations should also be found at larger 
redshift. However, the spectroscopic identification of SNe becomes difficult at large redshift be¬ 
cause they are faint. Additionally, optical contaminations from the host galaxy and the afterglow 
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Figure 5.2: The transposed Rc-band liglit-curve of GRB 090618 at the redshift of GRB 970828 
(green) compared with Rc-band (red) light-curve of GRB 090618. The purple and blue lines 
represent the limit given in the deep images for GRB 970828 by jT6j and m respectively. 


of GRBs reduce the signal-to-noise ratio. Thus, detection of SNe associated to GRBs are not 
expected to be successful above z ~ 0.7. Discarding all sources at larger redshift (10 LLA GRBs) 
and those that are associated to SNe leaves 12 LLA GRBs. Additionally, GRB 090417B (bright 
galaxy) and GRB 061210 which is in a galaxy cluster are removed to leave only 10 GRBs. 

Interestingly, three very nearby bursts, GRB 060505 |229] . GRB 060614 [MJ and GRB 100418A 
Em, are not associated to SNe. They are in high metallicity host galaxies (i.e. evolved ones) 
and were already described as long (Tgo) short GRBs, because no supernovae were found. 

To conclude, many LLA GRBs are firmly associated to SN (64% of all firm associations), 
the remaining being too far to give a firm conclusion. However, three nearby sources are not 
associated to SNe, despite their proximity, which questions their classification in the long GRB 
category, or questions the assumption that all long GRBs are associated to SNe. More observa¬ 
tions of nearby LLA GRBs are required to give a conclusion on the GRB-SN association, and 
might even lead to the creation of a new class of events: nearby, low-luminous and not associated 
to a supernova. 


5.4 Host galaxies of LLA GRBs 

The host of long GRBs are generally low metallicity, low mass and high star formation rate 
galaxies [230 '. 1231 . 232j . In this section, the differences between the host galaxies of LLA GRBs 
and those of normal lGRBs are studied. 

5.4.1 Mass and metallicity 

First, the mass and metallicity of host galaxies are studied. Figure [53] displays the position of the 
bursts (both LLA GRBs and normal lGRBs) in the mass-redshift plane. No large difference can 
be found between the two samples. Indeed, the mean value of the host masses is 10 9 ' 5±o ’ 32 Mq 
for normal lGRBs, while it is slightly smaller for LLA GRBs, the mean value being 10 9 ' 1±0 27 Mq. 
However, the masses of the host galaxies of GRBs are in average smaller than those of the galaxies 









5.4. HOST GALAXIES OF LLA GRBS 


81 


in the Sloan Digital Sky Survey (SDSS) sample, whose average mass is represented by the green 
line on Figure IQ1 (right ) (232: . 




Figure 5.3: Left: The mass-redshift distributions of LLA GRBs (blue filled circles) and lGRBs 
(red filled circles) host galaxies. The solid lines represent the mean value of the metallicity 
12 + log (O / H)kko4 = 8.7 (black) of SDSS galaxies [18^ and the mean value of the metallicity 
12 + log (0/L7)kko4 = 8.4 (blue) of LLA GRB host. Right: Same mass-redshift distribution 
rescaled to redhifts from 0 to 0.8. The green solid line represent the average binned mass of 
SDSS galaxies. 


In order to discuss the metallicity, a relation between the redshift, the mass of the galaxy 
and the metallicity is given: 


12 -(- log (O/LOkkcm = - 7.5903 + 2.5315 log M* - 0.09649 log 2 M* 

+ 5.1733 log t# — 0.3944 log 2 1# (5.3) 

— 0.4030 log tn log M*, 

where KK04 represents the metallicity scale of Kobulnicky & Kewley (2004) [ 233 ] . tn is the 
Hubble time at redshift z in Gyr and M* is the galactic stellar mass in unit of solar mass. The 
Hubble time at redshift z is given by: 


t H (z) = 


dz / 


Hi 


o Jz (1 + zd)\J H m (l + z/)^ Ha 


(5.4) 


Cosmological parameters are considered, fl m = 0.27, LIa = 0.73 and H 0 =71 km s 1 Mpc 1 in 
order to be consistent with usa. 

The mean value of the metallicity is close to 12 + log {0/H)kko4 = 8.7 [18], while the 
mean value for the host of LLA GRBs is 12 + log (O/H) kko4 = 8.4, to be compared to 12 + 
log (0/H)kko4 = 8.35 for the host galaxies of normal lGRBs. Thus, no significant difference 
between LLA GRBs and normal lGRBs is found for the metallicity. 

However, only a small fraction of LLA GRBs (17%) and of 1GRB (30%), are below the line 
representing a metallicity 12 + log (0/H)kko4 = 8.4 through Equation 15.31 obtained by fitting 
all SDSS galaxies. Interestingly, it implies that the host galaxies of GRBs have larger masses 
than those in the SDSS survey at a constant metallicity level, and that the normal 1GRB host 
masses are slightly larger than the LLA GRB host masses at a constant metallicity level. 
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It was proposed that the progenitor of long GRBs be massive metal-poor stars, which are 
rare in massive galaxies. The creation of such stars requires particular conditions, and especially 
low metal abundance in the galaxy (or at least in the star formation region). 

5.4.2 Metallicity and brightness 

Obtaining the magnitude of the host galaxies of GRBs in the optical band is extremely challeng¬ 
ing. Thus, the information is available only for a few bursts in my sample. Figure [5~fl represents 
the metallicity of the host galaxy as a function of its magnitude for LLA GRBs. 
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Figure 5.4: The distribution of metallicity-B magnitude LLA GRBs. The barred line present the 
cut-off limit of metallicity at 12 + log(O/H) KD02 = 8.5 [Hi. GRBs associated to broad-line Ic 
SNe are presented by black points while GRB 011121 which is associated to type Iln SN 2001ke 
and GRB 050826 (no association to SN) are presented by red and blue points respectively. 


Additionally, all bursts represented by the black dots are associated to broad-line Ic SN. A 
cut-off in metallicity was obtained to differentiate the SNe Ic associated to GRBs from those 
that are not 12 + log(O/H) KK04 = 8.5 [hi. Below this metallicity, all SNe are associated to 
a GRB, while above none of them is seen with a GRB. Interestingly in my sample, one burst 
(GRB 011121) is associated to SN and its host metallicity is larger than 12+log (O/H) KK04 = 8.5. 
However, this SN is of type Iln (hydrogen is found in its spectrum), which indicates a smaller 
initial main sequence star, around 20 Mg. Moreover, GRB 011121 is the most luminous LLA 
GRB. Thus, the energy budget required to produce both the SN and the GRB might not be 
fulfilled from a small star, indicating that the progenitor might be in a binary system, with the 
GRB and the SN produced by two different objects. 

In addition, GRB 050826 is also above 12 + log(O/H) KK04 = 8.5, however it is not associated 
to an SN, due to the lack of observations. 

Finally, I did not compare the hosts galaxy properties LLA GRBs with those of normal 
lGRBs, as the information required is difficult to find. This is one of the open question that I 
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wish to continue for a future project. 


5.4.3 Star formation rate and metallicity 

Figure [53] displays the star formation rate (SFR) against the metallicity for the hosts of normal 
and LLA GRBs. The mean values for the star formation rate are 2.94MQ.yr _1 for LLA GRBs 
and 6.29MQ.yr -1 for normal lGRBs [232] , The reason for this difference is that since less stars 
are formed as a whole, less very massive stars are formed as well. Thus, the possible progenitor 
of a GRB would be less massive, and as a result the GRB less energetic. 



Figure 5.5: The SFR-metallicity distribution of LLA GRBs (blue circles) and lGRBs (red circles) 
host galaxies. The black solid line represents the binned metallicity of SDSS galaxies in metallicity 
range Ba¬ 


it is also clearly seen that the SFR of GRBs hosts is larger compared to that of SDSS galaxies, 
as found in a previous study 15321 . 


5.5 Rate density of LLA GRBs 

5.5.1 Observed rate density 

Within 31 LLA GRBs, 28 GRBs were observed by the Swift satellite. The three remaining bursts, 
GRB 980425, GRB 011121 and GRB 031203, respectively observed by BATSE, INTEGRAL and 
Ulysses, before the Swift launch, are excluded from this study since each instrument has its own 
properties: flux threshold, redshift measurement efficiency, operation time of satellites, etc. 

Except one source GRB 101225A, called the “Christmas burst”, all other sources have a firm 
measured peak flux, see column 3 of Table 331 I excluded GRB 101225A of this computation 
even if it has a peak flux measurement after 1372 s [234 . 
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The observed rate density of one event is given by [535] : 

RGRB,obs = T7 -- 7fi (5-5) 

where S cov is the fractional sky coverage (0.17 for the Swift satellite), T is the time during which 
the satellite was observing. Here, it is the difference between the time of the first observation by 
the Swift and the time of the latest burst in the sample, i.e. 8 years. However, when I take into 
account the times when the Swift is not observing because of the South Atlantic Anomaly (~ % 
10 of all observations), satellite slewing and technical problems (~ %10 of all observations) this 
time is 6.4 years. So the telescope is actually performing the % 80 of all observations [2861 . 

Finally, z max is the maximum redsliift such that the fluence of a burst would be above the 
threshold detectability. Correspondingly, V max is the maximum volume, is given by m : 

pZmax 

V max = / — dz, (5.6) 

Jo d z 

where dV/dz is a volume element factor, which is given by [237] : 


dV ^ 4ttc dL 2 (z) 

dz H 0 (1 + z) 2 h(z) ’ 1 j 

where dL(z) is the luminosity distance, calculated by assuming a flat-A universe (Ok = 1 — Qm — 
Ha) with H 0 = 72 km.s -1 .Mpc -1 , Om = 0.3 and Ha = 0.7. For this precise computation, the 
fit given by Equation 13.71 is considered as It prevents the computing of a double integral. 

The h(z) factor appearing in Equation 15.71 is the normalized Hubble parameter: 


h(z) = H{z)/H 0 = [H m (1 + z? + H a ] 1/2 . (5.8) 

The maximum redshift z max appearing in Equation 15.61 corresponds to the maximum lu¬ 
minosity distance d max such that the burst would have been detected by satellites. z max is 
calculated using the luminosity distance given by Equation 14.51 and d max is given by m-- 


dm.n.x. — 


F p k(z) 


F l im k(z L 

im) 


dL(z), 


(5.9) 


where F p is the observed peak flux (given in the third column of Table 15151 for each burst), 
k(z) and k(zLim) the k corrections for a burst at the measured and maximum limiting distances 
respectively. Fj jirn = 0.4 ph.s _1 .cm~ 2 is the theoretical Swift /BAT sensitivity. 

I used the largest redshift value in my sample (z = 1.38 for GRB 050801 and detection peak 
flux is 1.46 ±0.14 ph.s -1 .cm -2 ) for zum because of the large uncertainty on the luminosity 
function of the sample. Actually, z max should be consistently computed by finding the solution 
to: 


dmax(z rnax ) 


F p k(z) 


FLim k(z rnax ) 


d L (z). 


(5.10) 


However, the computation is done using znm , since the error is expected to be small in any 
case. The first reason is that if z ma x was larger than ZLim > additional bursts should be observed. 
The second reason is that since LLA GRBs are intrinsically sub-luminous and they are close to 
the instrumental threshold, it would be difficult to observe them at a distance much larger than 

%Lim' 
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The k correction, which accounts for the downshift of 7 -ray energy from the burst to the 
observer’s reference frame, is computed using Equation 13.91 considering [ei,e 2 ] is 15 - 150 keV, 
which is the sensitivity band of Swift /BAT. Additionally, N(E ), the source photon spectrum is 
assumed to be given by the Band function, see Equation 11.11 using the average GRB spectral 
indexes a = —1, /3 = —2.3 and a rest frame power-law break energy at 511 keV. 

The rate density of each burst is given in column 5 of Table 15.21 The total observed rate 
density, obtained by summation of each single rate, is 0.37 Gpc” 3 .yr _1 . It does not take into 
account the two sources GRB 060218 (63.6 Gpc _ 3 .yr _1 ) and GRB 100316D (43.7 Gpc _ 3 .yr _1 ), 
since their peak flux is below the instrumental threshold as seen in the lower part of Table 15.21 
On the other hand, they have compatible values with the result of 15371 . when assuming to 
be 0.15 ph.s _ 1 .cm~ 2 . Moreover, from the summation of each GRB, two more GRBs are excluded: 
GRB 051109B (1.69 Gpc _ 3 .yr _1 ) and GRB 060505 (0.17 Gpc _ 3 .yr _1 ) because of their huge rate 
density which are higher than all others. The reason can be their small redshifts. 


Table 5.2: The rate density results of LLA GRBs sample. GRBs with their correspond¬ 
ing redshift (1, 2 respectively), photon peak Flux between 15 — 150 keV: BAT 1 - sec Peak 
(http://swift.gsfc.nasa.gov/archive/grb_table/) (3), isotropic luminosity (4), observed GRB rate 
(5), and intrinsic GRB rate (beaming angle taken into account) ( 6 ). 


GRB 

z 

Peak flux 

Liso 

R-GRB,obs 

R-GRB ,int 

name 


ph.s - 1 .cm —2 

ergs.s -1 

Gpc —3 .yr —1 

Gpc- 3 .yr— 1 

GRB 050126 

1.29 

0.71±0.17 

1.7x10 s1 

0.002379 

2.143 

GRB 050223 

0.5915 

0.69±0.16 

8.8xlO so 

0.00975325 

8.78575 

GRB 050525 

0.606 

41.704=0.94 

4.2 X10 s1 

0.0005637 

0.5077875 

GRB 050801 

1.38 

1.464=0.14 

6.2 X10 50 

0.00132513 

1.1936875 

GRB 050826 

0.297 

0.384=0.13 

4.9X10 49 

0.0858775 

77.35875 

GRB 051006 

1.059 

1.624=0.30 

1.5x10 s1 

0.00178275 

1.605875 

GRB 051109B 

0.08 

0.554=0.13 


1.693375 

1525.375 

GRB 051117B 

0.481 

0.494=0.14 

6.4x 10 49 

0.02108125 

18.99 

GRB 060505 

0.089 

2.654=0.63 

8.2 x10 48 

0.1696875 

152.8625 

GRB 060614 

0.125 

11.504=0.74 

2.2 X10 49 

0.01536125 

13.8375 

GRB 060912 

0.937 

8.584=0.44 

4.3x10 s1 

0.00077015 

0.69375 

GRB 061021 

0.3463 

6.114=027 


0.00393413 

3.543875 

GRB 061110 

0.758 

0.634=0.12 

2.8x10 s0 

0.00652538 

5.878 

GRB 061210 

0.4095 

5.314=0.47 

3.6X10 49 

0.00330775 

2.979625 

GRB 070419 

0.97 

0.204=0.10 

6.6X10 49 

0.0115675 

10.42 

GRB 071112C 

0.823 

84=1 


0.00092668 

0.83475 

GRB 081007 

0.5295 

2.64=0.4 

3.0x10 s0 

0.00377013 

3.396125 

GRB 090417B 

0.345 

0.34=0.1 

1.2 X10 49 

0.0783825 

70.6075 

GRB 090814 

0.696 

0.64=0.2 

8.4x 10 49 

0.00803413 

7.237125 

GRB 100418 

0.6235 

1.04=0.2 

2.4xlO so 

0.00628113 

5.658 

GRB 110106B 

0.618 

2.14=0.3 

4.8xlO so 

0.0034325 

3.092 

GRB 120422A 

0.283 

0.64=0.2 

4.4x 10 48 

0.0574575 

51.75875 

GRB 120714B 

0.3984 

0.44=0.1 

7.0X10 48 

0.04009 

36.1125 

GRB 120722A 

0.9586 

1.04=0.3 

4.0xl0 so 

0.00294663 

2.65425 

GRB 120729A 

0.8 

2.94=0.2 

1.3X10 50 

0.0018005 

1.621875 


GRB 060218 

0.0331 

0.254=0.11 

2.6xl0 4b 

63.64375 

45864 

GRB 100316D 

0.059 

0.14=0.0 

5.6X10 40 

43.73875 

39400 


5.5.2 Estimation of upper limit density 

The intrinsic rate density is obtained by multiplying Equation 15.51 with a term which takes into 
account that if GRBs are collimated, the number of observed bursts is reduced as some of them 
will not have their jet aligned with the detection. Then, Equation 15.51 becomes: 

RGRB,int. = 77 7 ; TP , 

Vmax &cov -*■ 


(5.11) 









86 


CHAPTER 5. ADDITIONAL PROPERTIES 


where B{9) = [1 — cos(0)] _1 and 6 is the beaming angle. In the LLA GRB sample, two 
sources have been shown to be collimated in jet, with a corresponding beaming angle 2.7° for 
GRB 120729A and 6.3° for GRB 060614. The smallest value (2.7°) is used in the computation 
of B(9) to give an upper limit to 9 for all bursts with an unknown beaming angle. 

The results are presented in the last column of the Table 15.21 The total intrinsic rate density 
obtained is 330.9 Gpc _3 .yr _1 following the same way as in subsection l5.5.1l 


5.5.3 Redshift effect on the rate density 

In order to consider all bursts without redshift, I should take into account the efficiency ( r ] z ) of 
obtained redshift, as applied by 2371 . So the Equation 15.51 becomes: 


Rgrb 


1 111 

I'm. ax Scov 7 Qz 


(5.12) 


where r) z is 0.3 for lGRBs observed by the Swift satellite. The total observed rate density is 
1.23 Gpc _3 .yr _1 . When the beaming of GRBs outflows is taken into account, the total intrinsic 
rate density is 1103.04 Gpc _3 .yr _1 . However, the value of r] z = 0.3 is biased towards low- 
luminosity bursts. Indeed, when the optical instruments are fast enough and the sources are 
luminous enough, the redshift measurements are highly efficient [ 238( . Conversely, if the sources 
are less luminous, such efficiency decreases, as it is difficult to get a spectrum of quality. As the 
bursts considered here are intrinsically less luminous (both during the prompt and the afterglow 
phase), it is expected that the redshift measurements be on average less efficient than for all 
other lGRBs. 


5.5.4 Evolution of the rate density with the redshift 

In the previous sections, I assumed that the rate density is constant with redshift. However, it is 
unlikely to be the case at high redshift. One way to deal with the problem would be to consider 
the observed redshift distribution of GRBs but it is known to be biased since GRBs are thought 
to be linked to the death of massive star. Another solution is to the consider star formation rate 
and to re-normalize it. In this case Rgrb(z) is given by: 


Rgrb(z) = p 0 


Bsf(z) 
Rsf{z = 0) 


Rsf(z) 


0.02 + 0.122 

1 + (z/3.23) 4 - 66 ’ 


(5.13) 


where Rsf(z = 0) = 0.02 is the local star formation rate, po is the re-normalization constant, 
interpreted as the local rate density of GRBs. 

For normal lGRBs, it was found by [ 2571 that po = 0.09 ± 0.01 Gpc Ayr 1 obtained by the 
best fit of the differential peak flux distribution of Swift. Indeed, taking into account the star 
formation rate, the number of bursts per unit of time, redshift and luminosity, is given by 2571 : 


dN dV(z) R G rb{z ) 

itdTT L = — TTT dz m 


(5.14) 


where </>(L) is the luminosity distribution function. Integrating over redshift and luminosity gives 
the rate at which bursts are observed normalized by po , which is found by a direct comparison to 
the observed value. As the number of sources considered in [237: is a few hundreds, the authors 
could self-consistently determine po and the luminosity distribution function of lGRBs. 
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However, as the LLA GRBs sample only contains 31 bursts, the computation cannot be 
reproduced since the luminosity distribution function of LLA GRBs cannot be constrained. 

5.5.5 Summary 

The total intrinsic and observational rate density of the LLA GRBs are presented in Table 15.31 
The effect of opening angle (made more than 100 times) and redshift efficiency (made more than 
3 times) are not negligible. 


Table 5.3: The total rate densities in each conditions. The observed and intrinsic total rate 
density are computed by using the parameters of 23 LLA GRBs. They are presented on the first 
two columns; the last two show the results of computation made by taking rj z into account. 


R-GRB,obs 

FLGRB.int 

RGRB,obs 

RGRB ,int 

Gpc 3 .yr 1 

Gpc 3 .yr 1 

Gpc 3 .yr 1 

Gpc 3 .yr 1 

0.37 

330.9 

1.23 

1103.04 


5.6 Discussion: comparison between the rate of low-luminosity 
GRBs and LLA GRBs 

The rate density of low-luminosity GRBs was already considered in many studies 230 235112401 
[2371 . These bursts are selected based on their faint prompt emission, while in my sample they are 
selected based on the luminosity of their afterglow. Most of these studies rely on GRB 980425, 

GBR 031203 and GRB 060218. However, the first two were not considered in my study as they 
were not observed by the Swift satellite. 

First, a study showed that the three GRBs were taken separately, form a different sub-class 
of low-energy close bursts [2301 . There, the rate was found to be 220 Gpc Ayr 1 . As the bursts 
were connected to SNe Ib/c, a comparison with the rate of these events showed that only 0.6% 
of the type Ib/c SNe might be associated to GRBs. 

Second, an independent study (235) based on the luminosity function of GRBs (they took 
into account bursts of luminosity L < 10 49 erg.s -1 ) and on a bandwidth of [1 keV - 10 3 keV] (in 
order to compare bursts detected by different satellites) showed that the local rate of these three 
GRBs is as high as 100 — 1800 Gpc^Ayr” 1 , compatible with the result of [230: . Note that in this 
computation, the possible jet structure is not taken into account, and thus not corrected for. In 
my computation, when beaming is not taken into account, the rate density for GRB 060218 is 
equal to 50.91 Gpc _3 .yr _1 , to be compared to 380 Gpc~ 3 .yr" 1 in (2351 . 

Such large differences between the results can be explained by the different energy range 
considered in their computation [1 keV - 10 3 keV]. The flux was transformed in this band by 
extrapolating the low-energy spectrum to high energies, which might not be adequate. Indeed 
the low-energy slope is extended up to 1 MeV, while the average peak energy of the spectrum is 
on average 214 keV E3- Note that the rate of normal lGRBs was also estimated to be as high 
as 100 Gpc _3 .yr _1 12351 . which has to be compared to the rate found in a more recent study: 

1 Gpc-Ayr- 1 (ML 

Finally, in a global study of low-luminosity GRBs and of normal lGRBs I237| . in which the 
luminosity function of GRBs is explicitly taken into account (see Equation 15. 71) . the rate of low- 
luminosity GRBs is found to be in the order of 200, close to my results. However, their analysis 
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consider only the extension of the distribution function of all GRBs up to low-luminosity, with 
an imposed cut-off at 10 49 erg.s _1 , in order to take into account only low-luminous bursts. In 
other words, the distribution function of low-luminosity GRBs is not fitted for. 

To conclude, I found that the under-luminous and the LLA GRBs have comparable local rate 
densities, while being selected on different criteria (prompt versus afterglow luminosity). How¬ 
ever, this is not surprising as the rate density in my computation is dominated by GRB 060218 
and GRB 100316D. This is because these two bursts are extremely close (z < 0.1) and their 
maximum flux is below the theoretical flux detection threshold. Interestingly, if I remove these 
bursts from the computation, I found that the local rate density of LLA GRBs is a factor two 
bigger than that of normal 1GRB. 

5.6.1 Comparison with Short GRBs rate 

Another comparison is performed between LLA GRBs and short bursts. This is important as 
some LLA GRBs are compatible with the short GRBs in the E p — Ei SO plane (see section fOl of 
chapter U). 

The rate of short GRBs was estimated to be 8+3 Gpc _ 3 .yr _1 if the emission is not collimated 
in a jet, and IIOCC 470 Gpc Ayr 1 if it is (1451 . The estimation was performed using 8 short GRBs 
observed by the Swift satellite. Interestingly, 5 GRBs in my sample are at the exact position 
of short bursts: GRB 050223, GRB 051117B, GRB 060505^ GRB 120422A and GRB 031203. 
However, this last burst does not take into account the computation of the rate density LLA 
GRBs as it was not observed by the Swift satellite. Summing the rate density of these GRBs 
(without taking into account the collimation and GRB 031203) gives the rate density of these 
bursts: 0.25 Gpc~ 3 .yr~ 1 , which is much smaller (by more than one order of magnitude) than 
that of short GRBs. 

Additionally, in the computation of the rate of short bursts, Fu m was set to 1.5 counts.s” 1 , 
while in mine, it was set to Fu m = 0.4 counts.s _1 . A larger flux limit reduces the maximum 
volume Vmax and thus increases the rate. However, looking at the peak flux of the four bursts 
under consideration shows that the peak flux is much below the experimental detection threshold 
of short bursts, which means that they would have not been detected if their duration was not 
long enough. Only GRB 060505 with a peak flux of 2.65 counts.s _1 would have been detected. 
Performing the computation for this burst only gives a rate of 0.91 Gpc _ 3 .yr _1 , comparable to 
the average rate (1 Gpc _ 3 .yr _1 ) of individual short bursts. Actually, this burst was already 
thought to be a short burst of duration slightly larger than 2 s ( Tg 0 =4 s) [2411 . 

Thus, I can conclude that the LLA GRBs (except GRB 060505) which are at the position of 
short GRBs in the E p — Ei SO plane cannot be considered as short bursts with a duration larger 
than 2 s, as their rate is too small compared to that of short GRBs. Concerning GRB 060505, 
its rate density is in agreement with a short burst of duration longer than 2 s. 


5.7 Discussion for possible progenitor 

LLA GRBs are found to be intrinsically different compared to normal lGRBs. They are closer 
as indicated by their redshift distribution, and their prompt emission is much less energetic. 

However, the rate density does not allow to give a firm conclusion. Indeed, it is clearly seen 
that LLA GRBs are separated in two groups: on the one hand GRB 060218 and GRB 100316D 
(to which GRB 980425 and GRB 031203 can be added ( 239) 1. whose rate density is much larger 
than the one of normal lGRBs, and on the other hands the other LLA GRBs whose rate density 
is similar to that of lGRBs. If it is assumed that this distinction is not relevant, then it is 
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possible to conclude that the rate density of LLA GRBs is much larger than that of normal 
lGRBs. Conversely, if such distinction is relevant and is an intrinsic property of these bursts, 
then no conclusion can be obtained. 

Finally, it has to be noted that no clear difference between LLA GRBs and normal lGRBs 
was found when studying their host galaxy. However, such study is delicate, as it is difficult to 
obtain optical measurements of the host galaxies of GRBs and reliable constraints on their mass, 
metallicity and star formation rate. 

Since clear differences can be found between LLA GRBs and normal lGRBs for some of their 
properties, two situations can be discussed. First, the progenitor of LLA GRB is different from 
that of normal 1GR.B. In this case, the differences have to be explained either by considering 
different emission mechanisms or properties of the outflow. Conversely, the progenitor of LLA 
GRBs can be different than that of normal lGRBs. 

5.7.1 Same progenitor 

Instrumental effect 

First, if I consider that the differences between LLA GRBs and normal lGRBs are only due 
to instrumental effects, I can say that the progenitor is assumed to be the same, as well as 
the emission mechanisms. It implies that the LLA GRBs are the low-luminosity part of the 
distribution of GRB luminosity. 

First, it is found that LLA GRBs are observed at closer distances than normal lGRBs. This 
is not really surprising as their intrinsic faint luminosity does not allow LLA GRBs to be above 
the threshold detection of the X-ray detectors. Additionally, it is argued that when comparing all 
bursts at redshift smaller than one, the cumulative distribution functions are different. However, 
such difference can be explained by the fact that only bursts with known redshift are taken into 
consideration in this study. As the redshift is taken in the late afterglow phase of GRBs (after 
several hours), the probability of obtaining the redshift is biased towards low-luminosity sources. 
Indeed, the LLA GRBs are by definition low luminous during the afterglow phase. It means that 
the probability of successfully obtaining the redshift for a burst at intermediate redshift (a few 
tens) is biased towards intermediate redshift. 

Additionally, the fact that some bursts in the sample do not fulfil the Amati correlation can 
be explained by the fact that the X-ray instruments are not designed to determine the prompt 
properties of very faint sources. 

However, it is difficult to explain the rate density of LLA GRBs (at least for GRB 980425, 
GRB 031203, GRB 060218 and GRB 100316D) only with instrumental effects. Indeed, if the 
LLA GRBs are only the low-energy extremity of the distribution function of all GRBs, it is 
reasonable to expect that the rate density as a function of luminosity be continuous. This means 
that the rate density should not increase towards low luminosity. However, this is not the case 
when considering bursts of intermediate luminosity. 

To conclude, it is difficult to explain all differences between LLA GRBs and normal lGRBs 
by instrumental effects only. 

Environment of bursts 

The burst environment can strongly impact the properties of a GRB. LLA GRBs, being weak 
events, might even be strongly affected. Firstly, I showed that the distributions of gas and dust in 
the host galaxy and in the Milky Way for LLA GRBs are not different from their counterparts for 
normal lGRBs. Thus, absorption and extinction do not favor the same progenitor assumption. 
Only one of the LLA GRBs (GRB 090417B) is a dark burst, which means that the soft X-rays 
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are strongly absorbed: its hydrogen column density is the second largest of the whole sample. 
However, one burst cannot explain the statistical differences. 

Secondly, by using the closure relation, constraints on the surrounding burst environment can 
be given. I showed that some LLA GRBs are compatible with a wind environment (for instance 
GRB 011121). However, no clear differences were found with normal lGRBs: most bursts are 
explained by a fireball expanding in the constant interstellar medium, while a few expand in 
a wind-type environment. Interestingly, it was found that two LLA GRBs are jetted, with an 
opening angle compatible to that of normal lGRBs. 

Lastly, the host galaxies properties (type, metallicity, mass and SFR) of LLA GRBs are 
not different compared to that of normal lGRBs. To conclude, it is unlikely that all differ¬ 
ences between LLA GRBs and normal lGRBs could be explained only by the properties of the 
environment. 

Geometry of bursts 

As attested by the identification of two jet breaks in the LLA GRBs, it is expected that the 
emission in GRB be collimated in a jet, mainly in order to reduce the energy budget. Such 
non-spherical geometry is displayed in Figure 1531 If the observer is in the emission cone (region 
labeled I in the figure), the effects of the jet are noticeable only during the late afterglow, and 
they manifest themselves by the jet break. However, if the observer is outside of the emission 
cone (region labeled III), the GRB is not observable. But the SN is visible as its explosion is 
almost spherically symmetric. 



Figure 5.6: Artistic view of a bright GRB occurring in a star-forming region. The possible 
positions of the observer: on the axis on which the jet beamed (I), on jet edge (II) and off-axis 
(III). Image Credit: NASA/GSFC. 


In the unlikely situation in which the observer is on the edge of the jet (region labeled II), 
the resulting observed GRB is expected to be weak, as the Lorentz factor is smaller at the edge 
than in the emission cone. However in this case, no clear jet break is expected, which is in 
contradiction with the two jet breaks observed. Additionally, the transition between the jet and 
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the interstellar medium is expected to be sharp. Thus, the rate density of events seen on the 
edge of the jet is expected to be smaller than that of events seen on-axis. Interestingly, I found 
that the rate of LLA GRBs is larger than that of normal lGRBs, which shows that LLA GRBs 
are not even seen on the edge of a jet. 

However, it is possible to consider stratified jets [242]: in this model, the Lorentz factor of 
the burst is assumed to be a smooth function of the angle between the direction of the observer 
the central axis of the jet. In this model the transition from the central part of the jet to the 
ISM is smooth and the central part of the jet is often assumed to be narrow ( 6j ~ 1/T). 

Interestingly in this model, a jet seen at angle larger than the jet opening angle is less 
luminous, as the Lorentz factor is smaller. Its peak energy is also expected to be lower, as the 
Lorentz boost is smaller. And finally, the rate density of GRBs seen at angles larger than 9j is 
expected to be larger than that of GRBs seen on axis, since the jets are often assumed to be 
narrow. However, it is expected that the rate density should be a smoothly varying function of 
the luminosity, which is not the case, thus challenging the assumption that LLA GRBs are seen 
off-axis from a stratified jet. 

5.7.2 Content of Ejecta 

In the classical fireball model, the dynamics is mediated by the dimensionless entropy r) = 
L/(Mc 2 ). In particular, the Lorentz factor of the outflow after the acceleration phase is equal 
to r). Interestingly, there is a correlation which links the luminosity to the Lorentz factor of 
the outflow To ~ 249 L? s 3 52 j243 l. As an example, the coasting Lorentz factor of GRB 980425 
deduced from the correlation is as low as rj = To = 7.14. 

Additionally, there are many correlations between the afterglow luminosity, its kinetic energy, 
prompt luminosity and Lorentz factor. Assuming that these correlations are also valid for LLA 
GRBs (as their progenitor is assumed to be the same as the one of normal lGRBs), it is possible 
to show that the luminosity of the afterglow correlates to the Lorentz factor, so to the baryon 
load. 

Indeed, the Lorentz factor correlates to the isotropic luminosity 12431 : 

r 0 — 249L?'q 52 , (5.15) 

which correlates to the X-ray luminosity of the afterglow at 1 day 15351 : 

i.Y.lday OC L° so 74 . (5.16) 

Combining these two equations leads to: 

To oc i^ 3 lday . (5.17) 

I have shown that LLA GRBs also have low-luminous prompt phase, which indicates that 
they fulfill (at least partially) the correlation given by Eauation l5.16l However, the high-baryon 
load argument does not allow to explain the sharp transition in the rate density between the 
LLA GRBs and the normal lGRBs. Actually, a smooth transition would be expected. 

5.7.3 Different progenitors 

Here, it is assumed that LLA GRBs and normal lGRBs originate from different progenitors by 
considering the result seen in the redshift distribution and E pl -E.i so plane in Figure 14.121 This 
could be the case as the rate of LLA GRBs is larger than that of normal lGRBs. It implies that 
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the progenitor of LLA GRBs be more common than that of normal lGRBs or that the explosion 
mechanism has more chance of success. 

It was shown by m that the W-R stars able to produce a GRBs should fulfill very specific 
conditions: they should be fast rotating (needed to create the accretion torus around the black- 
hole, which powers the GRB), have a massive core (energy reservoir) and have lost their hydrogen 
envelop (to allow the jet to expands easily once it has drill its way through the core of the 
progenitor). If initially the star is single, only a few percent of O and B stars might have such 
properties once they are about to collapse. However, these conditions are more easily fulfilled if 
the progenitor is in a binary system. For instance, the outer hydrogen layers can be removed in 
the binary interaction after the common envelop stage. 

In the case of LLA GRBs, the accompanying SN (when present) are of type Ib/c, which 
means that the progenitor has lost its hydrogen envelop. If it is further assumed that the LLA 
GRBs are intrinsically less energetic (the energy put in the relativistic outflow is smaller) than 
normal lGRBs, the difference might be explain by a smaller core before collapse, which means 
less matter to be accreted by the black-hole during fall-back. This implies a progenitor with a 
smaller initial mass. 

However, a lower mass progenitor would not necessarily lose all its hydrogen envelop. Such 
problem can be solved by considering that the progenitor is in a binary system. Actually, it is a 
reasonable assumption as 70 percent of massive stars are found in binary systems |245 l [246 . The 
lower-mass progenitor is also favored by the rate density of LLA GRBs that is larger. Indeed, 
by using the Salpeter mass relation, the initial mass function of stars can be fitted as: 


dN 

dm 


-2.25 


(5.18) 


Applied to a 4OM 0 progenitor for a normal 1GRB and a 15M 0 progenitor for a LLA GRB gives 
a ratio in the number of progenitors equal to 9. The ratio of progenitors can even be increased 
when considering the other requirements to have an efficient collapsar. 

However, if the progenitor of LLA GRBs is not in a binary system, it should be massive enough 
to lose all its hydrogen and a part of its helium envelops. Thus, the energy extracted during 
the accretion onto the black-hole and used to power the bursts is expected to be comparable in 
the case of normal lGRBs and LLA GRBs. Thus, LLA GRBs would not be intrinsically low- 
luminous, the difference being obtained by the efficiency of the radiative process. The correlation 
given by Equation 15.151 indicates that the lower the luminosity of the burst emitted during the 
prompt phase, the lower the Lorentz factor, i.e. the larger the baryon load. Assuming a constant 
energy conversion by the accretion, it means that the ratio of the core mass by that of the outer 
layers is larger in the case of LLA GRBs. Thus, it implies that the progenitor of LLA GRBs are 
initially more massive, than the progenitor of normal lGRBs. It could be a WO or WC W-R 
star. 

However, a more massive progenitor is more difficult to obtain (see the Salpeter mass function 
Equation 15.181) . and thus the rate of LLA GRBs should be smaller than that of normal lGRBs. 
This disfavors the model of a single star for the progenitor of LLA GRBs. So, the purposed 
progenitor for LLA GRBs can be a low mass CO core star in a binary system. 








Chapter 6 


Conclusion and Perspectives 


6.1 Conclusion 

My thesis deals with the classification of GRBs and their tight connection to SNe. I presented 
evidence that the LLA GRBs are a subclass of long GRBs. In the global context, the present 
association of some long GRBs to SNe indicate that (some) long GRBs originate from the core 
collapse of a massive star. This idea motivated the study of SNe in my thesis. 

By using the photometric and spectroscopic data of the type lib supernova 2004ex during 
the transition (30 days < t < 60 days) and nebular (t > 60 days past explosion) phases obtained 
by the 1.82 m Ekar, 3.5 m TNG and 2.2 m WFI telescopes, the light-curve of SN 2004ex was 
completed. By fitting the light-curve with a third order polynomial function, I found that the 
peak magnitude is —17.27, reached at 21.5 days after the explosion. 

The results were compared with the prototypes of type lib SNe: SN 1993J and SN 2008ax 
respectively. I found that the light-curve of SN 2004ex is very similar to the two prototypical 
light-curves. However, it is more similar to the one of SN 2008ax. On the other hand, its helium 
lines are still well-detected 4 weeks after the explosion. And the hydrogen (H a ) line at 6286 A 
is weaker than the one found for SN 1993J at a similar period. I confirmed that SN 2004ex is of 
type lib with the rapidly declining light-curve. 

By using the similarities to those prototypes, I assumed that the mass of the 56 Ni can be in 
the order of 0.07 — 0.15 Mq, while the mass of the ejecta is 2 — 5 Mq and its kinetic energy can be 
1 — 6 x 10 52 erg. The velocity calculation of the H a line (14000 km/s) has shown that the value 
is in the expected range for type lib SNe. 

An additional comparison was performed between type lib SN 2004ex and types lb and Ic 
SNe to find similar properties. This idea came from the fact that SNe might form a continuous 
sequence being explained mainly by the envelope mass. A similar behavior is seen in the light- 
curves of type lib and type lb SNe. 

The detailed classification of SNe and the increased number of observed GRBs motivated 
the research of a more precise classification of GRBs. Indeed, GRBs are poorly classified into 
two groups (short and long bursts) based on their duration. Thus, starting from the separated 
groups of long GRBs based on the X-ray afterglow luminosity El. 1 increased the number of 
faint sources by presenting and discussing the selection method of LLA GRBs. To prepare a list 
of bursts with known redshift observed before February 15 th , 2013, I removed the short bursts 
from the sample, considering Xgo < 2 s in general. This new list consists of 254 GRBs. Then 
the different corrections (energy and cosmological) were presented in details. In addition, I have 
precisely defined how only the late afterglow was selected from the light-curves of GRBs. And 
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finally, I showed how the LLA GRBs were selected from the remaining sample of bursts, by 
considering the faintness of their afterglow. The final sample consists of the 12% faintest X-ray 
afterglows from the total population of long GRBs with known redsliift. In addition, all LLA 
GRB afterglows were re-analysed manually. Their spectral and temporal characteristics were 
determined precisely and these results are presented in Table 14.11 

By making the statistical analysis of LLA GRBs both for their prompt and afterglow prop¬ 
erties, I found that they are on average closer and less energetic than normal lGRBs, possibly 
some evidences for a different progenitor. I also showed that this sample is not more affected 
from the dust and gas compare to normal lGRBs. 

I studied the environment of the progenitor using the closure relations both in X-ray and 
optical bands. Most of them can be described by a fireball in the slow-cooling regime expanding 
in the constant ISM. Some sources ( e.g. GRB 011121) can be described by a fireball in the 
slow-cooling regime expanding in a stellar wind environment which is expected in the case of 
W-R star progenitor. The few outliers (GRB 060614, GRB 120729A) can be accounted for by 
an early jet break. These two sources show achromatic jet breaks (seen in X-ray, optical, radio 
bands). This achromatic break is interpreted as the result of a jet. The time of the break can 
be used to measure the opening angle of the jet, which give information about the efficiency of 
a matter collimation. 

I showed evidence with the Amati correlation that LLA GRBs are also intrinsically fainter 
during their prompt phase. Actually, some events do not follow at all the -E pea k - Also correlation. 
This can be explained either by the viewing angle or the energy budget of the event. I thus 
concluded that LLA GRBs are intrinsically different from normal lGRBs. 

Additional properties of LLA GRBs were also considered. I firstly noted that LLA GRBs 
include many of the supernovae associated with GRBs (64%), including the best known asso¬ 
ciations, GRB 980425, GRB 031203, GRB 060218. This means that the conclusion drawn on 
the general GRB-SN association is based on a sub-sample of the LLA GRBs which might not 
be representative of such association. I stressed the need to confirm this point and the previous 
work on GRB-SN associations using different spectral and light-curve templates. I further gave 
reasons to explain why some LLA GRBs are not associated to SNe pointing the brightness of 
their host or the large redshift z > 0.6 making the observation difficult. 

Secondly, the host galaxy properties of LLA GRBs are discussed. In the general contents 
of GRBs, the host properties are an additional argument to differentiate between long bursts 
(whose hosts are associated to young, low metallicity and actively star forming galaxies) and 
short bursts (whose hosts are associated to elliptical, low or no star forming rate, metal poor 
galaxies). Here, in the case of LLA GRBs, I found that the difference with the host properties 
of normal lGRBs are small. However, since optical constraints on the host galaxy of GRBs are 
difficult to obtain, this reduces the sample to only very few LLA GRBs. With this result, GRBs 
associated to the broad line Ic SNe show lower metallicity than the other type of SN association 
[i.e. GRB 011121) and without SNe associations (i.e. GRB 050826). Thus, more observations of 
the host would give tighter constraints and allow to conclude on the differences (if any) between 
the host galaxies of LLA GRBs and normal lGRBs. 

Finally, I computed the rate density of LLA GRBs and it is found to be much larger than 
that of normal lGRBs. However, in the sample some of events have a very high rate density. 
Even removing them I found that the rate of LLA GRBs is still two times larger than that of 
long GRBs. This can be used as an argument to differentiate these two populations. 

As a conclusion, I briefly discussed the possible progenitors of LLA GRBs, trying to explain 
the observational differences either by environmental or instrumental effects (which can not fully 
explain them) or by an intrinsically different progenitor. By using the rate of LLA GRBs, I 
indicated that the initial mass function of stars points towards lower mass progenitor for LLA 
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GRBs. However, as their accompanying SN is of type I, they should be in a binary system in 
order to have lost their hydrogen envelope during the common envelope phase. 

In general conclusion, this LLA GRBs sample is a first step towards a more precise classifi¬ 
cation of GRBs, taking into account the luminosity of the afterglow in addition to the duration. 


6.2 Perspectives 

In my thesis, I investigated if there is a distinct subclass of low-energy GRBs. However, there 
are many perspectives for my work: multi-wavelengths comparison between LLA GRBs and 
normal lGRBs, additional constraints from the host galaxies with more observations, more precise 
determination of the rate density, etc. 

Firstly, since the spectra of the afterglow in the X-ray band are peculiar, additional differences 
may also be found at other wavelengths. During next years, I will proceed on following ways by 
selecting the sources with optical data and/or radio data, then by performing a statistical study 
of the LLA GRB properties at such wavelengths and at the end I might assess the consistency 
of the results with the constraints imposed by the X-ray emission in the context of the fireball 
model: slow-cooled fireball expanding in a constant interstellar medium or wind. This can help 
to constrain the properties of progenitor if it is a very massive star which can blow off its outer 
layers with a wind or if it is a binary system. This comparisons will be more important when 
early and late follow-up multi-wavelength observations of the afterglow emission will increase 
and improve with new missions. 

I additionally showed that the prompt emission properties are different. Actually, the Amati 
correlation is based on the band function which is not physically motivated. For instance, in 
recent years, a model with three components (Band, black-body and high-energy power-law) 
was used to fit the spectra of GRBs, and one may wonder if the properties of LLA GRBs 
remain incompatible with that of normal lGRBs. Using the unique capacity of Swift which 
enables observations between the prompt and the afterglow phases, I might find some clues in 
this particular time which could allow to further constrain the LLA GRBs. Thus, not only the 
study of the afterglow but also the study of the prompt emission proved to be of high interest in 
constraining LLA GRBs. 

Secondly, I completed the computation of the rate density of LLA GRBs, finding that it 
is consistent with that of low-luminosity GRBs. However, the computation assumed that the 
rate is constant over redshifts, which is unlikely. A more precise computation can be done by 
considering that the rate density of LLA GRBs follow the star formation rate history. However, 
additional difficulties appear as the luminosity distribution function of LLA GRBs should also be 
constrained to perform this computation. It is possible to constrain the luminosity distribution, 
together with the rate density of LLA GRBs by Monte Carlo simulation, as done in previous works 
for all long GRBs [2401 . The rate density calculation is important because it can first precise 
how many events per unite time unit volume are expected, second show if there are different 
kinds of GRBs, and third be compared to the rate of other extreme events ( e.g. supernovae). 
Additionally, it can be used to constrain the population type of the progenitor 12371 - However, 
since LLA GRBs are at low redshift, their progenitors are expected to be population I stars. 

Thirdly, many LLA GRBs are associated to SNe Ib/c, including the three prototypes of the 
GRB-SNe association (GRB 980425, GRB 031203, and GRB 060218), and one may wonder if 
it is a global property of LLA GRBs. From the previous comparison, it is known that GRBs 
are more concentrated in the brightest regions of their host galaxies than are core collapse SNe, 
suggesting that GRBs are formed from the most massive stars which formerly arose from metal- 
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poor environments. Moreover, GRBs host galaxies are fainter and smaller than SNe hosts. A 
comparison between the host galaxies of the GRBs associated to SNe Ib/c and those of the SNe 
Ib/c that are not associated to GRBs could constrain the progenitors and/or the environments 
of both GRBs and SNe. Especially, metallicity measurements can give more information about 
the stellar evolution. Moreover, it is thought that most of the long GRBs make a SN, which is 
out-shined by the afterglow emission (resulting from the deceleration of the relativistic afterglow 
by the ISM). This might be an explanation for the absence of optical bright SN for some LLA 
GRBs. However, this has to be contrasted with the results of GRB 060505 and GRB 060614 
which do not show any clue for SNe associations even with deep studies. 

The common question in GRBs area is what is the origin(s)/progenitor(s) of these events. 
Is there only one kind? or more? There are several proposed candidates for the progenitor of 
GRBs. The Wolf-Rayet (W-R) star is one of the example which is strongly supported by the 
collapser model. There are three possible subclasses of W-R star: with no hydrogen envelope, 
with a small hydrogen envelope (blue supergiant star) which is one candidate for ultra long 
GRBs, and with a large hydrogen envelope (red supergiant star) which is one candidate for the 
bursts characterized by a super-soft late time X-ray radiation, a large X-ray absorption, and 
an exceptionally long prompt 7 -ray emission duration [ 1248] . In this last model, the afterglow 
radiation is reprocessed by the surrounding material which was ripped out from the progenitor 
during its evolution. This kind of study, like the one performed in my thesis, provides tight 
constraints on both the afterglow model and on the progenitor of GRBs. 

In the case of ultra-long GRBs, the blue supergiant progenitor can explain the super-soft 
late time X-ray radiation and long-prompt 7 -ray emission duration but, it can not explain the 
large X-ray absorption since it is found to be compatible with the standard afterglow model. 
Moreover, this progenitor can explain why some stars produce GRBs but not SNe. For example, 
the thermal component on the afterglow of the ultra-long GRB 111209A is too weak compared 
to the expectations of a shock breakout (for instance GRB 060218 or GRB 100316D). Adding to 
this argument the lack of SN observation challenges the W-R progenitor for this event. The long 
duration might be explained by the accretion of the remaining hydrogen layers onto the central 
compact object powering the GRB for 10 4 seconds. The presence of the hydrogen layer can be 
explained by a single low metallicity star [ 79 ] 12491 . 

In the case of LLA GRBs, I propose that low mass CO core star (WC or WO stars) in a binary 
system be their progenitor. They can lose their outer layer by interaction with the companion 
star without losing their angular momentum. The idea is supported by the fact that type Ibc 
SNe form from relatively close binaries rather than single star [250] and considering that 30% of 
LLA GRBs are associated to SNe. Indeed, if it is a single star, the CO core can be created by 
burning its hydrogen and helium layers or by ripping out the outer layers with strong rotating 
and magnetic field. In this case, there is a risk of high metallicity (e.i. Z > 0.3 Zq) which can 
remove too much angular momentum or high baryon loading in the outflow which can prevent to 
have relativistic jet 9511 . Moreover, there are example that a single massive star might produce 
GRB with high metallicity ( e.g. GRB 980425, Z > 0.7 Zq). However, this issue is still doubtful 
and open issue to be work on. 

Moreover, another three possible progenitor scenarios proposed by [252] can be considered and 
applied to have different kind of GRBs: single star which has insufficient angular momentum, the 
merger of a compact remnant with its binary companion called He merger, and the merger of two 
helium cores called He star/ He star merger which can have too much angular momentum. The 
important differences between these three possible progenitors is the distribution of redshift. 
Because the last two models can occur at large redshift while single star progenitors can be 
found at lower redshift. Indeed, the number of single star progenitor decrease with their lower 
metallicity (so increased with redshift) because they can not blow off their hydrogen envelope 
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without metal. All these possible scenarios can be considered to explain the diversity of LLA 
GRBs. 

And finally, I could try to increase the number of LLA GRBs by considering the bursts 
without redshift. First of all, I am planning to constrain the range of possible redshift for a given 
burst and to establish if it can be a LLA GRB. The method depends on the computation of the 
flux at one day in the 2 — 10 keV band, and on the mean value of the flux at one day for each 
group. Furthermore, to determine if a burst is a LLA GRB, additional points of comparison 
which do not depend upon the redshift are needed. The first additional indication is the value 
of the spectral index, which should indicate a fireball in the slow-cooling state. The second 
one may be a possible clustering between the luminosity in optical and in X-ray wavelengths, 
Loptical ~ Lx -ray, as already studied in the context of short bursts [233 . 

The study of GRBs without redshift is important to contrast the population of long GRBs. 
Indeed, as the afterglow luminosity of LLA GRBs is low by definition, it is expected that the 
efficiency of measuring their redshifts be lower than for other GRBs. As a result, many bursts 
without redshift might be LLA GRBs. Thus, determining the group of GRBs without redshift 
would increase the number of sources in the LLA sample, allowing for stronger constraints of 
their properties. 
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APPENDIX A. THE TABLE OF SHORT ORBS 


Table A.l: The list of short GRBs removed from the global sample. Most sources are presented 
in [2Tj . 


Names 

Redshift 

T90 (s) 

^90, rest (s) 

ref. 

GRB 050406 

2.44 

5.40 

1.57 


GRB 050416 

0.653 

2.5 

1.51 


GRB 050922C 

2,198 

4.5 

1.41 


GRB 051221A 

0.547 

1.4 

0.9 


GRB 060206 

4.045 

7.6 

1.51 


GRB 060801 

1.131 

0.49 

0.23 


GRB 060926 

3.208 

8 

1.9 


GRB 061006 

0.4377 

1.04 

0.72 

T. Sokomoto 2011 

GRB 070506 

2.31 

4.3 

1.3 


GRB 070714B 

0.92 

2 

1.04 

GCN 6637 

GRB 070724A 

0.457 

0.4 

0.27 


GRB 071020 

2.142 

4.2 

1.34 


GRB 071227 

0.383 

1.728 

1.25 


GRB 080520 

1.545 

2.8 

1.1 


GRB 080905A 

0.1218 

1.0 

0.47 

GCN 8187 

GRB 080913 

6.7 

8.0 

1.04 


GRB 090205 

4.65 

8.8 

1.56 


GRB 090423 

8.0 

10.3 

1.14 


GRB 090426 

2.609 

1.2 

0.33 


GRB 090429B 

9.2 

5.5 

0.5 


GRB 090510 

0903 

0.3 

0.16 

GCN 9337 

GRB 090809 

2.737 

5.4 

1.45 


GRB 090927 

1.37 

2.20 

0.93 


GRB 100117A 

0.92 

0.3 

0.16 

GCN 10338 

GRB 100206A 

0.41 

0.12 

0.09 


GRB 100316B 

1.18 

3.8 

1.74 


GRB 100724A 

1.288 

1.4 

0.61 


GRB 100816A 

0.8035 

2.9 

1.61 

GCN 11111 

GRB 101219A 

0.718 

0.6 

0.35 

GCN 14164 

GRB 130131B 

2.539 

4.3 

1.21 
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Table B.l: The global sample list. 


GRB Name 

Redshift 

ECF 

log (T a ) 

Spectral index 

GRB970228 

0.695 

1 

0 

1.0 

GRB 970508 

0.835 

1 

0 

1.4 

GRB 971214 

3.42 

1 

0 

1.1 

GRB 980425 

0.0085 

1 

0 

1.0 

GRB 980613 

1.096 

1 

0 

1 

GRB 980703 

0.966 

1 

0 

1.77 

GRB 990123 

1.60 

1 

0 

1.0 

GRB 990510 

1.619 

1 

0 

1.2 

GRB 991216 

1.02 

1 

0 

0.7 

GRB 000210 

0.846 

1 

0 

0.9 

GRB 000214 

0.47 

1 

0 

1.0 

GRB 000926 

2.066 

1 

0 

0.7 

GRB 010222 

1.477 

1 

0 

0.97 

GRB 011121 

0.36 

1 

0 

1.45 

GRB 011211 

2.14 

1 

0 

1.3 

GRB 020405 

0.69 

1 

0 

1.0 

GRB 020813 

1.25 

1 

0 

0.83 

GRB 021004 

2.33 

1 

0 

1.01 

GRB 030226 

1.98 

1 

0 

0.9 

GRB 030328 

1.52 

1 

0 

1.1 

GRB 030329 

0.168 

1 

0 

1.0 

GRB 031203 

0.105 

1 

0 

0.8 

GRB 050126 

1.29 

28.67 

2.34 

0.7 

GRB 050223 

0.5915 

18.85 

0 

1.4 

GRB 050315 

1.949 

1 

0 

0.91 

GRB 050319 

3.240 

1 

0 

0.96 

GRB 050401 

2.9 

1 

0 

1.0 

GRB 050505 

4.27 

21.51 

4.39 

1.0 

GRB 050525 

0.606 

1 

0 

1.1 

GRB 050603 

2.821 

26.61 

4.83 

0.7 
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APPENDIX B. THE TABLE OF GLOBAL SAMPLE 


Table B.l Continued 


GRB Name 

Redshift 

ECF 

Log (Ta) 

Spectral index 

GRB 050730 

3.967 

30.12 

4.13 

0.62 

GRB 050802 

1.71 

1 

3.96 

0.81 

GRB 050814 

5.3 

26.66 

3.93 

0.7 

GRB 050824 

0.83 

1 

0 

0.82 

GRB 050826 

0.297 

39.18 

0 

1.1 

GRB 050908 

3.344 

1 

0 

0.65 

GRB 051016B 

0.9364 

1 

0 

0.91 

GRB 051109A 

2.346 

26.79 

3.93 

1.0 

GRB 051109B 

0.08 

41.03 

3.67 

0.7 

GRB 051111 

1.55 

22.78 

0 

1.1 

GRB 060108 

2.7 

1 

0 

0.93 

GRB 060115 

3.53 

1 

3.86 

1.3 

GRB 060123 

0.56 

14.88 

0 

1 

GRB 060210 

3.91 

20.67 

4.46 

1.0 

GRB 060218 

0.0331 

1 

0 

0.51 

GRB 060223A 

4.41 

24.02 

2.73 

1.4 

GRB 060418 

1.489 

28.43 

3.44 

0.8 

GRB 060510B 

4.94 

1 

4.55 

1.7 

GRB 060512 

2.1 

1 

0 

0.91 

GRB 060522 

5.11 

1 

2.86 

1.1 

GRB 060526 

3.221 

35.72 

3.84 

0.7 

GRB 060604 

2.1357 

1 

4.55 

1.3 

GRB 060605 

3.78 

1 

4.16 

1.0 

GRB 060607A 

3.082 

29.25 

4.75 

0.86 

GRB 060614 

0.125 

23.11 

5.00 

0.8 

GRB 060707 

3.425 

12.98 

3.58 

1.2 

GRB 060714 

2.711 

1 

3.20 

1.4 

GRB 060729 

0.54 

11.48 

5.11 

1.39 

GRB 060814 

1.9229 

0.543el2 

4.0532469599 

1.133 

GRB 060904B 

0.703 

0.6095el2 

3.7715299824 

0.99 

GRB 060906 

3.686 

0.6112el2 

4.1314970448 

1.01 

GRB 060908 

1.8836 

0.4535el2 

3.8975334343 

1.44 

GRB 060912A 

0.937 

0.6532el2 

3.0099962952 

0.63 

GRB 060927 

5.47 

0.5562el2 

3.7187572914 

0.97 

GRB 061007 

1.261 

0.5823el2 

3.3505464558 

1.011 

GRB 061021 

0.3463 

0.5648el2 

3.6227865016 

1.02 

GRB 061110A 

0.758 

0.5641el2 

4.4328432186 

0.42 

GRB 061 HOB 

3.44 

0.6004el2 

3.5675232493 

1.19 

GRB 061121 

1.314 

0.6363el2 

4.0298207737 

0.914 

GRB 061126 

1.1588 

0.6118el2 

4.1824847413 

0.971 

GRB 061210 

0.4095 

0.4696el2 

0 

1.86 

GRB 061222A 

2.088 

0.7206el2 

4.0310314327 

0.946 

GRB 061222B 

3.355 

0.5378el2 

3.4243226011 

1.1 

GRB 070103 

2.6208 

0.405el2 

4.0593104896 

1.34 

GRB 070110 

2.352 

0.4636el2 

4.0859957897 

1.122 

GRB 070125 

1.547 

0.5669el2 

5.0697655652 

1.03 
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Table B.l - Continued 


GRB Name 

Redshift 

ECF 

Log (Ta) 

Spectral index 

GRB 070129 

2.3384 

0.4375el2 

4.2584421857 

1.32 

GRB 070208 

1.165 

0.5754el2 

4.2082998908 

1.16 

GRB 070306 

1.4959 

0.7208el2 

4.2584129662 

0.951 

GRB 070318 

0.836 

0.6172el2 

3.8568756445 

1.13 

GRB 070411 

2.954 

0.4989el2 

4.032195801 

1.22 

GRB 070419A 

0.97 

0.5485el2 

2.9713147154 

0.8 

GRB 070419B 

1.9591 

0.7011el2 

4.0618597885 

0.599 

GRB 070521 

1.35 

0.7965el2 

3.8303392967 

1.03 

GRB 070529 

2.4996 

0.7026el2 

3.0613187686 

1.06 

GRB 070611 

2.04 

0.5533el2 

4.338033476 

0.83 

GRB 070714XA 

1.58 

0.7443el2 

2.9565630412 

0.6 

GRB 070721B 

3.626 

0.6799el2 

2.9853458163 

0.62 

GRB 070802 

2.45 

0.5441el2 

4.0060330287 

1.17 

GRB 070810A 

2.17 

0.6204el2 

2.6508042691 

1.15 

GRB 071003 

1.60435 

0.6344el2 

4.534876254 

0.91 

GRB 071010A 

0.98 

0.6197el2 

4.7217273524 

1.33 

GRB 071010B 

0.947 

0.5252el2 

3.8440873614 

1.04 

GRB 071021 

2.452 

0.4835el2 

4.0443553103 

1.13 

GRB 071025 

5.2 

0.4511el2 

3.5554040155 

1.2 

GRB 071031 

2.692 

0.2161el2 

3.0605638372 

1.69 

GRB 071112C 

0.823 

0.6617el2 

3.1851790301 

0.83 

GRB 071117 

1.331 

0.6174el2 

4.0192183178 

1.09 

GRB 071122 

1.14 

0.6675el2 

2.6026019902 

0.8 

GRB 080129 

4.349 

0.7908el2 

4.0039076255 

1 

GRB 080207 

2.0858 

0.3859el2 

4.0382671125 

1.46 

GRB 080210 

2.641 

0.3777el2 

4.0733785039 

1.64 

GRB 080310 

2.42 

0.5529el2 

3.8609640866 

0.95 

GRB 080319B 

0.937 

0.5794el2 

3.6972373233 

0.82 

GRB 080319C 

1.95 

0.7391el2 

3.3446106373 

0.61 

GRB 080330 

1.51 

0.4636el2 

1.9732429172 

1.25 

GRB 080411 

1.03 

0.618el2 

4.2267829447 

0.976 

GRB 080413A 

2.433 

0.5878el2 

2.4878096887 

1.02 

GRB 080413B 

1.1 

0.5646el2 

3.4313692834 

0.966 

GRB 080430 

0.767 

0.5758el2 

4.0568712478 

1.058 

GRB 080514B 

1.8 

0.5657el2 

4.5848295892 

0.88 

GRB 080515X 

2.47 

0.6162el2 

0 

0.77 

GRB 080603A 

1.688 

0.4593el2 

4.2058511056 

1.38 

GRB 080603B 

2.69 

0.5479el2 

3.7969714195 

0.85 

GRB 080604 

1.416 

0.4719el2 

3.9040295248 

1.18 

GRB 080605 

1.6398 

0.6634el2 

4.0155375741 

0.84 

GRB 080607 

3.036 

0.5385el2 

3.0003204137 

1.13 

GRB 080707 

1.23 

0.5928el2 

3.7949566769 

1.07 

GRB 080710 

0.845 

0.5394el2 

4.3152632845 

1.02 

GRB 080721 

2.591 

0.6025el2 

4.1104208036 

0.943 

GRB 080804 

2.2045 

0.5079el2 

3.8919994881 

0.96 

GRB 080805 

1.505 

0.6098el2 

3.7601262322 

1.09 
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Table B.l Continued 


GRB Name 

Redshift 

ECF 

Log (Ta) 

Spectral index 

GRB 080810 

3.35 

0.4696el2 

3.7178179322 

1.13 

GRB 080905B 

2.374 

0.6215el2 

3.8131821575 

1.03 

GRB 080906 

2.1 

0.5513el2 

4.0923362429 

1.09 

GRB 080916A 

0.689 

0.6877el2 

4.1299475192 

1.07 

GRB 080916C 

4.35 

0.6342el2 

0 

0.96 

GRB 080928 

1.692 

0.5205el2 

3.8278494596 

1.139 

GRB 081007 

0.5295 

0.6393el2 

4.1299852479 

1.04 

GRB 081008 

1.9685 

0.5619el2 

4.0729994937 

1.04 

GRB 081028 

3.038 

0.5057el2 

4.0889469511 

1.031 

GRB 081029 

3.8497 

0.5177el2 

4.2187489895 

0.97 

GRB 081118 

2.58 

0.3785el2 

4.7580541615 

1.46 

GRB 081121 

2.512 

0.555el2 

4.0471333686 

0.944 

GRB 081203A 

2.05 

0.5148el2 

3.7887993238 

1.04 

GRB 081221 

2.26 

0.6105el2 

3.814703017 

1.29 

GRB 081222 

2.77 

0.5057el2 

3.8580866173 

1.03 

GRB 081228 

3.4 

0.5125el2 

2.4061523322 

1.15 

GRB 081230 

2.0 

0.4672el2 

4.0145379775 

1.05 

GRB 090102 

1.547 

0.658el2 

3.0665015789 

0.787 

GRB 090113 

1.7493 

0.7088el2 

2.6966321612 

1.23 

GRB 090313 

3.375 

0.5577el2 

4.9087300153 

1.08 

GRB 090323 

3.57 

0.4941el2 

0 

0.98 

GRB 090328A 

0.736 

0.7362el2 

0 

0.65 

GRB 090407 

1.4485 

0.4106el2 

4.8448295739 

1.39 

GRB 090417B 

0.345 

0.7331el2 

4.0802294785 

1.34 

GRB 090418A 

1.608 

0.6169el2 

3.0009580801 

1.016 

GRB 090424 

0.544 

0.6525el2 

3.7459348683 

0.96 

GRB 090516 

4.109 

0.5184el2 

4.0767204034 

1.114 

GRB 090519 

3.85 

0.733el2 

3.3217838938 

0.61 

GRB 090529 

2.625 

0.6211el2 

4.0894524006 

0.67 

GRB 090530 

1.3 

0.5320el2 

4.0827279454 

1.23 

GRB 090618 

0.54 

0.6417el2 

3.9244392323 

0.92 

GRB 090715B 

3 

0.5081el2 

2.6819454853 

1.05 

GRB 090726 

2.71 

0.403el2 

0 

1.45 

GRB 090812 

2.452 

0.5321el2 

2.7811412941 

1.11 

GRB 090814A 

0.696 

0.4624el2 

3.3846541723 

0.8 

GRB 090902B 

1.822 

0.581el2 

0 

0.85 

GRB 090926A 

2.1062 

0.482el2 

0 

1.09 

GRB 090926B 

1.24 

0.6527el2 

3.3139980222 

1.05 

GRB 091003 

0.8969 

0.6413el2 

0 

0.77 

GRB 091018 

0.971 

0.5217el2 

4.2666547508 

1.09 

GRB 091020 

1.71 

0.5118el2 

3.8616918296 

1.106 

GRB 091024 

1.092 

0.8777el2 

3.730564648 

0.72 

GRB 091029 

2.752 

0.2289el2 

4.0740501254 

1.752 

GRB 091109A 

3.076 

0.5619el2 

2.8168354041 

0.98 

GRB 091127 

0.49 

0.6189el2 

4.3205476385 

0.79 

GRB 091208B 

1.063 

0.6663el2 

2.7696012151 

0.94 
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Table B.l Continued 


GRB Name 

Redshift 

ECF 

Log (Ta) 

Spectral index 

GRB 100219A 

4.6667 

0.7035el2 

4.0605956206 

0.69 

GRB 100302A 

4.813 

0.5558el2 

4.0658649322 

0.87 

GRB 100316D 

0.059 

0.3586el2 

2.866936236 

0.54 

GRB 100414A 

1.368 

0.6586el2 

0 

0.63 

GRB 100418A 

0.6235 

0.7839el2 

5.1833269981 

0.87 

GRB 100425A 

1.755 

0.5639el2 

4.0192053741 

1.03 

GRB 100513A 

4.772 

0.9701el2 

3.8596097877 

1.27 

GRB 100621A 

0.542 

0.7175el2 

4.0602318088 

1.41 

GRB 100728B 

2.106 

0.5473el2 

3.6972773917 

1.09 

GRB 100814A 

1.44 

0.7192el2 

5.0099537069 

0.44 

GRB 100901A 

1.408 

0.5579el2 

4.0820697543 

1.094 

GRB 100906A 

1.727 

0.6905el2 

3.85720183 

1.026 

GRB 101219B 

0.55 

0.6678el2 

4.8061895761 

1.28 

GRB 110106B 

0.618 

0.578el2 

4.0867643399 

0.95 

GRB 110128A 

2.339 

0.5421el2 

2.8943272861 

0.84 

GRB 110205A 

2.22 

0.4697el2 

3.2664328792 

1.117 

GRB 110213A 

1.46 

0.6721el2 

3.169913591 

0.991 

GRB 110213B 

1.083 

0.3251el2 

0 

1.74 

GRB 110422A 

1.77 

0.6329el2 

3.8718316614 

0.932 

GRB 110503A 

1.613 

0.5438el2 

4.274612385 

0.943 

GRB 110715A 

0.82 

0.8073el2 

4.8743891583 

0.84 

GRB 110731A 

2.83 

0.6599el2 

3.8544184927 

0.86 

GRB 110801A 

1.858 

0.5464el2 

4.2204990554 

1.08 

GRB 110808A 

1.348 

0.3578el2 

4.3816692333 

1.62 

GRB 110818A 

3.36 

0.533el2 

3.9012167041 

1.02 

GRB 110918A 

0.982 

0.4209el2 

0 

1.34 

GRB 111008A 

4.9898 

0.5459el2 

3.8555201095 

0.91 

GRB 111107A 

2.893 

0.5213el2 

2.9246505814 

1.25 

GRB 111209A 

0.677 

0.6168el2 

4.86 

0.81 

GRB 111211A 

0.478 

0.5431el2 

0 

1.1 

GRB 111228A 

0.714 

0.5868el2 

4.0364244956 

1.024 

GRB 111229A 

1.3805 

0.6459el2 

3.8157725195 

0.86 

GRB 120119A 

1.728 

0.7543el2 

3.2088038442 

0.59 

GRB 120326A 

1.798 

0.625el2 

4.6995329121 

0.864 

GRB 120327A 

2.813 

0.7006el2 

3.8480091823 

0.85 

GRB 120404A 

2.876 

0.51el2 

3.1307370462 

1.09 

GRB 120422A 

0.283 

0.5342el2 

4.9278600406 

0.42 

GRB 120711A 

1.405 

0.7138el2 

0 

0.849 

GRB 120712A 

4.1745 

0.4122el2 

3.1173024842 

1.37 

GRB 120714B 

0.3984 

0.5039el2 

2.0733622952 

0.8 

GRB 120722A 

0.9586 

0.8533el2 

3.6963001094 

0.17 

GRB 120724A 

1.48 

0.6491el2 

4.3783121399 

0.7 

GRB 120729A 

0.8 

0.652el2 

3.8010512275 

0.79 

GRB 120802A 

3.796 

0.3508el2 

4.0072350438 

1.11 

GRB 120811C 

2.671 

0.4633el2 

3.1677811811 

1.24 

GRB 120815A 

2.358 

0.6605el2 

4.3429908712 

0.79 
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APPENDIX B. THE TABLE OF GLOBAL SAMPLE 


Table B.l Continued 


GRB Name 

Redshift 

ECF 

Log (Ta) 

Spectral index 

GRB 120907A 

0.97 

0.6342el2 

4.0950336752 

0.81 

GRB 120909A 

3.93 

0.5606el2 

3.6971794351 

1.099 

GRB 121024A 

2.298 

0.5962el2 

4.0013035671 

1.01 

GRB 121027A 

1.773 

0.4523el2 

5.141031371 

1.45 

GRB 121128A 

2.2 

0.6575el2 

3.1466561085 

1.12 

GRB 121201A 

3.385 

0.4466el2 

3.755895154 

1.17 

GRB 121211A 

1.023 

0.5648el2 

3.6240103554 

1.06 

GRB 121229A 

2.707 

0.5948el2 

4.6602363951 

1.08 

GRB 050215B 

2.62 

0.645el2 

4 

0.62 

GRB 050318 

1.44 

0.542el2 

4 

0.91 

GRB 050801 

1.38 

0.643el2 

3.77815125 

0.79 

GRB 050819 

2.5043 

0.548el2 

4.477121255 

1 

GRB 050820A 

2.612 

0.595el2 

5.812913357 

0.89 

GRB 050822 

1.434 

0.508el2 

4.230448921 

1 

GRB 050904 

6.29 

0.615el2 

5 

0.85 

GRB 050915A 

2.5273 

0.494el2 

3.966000858 

1.01 

GRB 051001 

2.4296 

0.408el2 

4 

1.2 

GRB 051006 

1.059 

0.249el2 

4 

3.4 

GRB 051117B 

0.481 

0.789el2 

4 

0.32 

GRB 060111A 

2.32 

0.531el2 

3.62797999 

0.97 

GRB 060124 

2.296 

0.658el2 

4.785329835 

0.91 

GRB 060202 

0.783 

0.358el2 

4.477121255 

1.98 

GRB 060319 

1.172 

0.546el2 

4.184691431 

1.12 

GRB 060502A 

1.51 

0.627el2 

4.421603927 

0.91 

GRB 060505 

0.089 

0.763el2 

0 

0.4 

GRB 060708 

1.92 

0.487el2 

3.906011625 

1.02 

GRB 060719 

1.532 

0.643el2 

3.909983695 

1.23 

GRB 080825B 

4.3 

0.662el2 

0 

1.17 

GRB 081109 

0.9787 

0.699el2 

4.477121255 

1.00 

GRB 100424A 

2.465 

0.782el2 

4 

0.6 

GRB 100615A 

1.398 

0.893el2 

4.278753601 

1.23 

GRB 100728A 

1.567 

0.795el2 

4.477121255 

0.73 

GRB 101225A 

0.847 

0.711el2 

5 

0.7 

GRB 111123A 

3.1516 

0.412el2 

4.608526034 

1.38 

GRB 120118B 

2.943 

0.723el2 

3.431041945 

1.04 

GRB 120521C 

6 

0.448el2 

4.301029996 

1.3 

GRB 120922A 

3.1 

0.627el2 

3.385963571 

0.91 

GRB 050408 

1.2357 

0.655el2 

4 

1.06 

GRB 051022 

0.8 

0.869el2 

4 

1.06 

GRB 070508 

0.82 

0.741el2 

4 

0.83 


End of Table 
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Data analysis process 


C.l XRT modes 

The XRT instrument can operate in auto and manual states. The auto state is the normal 
operating mode in which XRT can automatically select the science mode according to the count 
rate of the source, while the manual state is used for calibration and for a given observation in 
which the science mode can be selected. 

There are four science modes: 

• Image long and short (IM): images are taken and the source position is calculated. 

• Low rate (LR) and Piled-up photodiode (PU), which is designed for very bright sources: 
fluxes up to 60 Crab, and high time resolution: 0.14 ms. 

• Window timing (WT) mode is to obtain by binning 10 rows in serial register and readout 
only the control 200 columns of the CCD. Its time resolution is 1.7 ms and its image capa¬ 
bility is one dimensional. This mode is used when the source flux ranges 1 to 600 mCrab. 

• Photon Counting (PC) works with full image and spectroscopic resolution. In this mode 
the time resolution is limited to 2.5 seconds. It is used for the source flux below 1 mCrab. 
It can be piled-up if there are more than 2 counts per bin. This is the mode that I used to 
extract the spectral parameters of LLA GRBs. 


C.2 Directory of data 

Here, I only introduce the information about the PC mode. However, more information about 
the PC mode and other modes can be found in the Swift XRT data reduction Guide 0 and at 
the UK Swift science data centre 0. The data are mainly taken from the Swift archive download 
aortal of Leicester 0 |1411 71, and sometimes from the NASA webpage for the faintest GRBs0. 

There can be one or several observations for one GRB depending on its brightness. Each 
observation has its own ID number and it contains several directories: 

• /auxil /xrt /bat /uvot /log /tdrss 

1 blip://swift.gsfc.nasa.gov/analysis/xrt_swguide_vl_2. pdf 

2 http: / / www. swift. ac. uk/analy sis/xrt / 

3 http: / / www. swift .ac.uk/swift _portal / 

4 http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pi 
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APPENDIX C. ONE EXAMPLE OF DATA ANALYSIS PROCESS 


• /event /image /hk /productions 

Inside each observation, there is an event fits file which contains 

• the primary header, 

• the event extension: the binary table which contains the observation time, position of 
sources, pulse-height information, 

• the good time interval extension 

• the bad pixel table 

The guide gives information about all parameters, however, the event extension is in the form 
of a binary table. It is called events, which contains the observation time, position of sources, 
pulse-height information. 


C.3 Structure directory in xrtpipeline 

The processing produces two different kind of files, Level 1 and Level 2, by following the structure 
of the diagram seen in Figure 1C. II The Level 1 file is the first file which is produced without 
information loss. Additional information is calculated and added to this file during the Stage 1, 
as it can be seen in Figure ILh II In this file, data is taken in the low-rate photodiode mode during 
the slews and the settling is stored. The intermediate level of files, called Level la, is used by the 
photodiode mode and the window timing mode by taking into account the timing mode. In this 
stage, bad pixels or bad columns, coordinates information, time tagging of events, computation 
of the PHA and PI values and elimination of the pile-up frames and partially exposed pixels were 
identified. However, in stage 2, Level 2 file is created by calibration and screening. This file is 
processed following standard procedures in the stage 3. After stage 3, an event file is obtained 
and used in XSELECT. It is read and processed to extract images, light-curves, and spectra of 
events. All stages can be processed in the xrtpipeline. 


C.4 XRT file name 

XRT file names for the event and image files are built as follows: 

• sw[obs_id]x[mm] [pp]_[lew] .ext 

• sw[obs_id]x[mm]_[lew].ext 

where sw indicates the mission name (Swift), obsJd is a 11 digit number which identifies the 
observation, x represents the XRT instrument, mm is the instrument operating mode (pc, for 
more information see guide), ww identifies the window settings of the CCD (between wl and w4 
for pc mode, for more information see guide), pp identifies if the event data is taken in pointing 
mode po , during a slew si or during a settling phase sd (I used pointing mode po), the lew is the 
file level which is cl for Level 2 event file. Finally, ext represent file extension: .evt or .img. 
Moreover, the housekeeping file name is : 

• sw[obs_id]x[hh].hk 

where hh can be hd which means that information from the science-packed data is stored in the 
header. 
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Stage 1 


Stage 2 


Stage 3 


Figure C.l: Diagram of the Swift data processing in pipeline. 


C.5 Example of XRT data analysis process 

The process to obtain high-level products from level 2 screened event file is explained step by 
step below for GRB 050801. The data analysis process is follows: 

• $ heainit 

• $ cd <path of data> 

• $ cd GRB050801 

• $ Is 00148522000 00148522000.tar swxpc0tol2s0.20070101v012.rmf 

where .rmf is the XRT calibration file. 

Note: Before running xrtpipeline , the position of the sources has to be prepared. Since the 
sources considered in the thesis have known redshifts, the positions obtained from optical 
observations have been chosen with known redshifts. it is better to use the positions 
obtained from optical observations as they are more precise, even if the differences are 
usually small compared to the positions obtained by Swift. Specially, the last observation 
from Swift shows that the error radius is around 4 arcsec for XRT, while it is around 2 
arcsec UVOT. 
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APPENDIX C. ONE EXAMPLE OF DATA ANALYSIS PROCESS 


— Swift pag(0: 

1. RA (J2000): 13h 36m 35.34s (204.14724) 

2. Dec (J2000): -21 ° 55’ 42.8” (-21.92854) 

3. 90% Error radius: 1.5” 

— Optical from GCN 3736: 

RA, Dec = 13:36:33.657, -21:55:29.55 


C.5.1 XRTPIPELINE 

All files are recreated, starting from Level 1 files to Level 2 files as well as light-curves, 
spectra, etc., for each different mode if there is an observation. 

• $ xrtpipeline 


Running SWIFT XRT pipeline Task: xrtpipeline Version: 0.12.6 Release Date: 2010-07-01 


1. Source RA position (degrees or hh mm ss.s) or POINT or OBJECT[110.80875] 204.14724 

2. Source DEC position (degrees or dd mm ss.s) or POINT or OBJECT[9.50541] - 
21.92854 

3. Target Archive Directory Path[./00158593000] ./00148522000 

4. Stem for FITS input files [i.e. sw00000000000][sw00158593000] sw00148522000 

5. Directory for outputs]./00158593000_xrt] ./00148522000_xrt 


xrtpipeline_0.12.6: Exit with no errors - Sun May 25 23:52:43 CEST 2014 


C.5.2 Photon counting mode in XSELECT 

I start from Level 2 files in pc mode to tighten the range of values for parameters which 
are used in the standard screening criteria. 

• $ xselect 

1. > Enter session name >[xsel4628] pc 

2. pc:SUZAKU > read event 

3. > Enter the Event file dir >[./00148522000_xrt//] ./00148522000_xrt 

4. > Enter Event file list >[sw00148522000xwtw2po_cl.evt] sw00148522000xpcw4po_cl.evt 

5. > Reset the mission ? >[yes] 


5 http: //www. swift .ac.uk/grb_region/ 
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• $ > extract image 

Total Good Bad:Time Phase Grade Cut 
5617 5617 0 0 0 0 


Image has 5617 counts for 0.2430 counts/sec 


C.5.3 Filtering 

The filter command allows to set the different types of filters: regions, time, intensity filter, 
grade, energy, phase. It can be used with the extract command while it effects are reset by 
the command clear. 


Region filtering 

Note: To create a region filter in ds9 , I need to choose the source and background regions 
by inserting right essential (RA) and declination (Dec) of the source obtained from optical 
data. The shape of the region should be a circle for a point-like source. For the source 
region, a 30 pixel radius is considered while a 60 pixel radius is chosen for the background. 
Eventually, the central part of the source region can be removed if the count rate is above 
0.6 counts.s _1 . 

• $ > plot image 

1. RA 13:36:35.34 

2. DEC -21:55:42.8 

3. WCS 

4. sexagesimal 

5. for sources: circle(497.31858,542.88596,30) 

6. for background: circle(445.31796,400.88717,60) 

7. save “source.pc.reg” 

8. save “back_pc.reg” 


Energy filtering 

Note: This process is applied to select the lower and upper channel boundaries. The 
energy to channel conversion is PI = 100 x E where E is the energy in keV. It is used to 
select the energy from 0.5 keV to 10 keV and to remove lowest energies. 

• $ > filter pha_cut 

1. > Lower cutoff for PHA >[30] 50 

2. > Upper cutoff for PHA >[1000] 
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Grade filtering 

The grade is the shape of spread charge, which is created by an event in a the CCD. ft is 
0-12 for pc mode, ft is suggested to use full grade for the spectral analysis. 

• $ > filter grade 0-4 

• $ > extract image 

Binning of data 

Note: The number of bins is important to determine the pile-up region. When the source 
is faint, 1 used 100 counts/bin, otherwise 500 counts/bin. 

• $ > filter region source_pc.reg 

• $ > clear grade 

• $ > set bin 500 

• $ > extract curve 

Total Good Bad:Time Phase Grade Cut 
790 703 0 0 0 87 

Fits light curve has 703 counts for 3.0410E-02 counts/sec 

• $ > cpd /xw Plotting device chosen: /xw 

• $ > plot curve 

• $ PLT> q 

Time filtering 

Note: There are different methods to create time filters. Since 1 am interested with the 
late afterglow, 1 need to chose a time interval of interest by considering T a . So 1 used the 
command filter time sec , where see is the spacecraft clock time, ft needs a start and stop 
time. First, TSTART and TRIGTIME are taken from the header of the event file of the 
observation which is the good time interval file. Second, they are subtracted from T a with 
half of the time step. 

- TSTART = 1.446137710726234E+08 / time start =144613771.0726234 

- TRIGT1ME= 144613681.728 / [s] MET TRIGger Time for Automatic Target 
— Tstart - Trigtime= 89.3446234 

— Ta - Tstart - Trigtime - half of the bin =7000 - 89.3446234 - 250 =6660.6553766 

• $ > filter time see 

1. PLT> q 

2. Enter start and stop times > i 6660.6553766,60000 

3. Enter start and stop times > x 

Writing timing selections to hie cp_keybd_gti001.xsl 
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C.5.4 Extraction of light-curve 

• $ > extract curve 

Total Good Bad:Time Phase Grade Cut 
790 183 591 0 0 16 


Fits light curve has 183 counts for 9.1812E-03 counts/sec 

• $ > save curve sourc:e_pc.lc 

Wrote FITS light curve to file source.pc.lc 

C.5.5 Extraction of spectrum 

• $ > clear pha_cut 

• $ > extract spectrum 

Total Good Bad:Time Phase Grade Cut 
790 199 591 0 0 0 

Spectrum has 199 counts for 9.9839E-03 counts/sec 

• $ > save spectrum source_pc.pi 

C.5.6 Extraction of the background 

• $ > clear reg 

• $ > filter reg back.pc.reg 

• $ > extract curve 

Total Good Bad:Tinre Phase Grade Cut 
242 190 52 0 0 0 


Fits light curve has 190 counts for 9.5324E-03 counts/sec 

• $ > save curve back_pdc 

• $ > extract spectrum 

Total Good Bad:Time Phase Grade Cut 
242 190 52 0 0 0 


Spectrum has 190 counts for 9.5324E-03 counts/sec 

• $ > save spectrum back_pc.pi 

• $ > exit 

• $ > Save this session? >[yes] 

Session saved, goodbye 
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C.5.7 Creation of the ancillary response file (ARF) 

Note: The ancillary response file is called ARF and it is used to make the conversion 
between photon numbers and energy. To prepare this file, the exposure map which is 
obtained after running xrtpipeline, and the spectrum which is obtained from the manual 
extraction, are used. The exposure map is used to correct the ARF file for hot columns, 
bad pixels, any lost of counts caused by using an annular extraction region. 

$ xrtmkarf expofile=00148522000_xrt/sw00148522000xpcw4pojex.img 

1. Name of the input PHA FITS file[source_pc.pi] 

2. Apply PSF correction (used if extended=no)?(yes/no)[yes] 

3. Name of the output ARF FITS file[source_pc_exp.arf] 

4. Source X coordinate (SKY for PC and WT modes, DET for PD mode) (used if 
extended=no): [-1] 

5. Source Y coordinate (SKY for PC and WT modes, DET for PD mode)(used if 
extended=no): [-1] 


Running ’xrtmkarf_0.6.0’ 


Input Parameters List: 

— Name of the input RMF file :’CALDB’ 

— Name of the input mirror effective area file :’CALDB’ 

— Name of the input filter transmission file :’CALDB’ 

— Name of the input arf file :’CALDB’ 

— Name of the input exposure map file :’00148522000_xrt/sw00148522000xpcw4poj3x.img’ 

— Name of the input vignetting file :’CALDB’ 

— Name of the input spectrum hie :’source_pc.pi’ 

— Name of the input PSF hie :’CALDB’ 

— Name of the output ARF hie :’source_pc_exp.arf’ 

— Source SKYX Al.OOOOOO’ 

— Source SKYY Al.OOOOOO’ 

— Source off-axis angle (arcmin) :’-99.000000’ 

— Extended source? : no 


xrtmkarf_0.6.0: Exit with success. 
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C.5.8 Binning the spectrum 

Note: To bin the spectrum and to write the header I used grppha. When I wanted to use 
the x 2 ; I had to bin the data such that there were at least 20 counts per bin in the 25 
group bin. If the flux too low, I used the cash statistic which requires at least one count 
per bin. 

Furthermore, I need to use at least 15 group bins to make x 2 statistic. Otherwise, I need 
to use the cash statistic which required 1 group bin so the spectrum would have one count 
per bin. 

With the grppha command, 30 channels are removed since XRT data should not be fitted 
below 0.3 keV. ARF, RMF, and background files have been set to make the spectrum 
suitable to be used in XSPEC. 

• $ grppha 

source_pc.pi 

source_pc_15bin.pi 


1. GRPPHA [bad 0-29] 

2. GRPPHA[exit] group min 15 

3. GRPPHA[group min 15] chkey backfile back_pc.pi 

4. GRPPHA[chkey backfile back_pc.pi] chkey ancrfile source_pc_exp.arf 

5. GRPPHA[chkey ancrfile source_pc_exp.arf] chkey respfile swxpc0tol2s0_20070101v012.rmf 

6. GRPPHA[chkey respfile swxpc0tol2s0_20070101v012.rmf] exit 

** grppha 3.0.1 completed successfully 

C.5.9 XSPEC 

XSPEC is used to fit the spectrum. It is part of FTOOL§^ which is a collection of utility 
programs to create, examine or modify data files in the FITS (Flexible Image Transport 
System) format. It is distributed with the HEASOFT packag(0- 

• $ xspec 

• XSPEC12 > data source_pc_15bin.pi 
1 spectrum in use 

— Spectral Data File: source_pc_15bin.pi Spectrum 1 

— Net count rate (cts/s) for Spectrum:l 7.635e-03 +/- 7.314e-04 (76.2 

— Assigned to Data Group 1 and Plot Group 1 

— Noticed Channels: 1-382 

— Telescope: SWIFT Instrument: XRT Channel Type: PI 


6 http: //heasarc.gsfc.nasa.gov/docs/software/ftools/ftools_menu.html 

7 http: //heasarc.gsfc.nasa.gov/heasoft/ 
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— Exposure Time: 1.985e+04 sec 
— Using fit statistic: chi 
— Using Background File back_pc.pi 
— Background Exposure Time: 1.985e+04 sec 

— Using Response (RMF) File swxpc0tol2s0_20070101v012.rmf for Source 1 
— Using Auxiliary Response (ARF) File source_pc_exp.arf 

• XSPEC12>cpd /xw 

• XSPEC12>setplot energy 

• XSPEC12>ignore bad 

ignore: 370 channels ignored from source number 1 

• XSPEC12>ignore **-0.3 10.0-** 

30 channels (1-30) ignored in spectrum # 1 
25 channels (358-382) ignored in spectrum # 1 

Galactic Nh calculation 

Note: I wanted to fix the Galactic Nh value so I used XSPEC to find the value from the 
following catalog: 

• XSPEC12>nh 
Equinox (d/f 2000) [2000] 

1. RA in hh nun ss.s or degrees[110.80875] 204.14724 

2. DEC in dd mm ss.s or degrees[9.50541] -21.92854 

> > Leiden/Argentine/Bonn (LAB) Survey of Galactic HI 
— LAB > > Average nH (cm**-2) 6.42E+20 —> 0.642 —>0.64 
— LAB > > Weighted average nH (cm**-2) 6.34E+20 

Note: So Nu,Gai is 0.64 and the redshift of the source is z = 0.606 

Source and background models for the spectrum 

• XSPEC12>mo wabs*zwabs*pow 

Input parameter value, delta, min, bot, top, and max values for ... 

1 0.001( 0.01) 0 0 100000 le+06 

— 1 :wabs:nH>0.64,-1 

1 0.001( 0.01) 0 0 100000 le+06 

— 2 :zwabs:nH> 

0 -0.01( 0.01) -0.999 -0.999 10 10 
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— 3 :zwabs:Redshift>0.606 
1 0.01( 0.01) -3 -2 9 10 

— 4 :powerlaw:PhoIndex> 

1 0.01( 0.01) 0 0 le+24 le+24 

— 5 :powerlaw:norm> 

Fitting process 

• XSPEC12>fit 

• XSPEC12>fit 1000 0.001 

Model Model Component Parameter Unit Value 

— 1 1 wabs nH 10 22 0.640000 frozen 

— 2 2 zwabs nH 10 22 2.29 9 06E-09 +/- -1.00000 

— 3 2 zwabs Redshift 0.606000 frozen 

— 4 3 powerlaw Pholndex 2.54628 +/- 0.268456 

— 5 3 powerlaw norm 1.69195E-04 +/- 3.24698E-05 


Chi-Squared = 26.630 using 12 PHA bins. 

Reduced chi-squared = 2.9589 for 9 degrees of freedom 
Null hypothesis probability = 1.610149e-03 

• XSPEC12>error 2 4 

Cannot do error calc: Reduced Chi 1 2 (= 2.95886) > maximum (2) 

X 2 statistic vs. cash statistic 

Note: Since the reduced chi-squared is too large, it should be 0.8 < x 2 < 1-1 as seen 
in Table El I had to change the number of bins which should contain at least 15 to 20 
counts each. Whenever not possible, I applied the cash statistic which requires only one 
count per bin. The disadvantage of this statistic is that it does not assess the quality of the 
fit. When the first observation was not enough to get statistical significance, I combined 
the other observations until a good statistic firstly with x 2 or even with cash statistic was 
obtained. 

• $ grppha 

source_pc.pi 

source_pc_cash.pi 


1. GRPPHA [bad 0-29] 

2. GRPPHA[exit] group min 1 






140 


APPENDIX C. ONE EXAMPLE OF DATA ANALYSIS PROCESS 


Percentage Points of the Chi-Square Distribution 


Degrees of 
Freedom 




Probability of a larger value of x 2 




0.99 

0.95 

0.90 

0.75 

0.50 

0.25 

0.10 

0.05 

0.01 

1 

0.000 

0.004 

0.016 

0.102 

0.455 

1.32 

2.71 

3.84 

6.63 

2 

0.020 

0.103 

0.211 

0.575 

1.386 

2.77 

4.61 

5.99 

9.21 

3 

0.115 

0.352 

0.584 

1.212 

2.366 

4.11 

6.25 

7.81 

11.34 

4 

0.297 

0.711 

1.064 

1.923 

3.357 

5.39 

7.78 

9.49 

13.28 

5 

0.554 

1.145 

1.610 

2.675 

4.351 

6.63 

9.24 

11.07 

15.09 

6 

0.872 

1.635 

2.204 

3.455 

5.348 

7.84 

10.64 

12.59 

16.81 

7 

1.239 

2.167 

2.833 

4.255 

6.346 

9.04 

12.02 

14.07 

18.48 

8 

1.647 

2.733 

3.490 

5.071 

7.344 

10.22 

13.36 

15.51 

20.09 

9 

2.088 

3.325 

4.168 

5.899 

8.343 

11.39 

14.68 

16.92 

21.67 

10 

2.558 

3.940 

4.865 

6.737 

9.342 

12.55 

15.99 

18.31 

23.21 

11 

3.053 

4.575 

5.578 

7.584 

10.341 

13.70 

17.28 

19.68 

24.72 

12 

3.571 

5.226 

6.304 

8.438 

11.340 

14.85 

18.55 

21.03 

26.22 

13 

4.107 

5.892 

7.042 

9.299 

12.340 

15.98 

19.81 

22.36 

27.69 

14 

4.660 

6.571 

7.790 

10.165 

13.339 

17.12 

21.06 

23.68 

29.14 

15 

5.229 

7.261 

8.547 

11.037 

14.339 

18.25 

22.31 

25.00 

30.58 

16 

5.812 

7.962 

9.312 

11.912 

15.338 

19.37 

23.54 

26.30 

32.00 

17 

6.408 

8.672 

10.085 

12.792 

16.338 

20.49 

24.77 

27.59 

33.41 

18 

7.015 

9.390 

10.865 

13.675 

17.338 

21.60 

25.99 

28.87 

34.80 

19 

7.633 

10.117 

11.651 

14.562 

18.338 

22.72 

27.20 

30.14 

36.19 

20 

8.260 

10.851 

12.443 

15.452 

19.337 

23.83 

28.41 

31.41 

37.57 

22 

9.542 

12.338 

14.041 

17.240 

21.337 

26.04 

30.81 

33.92 

40.29 

24 

10.856 

13.848 

15.659 

19.037 

23.337 

28.24 

33.20 

36.42 

42.98 

26 

12.198 

15.379 

17.292 

20.843 

25.336 

30.43 

35.56 

38.89 

45.64 

28 

13.565 

16.928 

18.939 

22.657 

27.336 

32.62 

37.92 

41.34 

48.28 

30 

14.953 

18.493 

20.599 

24.478 

29.336 

34.80 

40.26 

43.77 

50.89 

40 

22.164 

26.509 

29.051 

33.660 

39.335 

45.62 

51.80 

55.76 

63.69 

50 

27.707 

34.764 

37.689 

42.942 

49.335 

56.33 

63.17 

67.50 

76.15 

60 

37.485 

43.188 

46.459 

52.294 

59.335 

66.98 

74.40 

79.08 

88.38 


Figure C.2: The % 2 distribution. Finding the probability value for a chi-square of 2.9589 with 9 
degrees of freedom. First read down column 1 to find the row for 9 degrees of freedom row and 
then go to the right close to the value 2.9589. This corresponds to a probability of less than 0.99 
but greater than 0.95. 


3. GRPPHA[group min 15] chkey backfile back_pc.pi 

4. GRPPHA[chkey backfile back_pc.pi] chkey ancrfile source_pc_exp.arf 

5. GRPPHA[chkey ancrfile source_pc_exp.arf] chkey respfile swxpc0tol2s0_20070101v012.rmf 

6. GRPPHA[chkey respfile swxpc0tol2s0_20070101v012.rmf] exit 

** grppha 3.0.1 completed successfully 

• $ xspec 

• XSPEC12 > data source_pc_cash.pi 

Note: The processes are the same as above. The only difference is the goal of the cash 
statistic. 

Cash statistic 

• XSPEC12>stati cstat 

Default fit statistic is set to: C-Statistic 

• XSPEC12>fit 1000 0.001 


Model Model Component Parameter Unit Value 
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— 11 wabs nH 10 22 0.640000 frozen 

-22 zwabs nH 10 22 1.7 1 028E-07 +/- 6.11886E-02 

— 32 zwabs Redshift 0.606000 frozen 

— 43 powerlaw Pholndex 2.61753 +/- 0.184373 

— 53 powerlaw norm 1.96417E-04 +/- 3.12063E-05 


C-Statistic = 154.967 using 157 PHA bins and 154 degrees of freedom. 

Note: Since I cannot constrain Nu,int I fixed its value to “0.0” then I fitted the spectrum 
again to find the good spectral index. This Nn^int value is also compared with the value 
from the Swift web-page H : 

— PC mode. Mean photon arrival: TO+15995 s 
— NH (intrinsic): 0 (+4 x 10 20 , -0) cm-2 

• XSPEC12>mo wabs*zwabs*pow 

Input parameter value, delta, min, bot, top, and max values for ... 

1 0.001( 0.01) 0 0 100000 le+06 

— l:wabs:nH>0.64,-1 

1 0.001( 0.01) 0 0 100000 le+06 

— 2:zwabs:nH>0 

0 -0.01( 0.01) -0.999 -0.999 10 10 

— 3:zwabs:Redshift>0.606 
1 0.01( 0.01) -3 -2 9 10 

— 4:powerlaw:PhoIndex> 

1 0.01( 0.01) 0 0 le+24 le+24 

— 5:powerlaw:norm> 

Fit results 

• XSPEC12>fit 1000 0.001 


Model Model Component Parameter Unit Value 

— 1 1 wabs nH 10 22 0.640000 frozen 

— 2 2 zwabs nH 10 22 0.0 +/- 0.836612 

— 3 2 zwabs Redshift 0.606000 frozen 

— 4 3 powerlaw Pholndex 2.84846 +/- 0.337235 

— 5 3 powerlaw norm 2.31398E-04 +/- 7.81863E-05 


C-Statistic = 154.411 using 157 PHA bins and 154 degrees of freedom. 
Note: C-Statistic/d.o.f = 154.967/154=0.8600462962963= 1 


8 http:// www. swift .ac.uk/xrt_spectra/00148522/ 
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Error calculation 

• XSPEC12>error 4 

Parameter Confidence Range (2.706) 

4 2.31647 3.41198 (-0.531971,0.563545) 

Flux calculation 

• XSPEC12>flux 0.3 10 

Model Flux 8.5685e-05 photons (2.8001e-13 ergs/cm 2 /s) range (0.30000 - 10.000 keV) 

• XSPEC12>flux 2.0 10.0 

Model Flux 2.9385e-05 photons (1.6579e-13 ergs/cm 2 /s) range (2.0000 - 10.000 keV) 

Note: I calculated the energy correction factor from the hardness ratio since I use a flux 
light-curve. 

Note: If I had used the count light-curve, I would have had to use ’’Net count rate ” to 
find the ECF. 

• XSPEC12>show all 

Net count rate (cts/s) for Spectrum:l 7.685e-03 +/- 7.266e-04 (77.4 

Spectral data counts: 197 

Model predicted rate: 6.41970E-03 



